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Chapter 1 


Introduction 


It is now commonly accepted that quantum chromodynamics (QCD) is the underlying theory of 
strong interactions Q. Consequently, one tries to understand all hadronic phenomena in terms of 
quark-gluon processes. To large extent, the success of such a program has been already achieved in 
the deep inelastic lepton-hadron collisions due to the asymptotic freedom of QCD, in the high 
energy regime one can apply perturbation theory. On the other hand, description of low energetic 
processes in terms of QCD is very difficult, since in this case the theory exhibits a complicated 
nonperturbative structure, responsible for such important phenomena like, e.g., conhnement or 
spontaneous breaking of chiral symmetry. The most fundamental approach to low energetic QCD 
(energies smaller or comparable to the typical hadron mass, ~ 1 GeV) is based on the lattice 
simulations |^. However, due to the problems with, e.g., finite size effects or proper inclusion of 
fermions, such calculations are still not fully satisfactory. In this situation, the important role is 
played by various effective theories which maintain the basic features of QCD but, at the same time, 
they are much easier to deal with. One of such theories is the Nambu - Jona-Lasinio (NJL) model. 
The NJL model was introduced already in the early sixties as a theory of interacting nucleons |Q|. 
Later it was reformulated in terms of quark degrees of freedom. Numerous calculations demonstrate 
the success of the model in describing hadronic data. In this report we are going to discuss the 
properties of hot compressed hadronic matter as described by this model. 

The two very important characteristics of QCD are: the nonperturbative structure of its ground 
state, characterized by a nonvanishing value of the quark condensate, and the appearance of light 
pseudoscalar particles, which are identified with the (quasi) Goldstone bosons of the spontaneously 
broken chiral symmetry. There is another feature of QCD which concerns hot (dense) systems. Due 
to the asymptotic freedom, at high temperature (density) the quarks and gluons are not confined 
to hadrons anymore but form the so-called quark-gluon plasma (QGP). Moreover, one expects that 
during such a deconhnement phase transition the chiral symmetry is additionally restored Q . The 
successful effective theory of QCD should incorporate these properties. 

An essential feature of the NJL model is its chiral invariance, spontaneously broken in the true 
ground state of the theory. This phenomenon leads directly to the appearance of the Goldstone 
bosons. The relative simplicity of the model allows for explicit calculations in this case. The model 
includes also the chiral symmetry restoration phase transition. This fact resembles us the case of 
QGD where an analogous phase transition is expected. In consequence, the NJL model seems to 
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be a very interesting tool to study the temperature and density dependence of various physical 
quantities. 


There are also other interesting properties of the NJL model, which deserve some comments. 
The addition of small (a few MeV) current quark masses to the NJL Lagrangian breaks explicitly 
chiral invariance. In this case, the latter becomes only an approximate symmetry. Moreover, the 
difference between various current quark masses leads to the mass splittings among the mesons 
(baryons) with different flavour content. We note that this type of behaviour is typical for QCD, 
where the small current quark masses play a similar role. The NJL model allows us also to 
understand the success of the constituent quark picture. Due to the strong four-fermion pointlike 
interaction, the light (current) quarks appearing in the NJL Lagrangian gain large (constituent) 
masses. In this way, we can find connections between the concepts of the spontaneous breaking of 
the chiral symmetry and the concepts of the nonrelativistic quark model. 

The shortcomings of the NJL model are: i) the pointlike character of the quark-antiquark 
interaction causes that the theory is not renormalized, ii) the NJL type of the quark interaction 
does not explain the confinement, moreover iii) there are no gluons in the model. The fact that 
the model is not renormalizable is not so important as long as we do not want to treat it as 
a fundamental theory. The lack of confinement is a more serious problem. This may lead to 
unphysical processes like, e.g., decays of heavy mesons into quark-antiquark pairs. The treatment 
of such cases in the NJL model must be done with a special care. Similarly, the fact that there are 
no gluons as dynamic degrees of freedom may lead to problems, especially, if one tries to apply the 
model at relatively high energies. Because of these disadvantages the model can be easily criticized. 
Nevertheless, in our opinion, the NJL model is, in spite of its deficiencies, a very interesting tool 
for studying phenomena related to the chiral symmetry and the chiral phase transition. 


Through the use of the NJL model one can evaluate the temperature (density) dependence of 
various physical quantities like, e.g., the quark condensate (constituent quark mass) ||, 0 , i, the 
pion decay constant , the meson masses @, I, n, [ 13 , n, n, the meson-meson and quark- 
meson coupling constants i 0 - In the framework of the model, the calculations of the critical 
temperature and the critical baryon chemical potential have been performed i 0 , and the order of 
the chiral symmetry restoration phase transition was studied [O]. Moreover, the thermodynamic 
properties of the quark-meson plasma were investigated |]T, |l^. Many results obtained from the 
NJL model (e.g., the temperature dependence of the quark-number susceptibility |17|, or of the 
meson screening masses [^]) are in good qualitative agreement with the lattice simulations of QCD. 


In contrast to the abundant studies of the equilibrium situations, the investigations of the non¬ 
equilibrium cases in the NJL model are relatively rare. Initially, a transport theory for the model 
has been derived in the framework of the Keldysh closed-time-path formalism combined with the 
effective action method [^]. Later, another derivation (in the mean-field approximation) has been 
presented [^], which leads to the kinetic equations whose form is explicitly chirally invariant and 
which include the spin dynamics. A certain class of solutions of the classical transport equations 
found in |]l^ has been studied in both [^] and The aim of a series of publications |^] was the 
evaluation of the in-medium cross sections. Such cross sections can be used in the phenomenological 
kinetic equations describing non-equilibrium phenomena in hadronic matter. In particular, the 
singular cross sections (giving the critical scattering at the chiral phase transition) were used in the 
rate equations describing the hadronization of the quark plasma [p^ ]. 
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In this report we give the review of the recent results obtained in the NJL model, describing 
the properties of hot compressed hadronic matter. The first large class of problems concerns the 
behaviour of static meson correlation functions. In particular, this includes the investigation of 
the screening of meson fields at finite temperature or density. Another wide range of problems 
presented in our report concerns the formulation of the transport theory for the NJL model and its 
applications to the description of high energy nuclear collisions. Recently, several reviews of the 
NJL model have been published |^]. The aim of the present article is not to give yet another wide 
and complete review, but rather to summarize a series of the new results. Our report is based on 
the following original articles: 


I W. Florkowski and B. L. Friman: Spatial Dependence of Meson Correlation Functions at 
High Temperature, Zeit. fiir Phys. A347 (1994) 271-276 [^] , 


II W. Florkowski and B. L. Friman: Meson Screening Masses in the Nambu - Jona-Lasinio 
Model, Acta Phys. Pol. B25 (1994) 49-71 []^|, 


III 


W. Florkowski and B. L. Friman: Oseillations of the Statie Meson Fields at Finite Baryon 
Density, Nucl. Phys. A in print [^], 


IV W. Florkowski and W. Broniowski: Melting of the Quark Condensate in the NJL Model with 
Meson Loops, Phys. Lett. B386 (1996) 62-64 |^], 

V W. Florkowski, J. Hiifner, S.P. Klevansky and L. Neise: Chirally Invariant Transport Equa¬ 
tions for Quark Matter, Ann. Phys. (NY) 245 (1996) 445-463 [^|, 

VI W. Florkowski: Large Time-Seale Fluetuations of the Quark Condensate at High Temperature, 
Phys. Rev. C50 (1994) 3069-3078 HI, 


VII 


J. Dolejsi, W. Florkowski, and J. Hiifner: Critieal Seattering at the Chiral Phase Transition 
and Low-pT Enhancement of Mesons in Ultra-Relativistie Heavy-Ion Collisions, Phys. Lett. 


B349 (1995) 18-22 IJ. 


Let us now outline the organization of this paper. Since our studies within the NJL model 
have been motivated by the physics of the relativistic heavy-ion collisions. Chapter || serves as a 
brief survey of high-energy nuclear reactions. This gives us a general background for the discussion 
of the properties of hot and dense hadronic matter. Chapter |3| is the introduction to the NJL 
model, where the basic concepts are presented and the connections between different formulations 
of the model are clarified. In Chapters H — ^ [EQUILIBRIUM ENSEMBLES] we discuss the phys¬ 
ical systems in thermodynamic equilibrium, whereas in Chapters [NON-EQUILIBRIUM 

PHENOMENA] the cases out of equilibrium are described. In order to make our report more 
self-contained, Chapters ^ and 0 are the short introductions to the imaginary-time and real-time 
formalisms. (We note that the imaginary-time formalism is usually used to describe the equilibrium 
ensembles, whereas the real-time formalism is used to study the systems out of equilibrium. We 
follow this tendency here.) More details about the structure of this article can be found in the 
Table of Contents. 
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Chapter 2 


High Energy Nuclear Collisions 


2.1 Relativistic Collisions 

In the relativistic heavy-ion collisions, the energies of the colliding nuclei are of the order of a few 
GeV per nucleon. In this energy range, the basic properties of the nuclear equation of state can 
be tested, which has important astrophysical relevance to neutron stars and supernova explosions. 
Here one encounters many interesting and well established phenomena like, e.g, collective flows or 
subthreshold production of particles. In general, during heavy-ion collisions, large systems of hot 
and dense hadronic matter are produced (central collisions of symmetric heavy ions at 1 GeV per 
nucleon are likely to yield about 3 times normal nuclear matter density). The particles inside such 
a hreball do not propagate completely freely: their Gompton wavelength may be comparable with 
their mean free path. In this situation, we expect that some of the particle properties (e.g., their 
masses, widths or coupling constants) can be changed. These in-medium modifications can lead to 
the experimentally observed phenomena. For example, the change of the p meson mass in dense 
matter can influence the measured dilepton spectrum |2^. Nowadays, one attempts to connect 
in-medium modihcations of hadron properties with the partial restoration of chiral symmetry. If 
this idea turns out to be correct, we can treat the change of hadron properties as a signature of 
the rearrangement of the QGD vacuum. 

The experiments with relativistic heavy-ions are intensely carried out in many places, for exam¬ 
ple, at the Gesellschaft fiir Schwerionenforschung in Darmstadt (GSI), at the Grand Accelerateur 
d’lons Lourds in Gaen (GANIL), at the Joint Institute for Nuclear Research in Dubna (JINR), and 
at the Lawrence Berkeley Laboratory in Berkeley (LBL). 


2.2 Ultra-Relativistic Collisions 

If the beam energies are larger (the energy per nucleon exceeds 10 GeV) we talk usually about the 
ultra-relativistic heavy-ion collisions. The ultimate goal of such collisions is the observation of the 
phase transition from hadronic matter to the quark-gluon plasma p7| , |^. The search for 

signatures of this phase transition is one of the most challenging problems of high-energy nuclear 
physics, both from the experimental and theoretical points of view. Several plasma signatures have 
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been proposed so far. They can be grouped into the following categories: photons and lepton pairs, 
strangeness and antibaryon enhancement, J/ip suppression, transverse flow and thermodynamic 
variables measuring equation of state (for a complete review see |31|). 

The first experiments with ultra-relativistic heavy-ion collisions took place at the Brookhaven 
National Laboratory (BNL) and CERN in 1986. The Alternating Gradient Synchrotron (AGS) 
accelerated beams up to ^*^Si at 14.5 GeV per nucleon. At GERN, the Super Proton Synchrotron 
(SPS) accelerated at 60 and 200 GeV per nucleon in 1986 and ^^S at 200 GeV per nucleon in 
1987. In 1990 a long-term project on heavy-ion physics was realized at CERN with several weeks 
of ^^S beams. Initially, no special heavy-ion machines were constructed for ultra-relativistic heavy 
ion collisions, but rather old existing accelerators had been upgraded. Only now, the completely 
new experiments take place at CERN with ^^®Pb beams. These are for the first time really “heavy” 
ions, providing large volumes and lifetimes of the reaction zone. The future of the field is connected 
with the construction of the colliders: RHIC at BNL (Au on Au reactions at ^/s = 200 GeV per 
nucleon), and LHC at CERN (Pb on Pb reactions at y/s = 6.3 TeV per nucleon). 

The present experimental evidence indicates that in ultra-relativistic heavy-ion collisions an 
extended and very dense system of hadronic matter is indeed formed. It differs in many aspects 
from the systems formed in elementary hadron-hadron reactions. However, it has also become 
clear that our search for QGP will not be an easy exercise. Although a number of “plasma signals” 
have been observed, at the same time the conventional theoretical models have been improved to 
such a level of the agreement with the data, that it does not require any radically new physics. 
Nevertheless, the physics of the ultra-relativistic heavy-ion collisions only now comes into its mature 
age. With Pb on Pb collisions at LHC, which offer the initial energy density 50 to 100 times larger 
than that of normal nuclear matter, the “new physics” of one kind or another should certainly 
appear. 


2.3 Theoretical Concepts 


The physics of the (ultra)relativistic heavy-ion collisions is an interdisciplinary field which combines 
different methods and ideas from the particle and nuclear physics. Since we deal with macroscopic 
objects containing a large number of particles, the methods of statistical physics and thermodynam¬ 
ics are especially useful. In fact, many estimates are done on the basis of purely thermodynamic 
considerations. Nevertheless, the systems produced in these collisions are not static. The need 
for the dynamical description involves rich applications of hydrodynamics. Eurthermore, since the 
matter produced in high energy nuclear collisions lives only for a short while, it is natural to ex¬ 
pect that its space-time evolution proceeds far away from equilibrium. Gonsequently, there exists a 
growing interest in applying and developing transport theories which are suitable for the description 
of non-equilibrium processes. Eor the relativistic heavy-ion collisions one uses transport equations 
of the BUU (Boltzmann-Uehling-Uhlenbeck) type, which evolve phase-space densities of nucleons, 
pions and other hadrons. With the ultra-relativistic collisions in mind, a quark-gluon kinetic theory 
based directly on QCD has been formulated |]3^ , ^]. However this approach, with a few exceptions 
[341, is still far from practical applications. 

To a large extent the successful description of ultra-relativistic reactions has been obtained 
on the basis of the microscopic Monte-Carlo simulations. These models are in some sense an 
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extrapolation of low energy hadron-hadron string models |^]. At energies of 200 GeV per nucleon or 
higher, the density of strings becomes high and they start to overlap. Consequently, the interaction 
and fusion of strings must be taken into account in most cases |36|. For extremely high energies, the 
concept of strings and flux tubes can be even abandoned and replaced by the picture of a parton 
cascade developing in the perturbative vacuum [^|. 


Last but not at least, the physics of the ultra-relativistic heavy-ion collisions has triggered the 
fast development of the quantum theory of fields in and out of equilibrium. Using the methods of 
held theory one can study the in-medium properties of particles. Moreover, this approach allows 
also for the formulation of the kinetic equations satished by the particle distribution functions. 


A description of the whole complexity of high-energy nuclear collisions in terms of a simple model 
is impossible. In particular, effective chiral theories cannot explain the mechanism of all processes 
happening at different energy scales. However, many important problems discussed in the context 
of high-energy nuclear collisions can be successfully analyzed within the NJL model. Let us mention 
a few of them: rearrangement of the QCD vacuum under extreme conditions, modihcation of the 
properties of the particles propagating through a hot and dense hadronic medium, structure of 
the chiral phase transition, and non-equilibrium dynamics of hadronic matter. It is an attractive 
feature of the NJL model that these problems can be addressed in a unified manner. 
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Chapter 3 


Nambu — Jona-Lasinio Model 


This Chapter is the introduction to the NJL model. We define here our notation and conven¬ 
tions. The relationship between the model and QCD is claryfied. The presentation starts with the 
discussion of the simplest, one-flavour version. Later on, more realistic cases are considered. 


3.1 One-Flavour Case 

The simplest variant of the NJL model is its one-flavour version based on the following form of the 
Lagrange density 


£ 


Ip [i ^ — m) il) + G {ipip)"^ -b {ipi'y^'ip) 


(3.1) 


Here ^p is the Dirac field describing quarks, m is the current quark mass, and G is the coupling 
constant |^. The latter has a dimension (mass)“^ and is assumed to be positive (in this case the 
forces between quarks and antiquarks are attractive). The field pj can carry additional colour 
degrees of freedom. In such a case, the notation in ( p.ip assumes summation over colour indices, 
for example, 'ip'ip = 'ip°‘ip'^, where Nc is the number of colours. 

Expression ( |3.1|) is an effective Lagrangian describing the local quark-antiquark interaction. The 
pointlike character of this coupling is usually explained by the assumption that the gluon degrees 
of freedom can be absorbed in the effective quark-antiquark coupling constant G. This type of 
reduction is considered to be a suitable approximation to QCD in the low-energy (long wavelength) 
limit. 

In the case m = 0, Lagrangian o is invariant under 11 ^ 4 ( 1 ) transformations 

ip^ Ip' = exp (^-*750 Ip- (3.2) 

^In our report we shall always work in 1-1-3 dimensions. The formulation of the theory in l-fl dimensions has 
different aspects, since in this case the model is renormalizable. A 2-dimensional massless theory is commonly called 
the Gross-Neveu model j^. 
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Consequently, we find that the axial current is conserved 


d^A^"{x) = 0, A^"{x) = 'ip{x)'y^j5'ip{x). (3.3) 

If the current quark mass does not vanish, m ^ 0, the equality ( |3.3D is not fulfilled anymore. In 
this case we find 


d^A^{x) = 2im'ip{x)'-f^ip{x). 


(3.4) 


Since m is small, Eq. (3.4) is usually called the partial conservation of the axial current (PCAC). 

(3.5) 


For arbitrary m Lagrangian (^) is invariant under Uv'(l) transformations, 

^ = exp V’, 


which leads to the conservation of the baryon current, 

d^V^^ix) = 0, E^(x) = xl;{x)j>^^l;ix). 


(3.6) 


3.2 Two-Flavour Symmetric Case 


A straightforward generalization of the one-flavour model leads to the following Lagrangian 


C = Ip {i ^ — rfi) Ip + 


Gs r 


i=0 




(3.7) 


In this case, the field ^ has an extra flavour index and Tj are the isospin Pauli matrices (with 
To = 1). The quantity m is a two by two diagonal matrix containing the current quark masses, i.e., 
m = diag(m„,mrf). 


In the chiral limit, m —> 0, Lagrangian (^) is invariant under Uyi(I) transformations defined 
by the rule (|3.2| ). This fact gives the conservation of the axial current, which we already know from 
the consideration of the one-flavour model, see Eq. (3.3). Moreover, also for m = 0, Lagrangian 
(3.7) does not change if we make SUa( 2) transformations 


= exp (^-* 75 ^^^ y, (3.8) 

where r = (Ti,r 2 ,r 3 ). This symmetry leads us to the conservation of the chiral current 

9^A^(x) = 0, A^(x) = 'ip{x)'y^-f5Tip{x). (3.9) 

Similarly to the one-flavour case, if m 7 ^ 0 instead of Eq. (|3.9D one finds 
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dfj_A^(x) = i'ip{x) {rrh + mr) 'y 5 'ilj{x) 


(3.10) 


In the case ruu = rrid = m, one finds that Lagrangian Q is invariant under SUv(2) transfor¬ 
mations 

= exp (3.11) 


which gives the conservation of the isospin current 

9^V^(x) = 0, V^(x) = 'ip{x)j^T'ip{x). 


(3.12) 


Finally, for any value of rh, one can check that the Lagrange density o is invariant under Uv"(l) 
transformations, such that the baryon current defined by Eq. ( |3.(j ) is conserved. 

Let us now make a few comments about our terminology. Both Uyi(l) and SUa( 2) symmetries 
will be later called axial or chiral ones. Similarly, and can be called axial or chiral currents. 
However, the name axial current will be usually reserved for A^, and the name chiral will refer to 
A^. In the analogous way, both baryon and isospin currents, and V^, can be shortly called the 
vector currents. Finally, since the theory based on Lagrangian ( p.7[ ) is for m = 0 invariant under 
both Uy 4 (l) and SUa(2 ) symmetry groups, we shall call this version of the model symmetric. 

Discussing the properties of our two Lagrangians, we have introduced four types of symmetries, 
and associated with them four conserved currents. It is well known that these symmetries charac¬ 
terize also the QCD Lagrangian. In this sense, the NJL model has symmetry features common with 
QCD. Symmetries connected with the unitary transformations Uy(l) and SUy(2) have a simple 
manifestation in nature. The first one leads to the baryon number conservation, whereas the second 
one is responsible for the organization of hadrons into the isospin multiplets. On the other hand, 
the axial symmetries U^(l) and SUa( 2) do not have any direct realization. In particular, they do 
not lead to any degeneracy of the particle spectra; If the Uyi(l) symmetry were naively realized 
in nature, each hadron would have an opposite parity partner. Similarly, a naive realization of the 
SU/i(2) symmetry would require that each isospin multiplet has a partner multiplet which groups 
the particles with opposite parity. Such a situation, however, does not take place. 

The problem with the SUa( 2) symmetry is resolved by the assumption that this symmetry is 
realized in the Goldstone mode, i.e., it is spontaneously broken. The phenomenon of the spon¬ 
taneous symmetry breaking is always accompanied by the appearance of the so-called Goldstone 
bosons. Hence, in the case of SUa( 2) we expect that there exist isovector pseudoscalar massless 
particles. Having in mind the smallness of the pion mass, we can state that such particles indeed 
exist. Speaking more precisely, pions should be regarded as quasi Goldstone bosons, i.e., particles 
which appear after the breaking of an approximate symmetry. 


The situation with the Ua( 1) symmetry looks different. Since we do not observe any Goldstone 
bosons which could be associated with Ua(1 ), we conclude that this symmetry cannot be sponta¬ 
neously broken. In fact, as it was shown by t’ Hooft |3S], the U/i(l) symmetry is broken due to the 
instanton effects. 


The two versions of the NJL model, which have been discussed so far, exhibit the spontaneous 
breaking of the Ua( 1) symmetry. In addition, the two-flavour symmetric model breaks sponta¬ 
neously the SUyi(2) invariance. Unfortunately, as we have just learnt, the Ua( 1) symmetry is not 
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spontaneously broken in the real world. Consequently, these two models must be still improved, 
if we want to use them for the description of real particles. In fact, this can be easily done by 


the modification of Lagrangian (3.7). One can add a term to expression (3.7), which simulates 
instanton effects and breaks explicitly the Uy 4 (l) invariance. At the same time, this term leaves the 
SUy4(2) symmetry intact (including its spontaneous breaking). This procedure will be described in 
detail in the next Section. 

Although the study of more realistic situations requires the elimination of the Uy 4 (l) symmetry, 
the two simplest forms of the NJL model can be still useful. The point is, that in the discussed 
versions of NJL model, the mechanism of the symmetry breaking is the same for both Uyi(l) and 
SU/i(2) groups. Consequently, in many problems where the flavour structure is not so important, 
we can study simplified theories. In this report, we shall investigate the properties of particles in 
hot and dense medium. In such a case, the extra complexities resulting from the exact consider¬ 
ation of the flavour content of particles are very often not so important. Therefore, many of our 
investigations will be based on simplified models defined by Lagrangians (^) and ( 0 )- 


3.3 Two-Flavour Standard Case 


The Ua( 1) breaking interaction, accounting for the axial anomaly in QCD, can be described by 
the following Lagrangian [^] 


Ca = Ga {det [V'i(l + 75)V’i] + det [V’i(l - l5)ipj] } 


(3.13) 


Here det stands for the determinant with respect to the flavour indices, and Ga is the coupling 
constant describing the strength of this interaction. (The last equality in ( 3.131) is valid only for 
the two-flavour case.) One can check that expression ( |3.15 ) is invariant under SUyi(2) transfor¬ 
mations. The easiest way of proving this property is the consideration of an infinitesimal SUyi(2) 


transformation. The formula (3.8) for very small x (denoted later by 6x) gives 




(3.14) 


'ijj'rtp' = il^TTp — 0^750(5%, = 0i75r0 -|- 00(5x. 


(3.15) 


One can notice that the first two terms on the right-hand-side (RHS) of the second line in ( p.l3| ) 
and the last two terms (also in this line) are separately SU^(2) invariant. On the other hand, 
Lagrangian (|3.13|) is not U/i(l) invariant. Thus, adding formula (|3.13|) to Lagrangian (|3.7D , we 
obtain a theory which is only SUa( 2) invariant. In the case Ga = Gs = G one finds 


C = Ip {i ^ — rh) ij) -\- G (00)^ -|- (0i75r0) 


(3.16) 
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The first Lagrangian having exactly the same form as expression ( |3.16 ) was written by Nambu 
and Jona-Lasinio (the second paper published in 1961 0). Since at that time quarks and gluons 
were not known, the field ^|J was used to describe nucleons (protons and neutrons). However, 
shortly after its formulation, the model of Nambu and Jona-Lasinio was abandoned because its 
non-fundamental character (reflected mainly in the nonrenormalizability of the theory) became 
clear. Later on the model was reinterpreted in terms of the quark degrees of freedom. 


3.4 Mean-Field Approximation 


Let us now come back to the discussion of the one-flavour model. By introducing the fields a and 
vr defined as 


a = — 


(3.17) 


we can recast Lagrangian (H into the form 


C = — m)^|: — a '0'i/’ ~ ~ 


< 7^+712 

AG ' 


(3.18) 


The equivalence of ( |3.1| ) and (3.18) becomes clear if one checks that the variation of (3.1) with 
respect to gives the same equation of motion for ip as the variation of (3.18) with respect to ip, a 
and TT, namely 


[i^ — m — (t(x) — Z757 r(x)] = 0. (3.19) 

In the mean-field approximation, the operators a and vr are replaced by their mean values, 
namely 

(T —> (T = {a) = Tv{pa), TT —> TT = (vr) = Tr(/57r). (3.20) 

This leads to the self-consistent mean-field equations: 


[i^ — m — a{x) — i75'7r(x)] 'ip{x) = 0, 


(3.21) 


— 2GTr [pip{x)ip{x)\ = cf{x), 


(3.22) 


— 2GTh:[p'ip{x)i'^^'ip{x)\ = 7r(x). 


(3.23) 


In Eqs. ( IT^ - (|^ ) p is the density operator and Tr denotes the trace over all physical states 
of the system. The particular form of p is not important for our considerations. We only note 
that in typical situations, p is either an operator projecting on the parity invariant ground state 
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or it describes parity invariant statistical ensembles. In these two cases (vr) = 0. In more general 
situations (e.g., when the system is out of equilibrium) we cannot exclude the case (vr) 7 ^ 0 

The U/i(l) transformation of the quark field ijj induces the change of the fields a and tt. Alto¬ 
gether, the fields tp, a and tt obey the following transformation rule 


ip - 

ip'= exp{-i'y 5 -)ip, 

(3.24) 

a - 

-> cj'= d cosy; — fr sin X, 

(3.25) 

TT - 

-> tt' = TTCosx + dsiny. 

(3.26) 


In the vector space spanned by the fields a and vr, the chiral symmetry manifests itself simply as 
a rotation. 

It is also important to realize that in the mean-field approximation the conservation of the axial 
current holds only for the mean value. Speaking more precisely, in the limit m = 0 one finds 

d^A^{x) = 2a{x)'ip{x)i^5'ijj{x) — 2TT{x)'ip{x)'ijj{x), (3.27) 

which after calculation of the trace gives 

d^{A>^{x)) = d^Tv{pA>^{x)) = 0. (3.28) 


^Generally speaking, in the space-time regions where (tt) 7 ^ 0 the qnark condensate is chirally rotated from its 
usual orientation in the isospin space. Such a piece of wrongly oriented condensate is called the disoriented chiral 
condensate Its production may explain rare events with a deficit or excess of neutral pions observed in cosmic 
ray experiments. 
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Part II 

EQUILIBRIUM ENSEMBLES 
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Chapter 4 


Imaginary-Time Formalism 


Chapters 5 — 9 of our report will be devoted to the study of systems in thermodynamic equilibrium. 
In order to calculate different physical quantities at finite temperature or at finite baryon chemical 
potential, we shall use the imaginary-time formalism [^, |^, ^|. In this Chapter, we give a short 


introduction to this formalism and discuss several computational rules used at finite T or /r. 


4.1 Temperature Green’s Functions 

While studying many particle systems, it is convenient to use the grand canonical ensemble. Intro¬ 
ducing the operator K = H — fj,N {H is the Hamiltonian and N is the baryon number operator) we 
can define the partition function and the statistical density operator by the following expressions 

Z = Tie-f^^ = e-f^^ ( 4 - 1 ) 

and 

p = Z-^e-^^, (4.2) 

where (3 is the inverse temperature and Q. is the thermodynamic potential [ks = !)■ 

The operator K can be interpreted as the grand canonical Hamiltonian. Thus, for each operator 
in the Schrodinger picture, Os{^), we can define its partner in the modified Heisenberg picture, 
namely 


x) = e^'^ 05 (x)e . (4-3) 

Using this prescription, we introduce the quark temperature Green’s function 

^afe(r,x;r',x') = -Tr |pfv [V’i^a('r,x)V;Kb('r',x')] | • (4.4) 

Here 'ipKa{T,x) is the quark field operator, a and b are the indices describing internal degrees of 
freedom (spinor, flavour and colour), and T^- is the r-ordering operator. Of course, an analogous 
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expression can be given for other types of fields. Therefore, from now on we shall continue our 
discussion treating simultaneously the case of fermions and bosons. 

In the typical situations {H= const, the system is uniform) x; r', x') is a function of 

the difference of its arguments, i.e., we can write x; r',x') = Qab{'^ ~ — x'). Using the 

fermion (boson) commutation relations and the cyclic property of the trace, we find that the Green’s 
function is an anti-periodic (periodic) function of time 

Qahij x) = =F Qah{T - t' X - x' 

allows us to write the following representation 

Qabir., X - X ) = T ^ X - x') 

n 

where 

rP 

Gabii^n-, X - x') = / dTe'’‘^"'"GabiT, X - x'), 

Jo 

and ojn are the so-called Matsubara frequencies: 

ujn = (2n -|- l)7rT (for fermions), u)n = 2n7rT (for bosons). (4.8) 

The temperature Green’s functions are very useful for the quantitative calculation of the ther¬ 
modynamic properties of the system. On the other hand, in order to study the frequencies and 
lifetimes of excited states at finite temperature, we should study the structure of the real-time 
Green’s functions. They are defined by the expression 


Equation (4.5 


). (4.5) 

(4.6) 

(4.7) 


iGabit,X-,t',x) = Tt{pT [tjjKait,x)GKb{t',x)]} , 
which is a straightforward generalization of the T = p = 0 case 


(4.9) 


iGab{t,X-,t',x) = (0|r [V’i^a(t,x)V>X6(t',X )] |0). 


(4.10) 


In Eqs. (^) and (|4.1Cl|) T is the t-ordering operator and the field operators are written in the 
standard Heisenberg picture, i.e., the following notation is used here 


OK{t,x) = e*^*05(x)e 


-iKt 


(4.11) 


We note that the real-time Green’s function ( 4.10D is defined as the expectation value of the operator 
product in the ground state |0), and that the grand canonical density operator p is reduced to the 
projection operator onto this state in the case T = p = 0. 

Besides the formula (|4.9|), one employs also the retarded and advanced Green’s functions 
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x; t\ x) = e{t - t')Tv [ipKait, yi),'il)Kb{t\ x')] (4.12) 

and 

^G'^(^>x;i',x') = -e{t' -t)Ti\^p[i)Ka{t,yi),'4>Kb{t', (4.13) 

where [ , ]± denotes the anticommutator for fermions and the commutator for bosons. These two 
functions have the following spectral representations 


and 


Gabi<^,(i) = 

r dJ Pab{bj',ai) 

J 2Tr oj — uj' -\- ie 

— OO 

(4.14) 

Gab{<^:^) = 

OO 

r duj' Pab{uj',q) 

J 271 UJ — uj' — ie 

(4.15) 


— OO 


In Eqs. ( 4.141) and ( 4.15|) the function pab{bj', q) is the spectral density (the same in the two cases) 
and e is an infinitesimally positive constant. 

The crucial point in the imaginary-time formalism is that the temperature Green’s function has 
the same spectral representation, i.e., we can write 


OO 

Qabij'^ni q) — J 


duj' Pab{uj',(i) 

2tT itOn — bj' ’ 


(4.16) 


where again the same spectral function /9afe(a;',q) appears. Using Eqs. ( 4.14 ) - ( 4.11 ) and the 
requirement of the correct asymptotics of the Green’s function in infinity, one can show that the 
temperature Green’s function (defined only at complex and discrete frequencies iujn) has a unique 
analytic continuation to the whole complex plane, where it coincides either with GQfe(a;,q) (for 

upper half-plane) or G„b(a;,q) (for lower half-plane) |Q. This allows us for the calculation of 
_ ^ ^ 

and G^^ from the knowledge of the temperature Green’s function Qab{^bJn■,^)■ Furthermore, 
knowing and G^{, we find Gab{bJ,(i) through the relation 


GabibJ, q) = [1 ± G^U^, q) + [1 ± G^,{u;, q). 


(4.17) 


The importance of the discussed connections between the temperature and real-time Green’s func¬ 
tions is related to the fact that the Feynman perturbation theory applies only to the temperature 
Green’s functions; at finite temperature the real-time Green’s functions do not have such a dia¬ 


grammatic expansion, unless one introduces an extended Hilbert space |42|. Thus, when studying 
excited states at finite temperature we first calculate the temperature Green’s function (in this case 
we can use all known techniques for Feynman diagrams), and only afterwards we do the analytic 
continuation to real frequencies. 
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4.2 Sums over Frequencies 


The Feynman rules for the temperature Green’s functions coincide with the T = fi = 0 rules 
with the replacement 






d^p 


where the continuous variable p^ should be replaced by discrete complex energies |^: 


and 


p^ —S’ iojn + ^^ = (2n + l)7rzr + p (for fermions) 


p^ —>■ ioJn = 2mTiT (for bosons). 


(4.18) 


(4.19) 

(4.20) 


c 

M 


P 


0 


Figure 4.1: The sum over the Matsubara frequencies can be calculated as a contour integral in the complex 
energy plane. 

Very often, it is convenient to convert the frequency sums to contour integrals. In the case of 
fermions, one can write [^] 


iT f{iujn + p) 

n 


J_ f dp°f{p° + p) 
2tt J exp(/3p‘’) + 1 ’ 


1 


In this report we consider the case when the bosonic chemical potential is zero. 


(4.21) 
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where Cm denotes the integration contour in the complex space, see Fig. [4.1]. Changing the 
integration variable, we obtain 


iT ^/{iun + p) = - 


+ioo+fi+e 

1 f dp^fip^) 


+ioo+fi-e 

1 f dp^f{p^) 


2tx J exp[(5{p^ — ^)] + 1 2tt J exp[/9(/x — p^)] + 1 

— 2 CO+/i+e 

+200 




(4.22) 


where C is another integration contour in the complex energy plane, shown in Fig. [4.2]. The 
position of C is fixed by the value of the chemical potential p. 



Figure 4.2: Integration contour in the complex energy plane used for the evaluation of the matter parts. 

There are two special cases which will be particularly interesting for us, i.e., the case of vanishing 
chemical potential and the case when the temperature is zero. For p = 0 and T ^ 0, Eq. (4.22) 
can be simplified to the form 


+200 


+200+e 


iT^niu„) = ^ I dpOf(pO)-^ I 

—200 —200+ e 

whereas in the limit T —> 0 and p ^ 0, we can use the rule 


dpO [/(pO) -p /(-+)] 
exp(/3+) + 1 


(4.23) 
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(4.24) 


+ 200 


+ J dp^f{p°) + ^ ^ dp^f{p°). 


In the case when the integrand depends additionally on the external frequency, e.g., / = f{p^, p, a;), 
the integrals on the RHS of expressions (|4.23 ) and ( 4.24 ) should be evaluated for purely imaginary 
values of tv and subsequently analytically continued to real frequencies. Similar prescriptions for the 
conversion of the frequency sum into the contour integrals can be given for bosons. We skip their 
presentation here since our later considerations will be restricted mainly to the fermionic sums. 


The meaning of the decompositions ( 4.23 ) and ( 4.24 ) is that we can separate the so-called 
vacuum part from the medium part. The vacuum part of a physical quantity does not explicitly 
depend on the occupation of phase space and reduces at T = /r = 0 to its vacuum expectation 
value. On the other hand, the medium part depends explicitly on the occupation of phase space 
and consequently vanishes in vacuum. 


4.3 Meson Correlation Functions 

Let us now consider again the special case when is the quark field operator. In this 

situation, we define the meson correlation function by the expression 

Xafe(Lx) = iTr|pf6if„(t,x)OK6(0,0)| , (4.25) 


where 


6 xa(t,x) = V'x(Lx)MaV'x(Lx) - Tr |pf'0i^(Lx)MaV'x(Lx)| . 


(4.26) 


The quantities Ma in Eq. (4.26) are products of the Dirac spinor matrices T^ and the Pauli 
isospin matrices Tj (in the case of the SU(3) flavour group the Pauli matrices should be replaced 
by the Gell-Mann matrices Xa). In writing Eq. ( 4.25 ) we have assumed that the Hamiltonian is 
independent of time and that the system is uniform (this allowed us to make a shift of one of the 
arguments to zero). 

The meson correlation function (4.25) describes the real-time propagation of a meson which con¬ 
sists of a qq pair. Similarly to the case of one-particle Green’s functions discussed before, in order 
to calculate the real-time correlation function ( 4.25| ) we initially calculate the temperature correla¬ 
tion function and subsequently do the analytic continuation. The temperature meson correlation 
function is defined as 


Xab('r,x) = Tr |/5f.rOxa('r,x)di^b(0,0)| , 


(4.27) 


where the held operators are written in the modihed Heisenberg picture (O) (we do not intro¬ 
duce the special notation for the temperature correlation function, it can be distinguished by its 
imaginary time argument r). 
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Chapter 5 


Quark Self-Energy and the Gap 
Equation 


In the NJL model, one calculates the quark self-energy from the self-consistent Schwinger-Dyson 
equation. This is usually done either in the Hartree or Hartree-Fock approximation. Furthermore, 
different regularization schemes can be used to define divergent integrals. In order to illustrate these 
possibilities, we come back to the separate discussion of the one-flavour, two-flavour symmetric, and 
two-flavour standard models. The one-flavour case will be treated in the Hartree approximation 
and the 3-dimensional cutoff will be used as a regulator. The two-flavour symmetric case will be 
considered in the Hartree-Fock approximation and we shall regularize it using a version of the 
Pauli-Villars method. Finally, we shall discuss the standard version of the two-flavour model in the 
Hartree-Fock approximation, regularizing it with the help of the Schwinger proper-time method. 

We analyze the differences in the formulation of the theory in more detail, since several problems 
of interest for us will be studied later in different versions of the model — it often happens that one 
version is more suitable for the analysis of a given problem than other versions. In particular, the 
in-medium dependence of the dynamic and screening masses of mesons (Chapters 6 and 7) as well 
as the oscillations of the static meson fields at finite baryon density (Chapter 8) are studied within 
the two-flavour symmetric model. The methods going beyond the Hartree-Fock approximation 
(Chapter 9) are introduced on the basis of the two-flavour standard model. On the other hand, the 
transport theory for the model (Chapter 11) is constructed for the one-flavour model. 

Using the methods introduced in Chapter 4, we calculate the quark self-energy at finite tem¬ 
perature or density. In this way we can determine the in-medium dependence of the constituent 
quark mass. Illustrative calculations done in different versions of the model lead to the same qual¬ 
itative conclusion; with increasing T (n) the constituent mass decreases, and for sufficiently large 
temperature (density) it approaches the value of the current mass. 


5.1 One-Flavour Case 

In the one-flavour model, using the Hartree approximation we write the self-consistent Schwinger- 
Dyson equation in the following form 
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(5.1) 


T. = 2iG J |tr [G(p)] + z 75 tr [i--(^G{p)] | , 

where tr is the trace over spinor (and colonr) indices and G{p) is the quark propagator 

G{p)~^ ='^— Tj — m + ie. (5.2) 

The two terms on the RHS of Eq. (|5.1|) appear since Lagrangian ( |3.1| ) describes both scalar and 
pseudoscalar interactions, see Fig. ||5.1| . However, parity implies that the second term in (|5.1| ) 
vanishes. 



Figure 5.1: Diagrams used for the calculation of the quark self-energy in the Hartree approximation (one- 
flavour model). The solid line describes the quark propagator, whereas the dashed line is introduced to 
describe the point-like quark-antiquark interaction. Only the first term gives a non-zero contribution to the 
self-energy. 

Using the general formula for the quark propagator, Eq. o, we define the quark condensate 
through equation 


im) = 


Rewriting Eq. (5.3) in the momentum space, 


-itiG{x = 0 ). 
we find that 


{qq) = -ANc i 


d'^p + m 

(27r)^ p'^ — {T, + mY + ie' 


(5.3) 


(5.4) 
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The sum S + m can be regarded as the constituent (effective or dynamic) quark mass M, whose 
origin are the quark-antiquark interactions. Eqs. (0) and (|5.4f ) lead to the simple relation 


M = m — 2G{qq), 


(5.5) 


which connects the constituent quark mass directly to the condensate. Eq. (|5.5|) is called the gap 
equation since it determines the energy gap in the energy spectrum of one-particle excitations. 


In writing expressions (0) - (0 we have used the Feynman rules for T = = 0. As it 

was pointed out in Chapter 4, the Feynman rules at finite T and fi have the same structure. 
Consequently, we can always start our discussion with the consideration of the standard Feynman 
diagrams. Later on, the expressions valid at finite T and ^ are obtained by employing the prescrip¬ 
tions: (4.18), (4.23) and (4.24). In particular, it turns out that the form of formula (5.5) is valid 
also at finite T and p provided one writes (qq) as a sum of two contributions 


{qq) = {qq)vac + {QQ)med- (5-6) 

For the moment, let us concentrate in more detail on the vacuum part of the condensate {qq)vac- 
Of course, it can be evaluated directly from Eq. (^.4[) . For further use, it is convenient to introduce 
the function /i,vac(Af^) defined as 


h,Yac{M‘^) — 8iNc 


ddp 1 

(27r)‘^ -b ie 


Using Eqs. (p.4|) and (|5.7|) we write 


(5.7) 


{QQ)vac — —-^h^vaciM^)- 


Integration over energy in (5.7) leads to the following expression 


(5.8) 


/i,vac(M2) = AN, I Ep = ^p^ + MT (5.9) 

The 3-dimensional integral appearing in Eq. ( |5.9| ) diverges. The simplest way of its regularization 
is the introduction of a cutoff A, confining the region of integration to the sphere |p| < A. In this 
case one finds 


jR,3D 

-'i,vac 


(M^) 



/M\2 a -b VA2 -b M2 

Uj - M - 


(5.10) 


Using now Eqs. (^), (^) and (5.1C) we obtain the following explicit form of the gap equation 


M 


1 - 


m 

M 


= M 


NrGA^ 


TT^ 



/M\2 a- b VA2 -b M2 

U) - M - 


(5.11) 
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It is interesting to look at some consequences of Eq. (^.11 ). At first, let us consider the case of 
the chiral limit m = 0. The expression in the square bracket on the RHS of formula (|5.11j ) is positive 
and smaller than 1 for M > 0. Thus, the nontrivial solution exists only if 0 < nyiNcGA^) < 1. See 
Fig. [5.2], where Eq. ( ^.11| ) has been plotted for two different values of G. As GA^ increases over 
the critical value vr^/A^cj AI starts rising from zero. Consequently, the form of Eq. ( |5.11 ) indicates 
that the force between quarks and antiquarks must be positive {G > 0) and strong enough to 
cause the formation of the condensate. In the analogy to the BCS theory of super-conductivity, 
one expects that the nontrivial solution of Eq. ( 5.11 ) corresponds to the true ground state of the 
model Q This is indeed so, and this fact can be verified by performing a variational calculation 
(for more details on this subject see the original paper by Nambu and Jona-Lasinio @| or the last 
review listed in |^5[). If the current quark mass m does not vanish, one can find out that there 
is always a nontrivial solution to Eq. ( 5.1lj ). Moreover, it corresponds to the ground state of the 
theory, similarly as in the special case m = 0. 



Figure 5.2: Graphical representation of Eq. ( 5.11 ) in the limiting case m = 0 and for tw o different values of 
the coupling constant G. The common factor M, appearing on both sides of Eq. ( ji.lll ), has been omitted. 
The LHS is represented by the solid line, whereas the RHS is represented by the two dashed lines. Only if 
GA^ is greater than the critical value there exists a nontrivial solution. 


Of course, formula (5.11) is valid only for T = fi = 0. If the system is in contact with a heat 
bath or with a reservoir of particles, Eq. (|5.11 ) has to be generalized. To hnd the in-medium form 
of the gap equation we calculate the in-medium part of the condensate (9Q')ixLed' Similarly to Eq. 


^In the Hartree approximation the true ground state of the theory can be described as consisting of the Dirac sea 
of the constituent quarks — complicated vacuum of strongly interacting light current quarks is approximated by a 
simple vacuum of the heavy constituent quarks. 
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(5.8) we can write 


(99)med 2 )• 

In the special cases (i.e., for T > Oand/U = 0 or for T = Oand/U > 0) we have: 


(5.12) 




200 + e 


— 16iA^c 


—'ico+e 


oo , 

ANc f dpp^ 

vr^ J Ep 

0 


dp^ r d^p 1 

27r J (27r)3 Pq — 


1 

e^p/T -\- 1 


1 

M2 eP°/^ + 1 


(5.13) 


and 


A,mat(^') - 


-p2-M2 


2Nc f dpp'^ '■ '^--2 

0 0 


2Nc f dpp^ 


TT^ 


En 


(5.14) 


The form of Eqs. (5.13) and (5.14) follows from the application of the replacement rules (4.18), 
( ^-231) and ( 4.24 ) in the expression on the RHS of Eq. (^]^) . 


It is important to realize that expressions ( 5.13| ) and ( ^.141) are well defined converging integrals, 
and they do not have to be regularized, which is in contrast to the calculation of the vacuum part 
(5.7). Note that the sum 7i,vac(-^^) 


+ med(-^^) should vanish in the limit T —> oo or > oo. 

This condition follows from the general property of the fermionic Matsubara sums. Since we have 
introduced the finite cutoff in the calculation of Ii^Yac{M‘^), the discussed condition may not be 
fulfilled. For instance, at finite temperature we find the expression 


jR,3D 

^i,vac 


(M^) + 





dpp'^ 2 
Ep e^p/^ + l 


(5.15) 


which does not vanish at T ^ oo. However, if the function tem(-^^) i® regularized in the same 
way as the function /i,vac(A7^), the sum desired high temperature 

limit. This example leads us to the simple conclusion: if the cutoffs used for the regularization of 
the vacuum parts are finite, the temperature (matter) parts must be additionally regularized. This 
ensures the correct high temperature (density) behaviour of the complete expression. 

With these conditions in mind, we write the finite temperature gap equation in the form 


A 


M 


1 - 


m 

M 


= M 


2NcG f dpp^ 


TT^ 


En 


1 - 


gEp/T _j_ \ 


(5.16) 
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Since the RHS of Eq. (|5.16D vanishes in the limit T ^ oo, the constituent quark mass M drops 
to the current one m (with increasing T). In the special case m = 0, there exists a critical 
temperature Tc such that for T > Tc formula ( |5.16| ) has only a trivial solution M = 0. This 
fact implies the restoration of chiral symmetry at high temperature. See Fig. [5.3], where the 
temperature dependence of the constituent quark mass has been plotted. 



Figure 5.3: Temperature dependence of the constituent quark mass M in the limiting case m = 0. If 
temperature is greater than the criti cal o ne (in this case Tc turns out to be 190 MeV) there exists only a 
trivial solution to the gap equation ( ^.16| ), i.e., M = 0 for T > Tc. This fact indicates the restoration of 
chiral symmetry. In the present calculation the 3-dimensional cutoff regularization was adopted and we used 
the following parameters: A = 0.65 GeV, G = 0.4 fm^ and iVc = 3 (this gives M = Mq = 313 MeV at 
T = 0). 

In the analogous way, for finite densities but T = 0, one finds 


M 


1 - 


m 

M 


= M 


2NcG f dpp^ 


TT^ 


Er> 


(5.17) 


PF 


where pf is the Fermi momentum of quarks, i.e., pf = Calculating explicitly the 

integral on the RHS of expression (5.171) one finds 
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M 


1 - 


m 

M 


= M ■ ^ 



fM\‘^ fM 

'nx) -(t 


2 A + VA2 + M2 
In — 


M 


\Pf) M J' 


(5.18) 


Equation ( ^.17 ), or equivalently Eq. ( |5.18| ), leads again to the chiral restoration. See Eig. [5.4], 
where the constituent quark mass (being the solution of Eq. ( ^.17 ) for the case m = 0) has been 
plotted as a function of pp- Now the phase transition takes place at finite baryon density. Similarly 
to the finite-temperature case, there is a critical value pp, such that for pf > Pp there is only a 
trivial solution to the gap equation. 


> 

CJ 

s 

s 



Figure 5.4: Density dependence of the constituent quark mass M (for the 3-dimensional cutoff regularization 
and in the chiral limit m = 0). The mass is plotted as a function of the quark Fermi momentum pp. In this 
case, there is also a chiral restoration phase transition happening at pp = 314 MeV. The numerical values 
of the parameters are the same as in the previous figure. 


To our knowledge, the first evaluation of the in-medium dependence of the constituent quark 
mass in the NJL model was given by Hatsudo and Kunihiro Q. Afterwards, many authors repeated 
this kind of calculation. The results shown in Figs. [5.3] and [5.4] are obtained for the set of 
parameters equivalent to that used in [1T|. The decrease of M with increasing temperature, see 
Fig. [5.3], implies the decrease of the quark condensate (strictly speaking of the absolute value of 
{qq)). We note that, a similar tendency in the behaviour of (qq) has been observed in the lattice 
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simulations of QCD as well as in the calculations based on the chiral perturbation theory. A more 
quantitative comparison of the NJL results with other estimates will be given in Chapter 9. 


5.2 Two-Flavour Symmetric Case 

We shall consider the two-flavour symmetric model in the Hartree-Fock approximation. The ap¬ 
propriate form of the Schwinger-Dyson equation, see Fig. [5.5], is now 


S 


l^tr [TiG{p)] - TiG{p)Ti 
+ [iri75G(p)] - ira5G{p)iTij5 \ . 


(5.19) 


Because the second term with the trace in Eq. ( ^.19 ) vanishes and the terms without the trace 
cancel each other, expression (5.19|) can be rewritten in a simpler form, namely 


^ = GsiJ2 

i=o • 


d'^p 


Tjtr [TiG{p)]. 


(5.20) 



Figure 5.5: Diagrams contributing to the quark self-energy in the Hartree-Fock approximation (two-flavour 
symmetric model). We note that the introduction of the dashed lines, describing the point-like quark- 
antiquark interaction, allows us to distinguish graphically the Hartree and the Fock terms. 
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We note that the quark propagator has now extra isospin indices (in comparison to the one- 
flavour model discussed before) and tr includes now the trace over the isospin space. Consequently, 
in the present case we have 2(gg) = —iivG{x = 0“). This is due to our convention that (^q) de¬ 


scribes the contribution to the quark condensate from one flavour only, see Eq. (5.4). Nevertheless, 
one can check that the form of Eq. ( ^) remains unchanged in the present case — the replacements 
( 99 ) 2 (^( 7 ) and G —> G 5/2 in (^) compensate each other giving M = m — 2Gs{qq)- 

The 3-dimensional regularization scheme (introduced in the previous Section) has several draw¬ 
backs, e.g., it is not Lorentz invariant. Therefore, we shall now discuss another method of regu¬ 
larization. As we shall later see, it can be regarded as a version of the Pauli-Villars subtraction 


scheme. At first, let us consider the function /i,vac(-^^) defined by Eq. (^). By making the Wick 
rotation to the variable 734 = —ip^ we can write 


h,vac{M^) = 8N, I 


d'^PE 


(27r)4 p\ + M‘^' 

where d'^pE = dp^d^p and p% = p\ + p^. We replace now /i,vac(-^^) by the series 


(5.21) 


N 

- /fvac = E ^*A,vac(A?). (5.22) 

j=0 

Here Aq = 1, Aq = M and N is the number of the so-called subtractions. The regulating masses 
Aj and the coefficients Ai [i > 0) should be chosen in such a way as to provide the finite result for 
This requirement leads to the following set of constraints 


N N N 

Y.A, = 0, 1:a,A? = 0, ... 

2=0 2 = 0 2=0 


Using Eq. (p.23| ) we find 


(5.23) 


jR,PV 

-'i,vac 


(M^) 


N ^ 
2=0 


(5.24) 


As it was discussed above, the temperature (matter) parts should be regularized as well. We 
do it again by making the replacement 


I. 


d(M^) - = 


i,mei 


N 

E ^*A,me' 

i=0 


d(A? 


(5.25) 


where med(''^i) defined by Eqs. ( 5.13|) and ( 5.14 ). 


In Eigs. [5.6] and [5.7] we show the in-medium dependence of M obtained from the two-flavour 
symmetric model regularized with the help of the Pauli-Villars method described above. Of course, 
these two figures are analogs of Figs. [5.4] and [5.5] from the previous Section. In the present case 
we do not limit ourselves to the limit m = 0. As the result we find that M does not drop suddenly 
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to zero, but it goes smoothly to the value M = m for T ^ oo or ^ oo. Moreover, the method of 
subtractions “switches off” large momenta more gradually than the 3-dimensional cutoff procedure. 
Thus, the temperature and density dependence of M shown in Figs. [5.6] and [5.7] looks weaker 
than that shown in Figs. [5.4] and [5.5]. Nevertheless, the qualitative behaviour of M in the two 
versions of the model is the same. In particular, vanishing of the condensate at large T or ^ in 
both cases indicates the chiral symmetry restoration phase transition. 



Figure 5.6: Temperature dependence of the constituent quark mass. Pauli-Villars regularization scheme 
(as described in the text) was employed here with three regulating masses: Ai = 0.68 GeV, A 2 = 2.1Ai and 
A 3 = 2 .IA 2 . The coupling constant Gs = 0.75 fm^, and the current quark mass m = 8.56 MeV. 

We note that at T = 0 all values of the chemical potential in the range 0 < /x < M are equivalent 
and correspond to the physical vacuum. The reason is that the chemical potential jj, is the energy 
needed to add a quark to the system. Clearly ^ must exceed the minimal energy, the vacuum rest 
mass of a quark M, before the Fermi sea of quarks begins to be populated. 


35 




350 


^ 250 

> 

(U 

s 

S 150 


50 


400 450 500 

11 [MeV] 


Figure 5.7: Constituent quark mass as a function of the quark chemical potential /x. Pauli-Villars regular¬ 
ization procedure was used here with the same parameters as in Fig. [5.6]. 


5.3 Two-Flavour Standard Case 

Let us turn to the discussion of the two-flavour standard model. Calculating the quark self-energy 
in the Hartree-Fock approximation one finds 


S = 2Gi 


d^p 


tr[G(p)] - G{p) + irystr [ir^i^Gip)] - iT'y 5 G{p)iT'y 5 \ , 


. (27r)4 

which after calculation of the trace over spinor indices gives 


(5.26) 


M = m + 4iMG{2NcNf + l) J 


d^p 


(27r)^ _ jif 2 _|_ jg 


(5.27) 


In contrast to the symmetric version of the model, in this case the Fock terms (see Fig. [5.8]) 
do not cancel each other. The extra “1” in ( ^.27 ) comes just from the inclusion of the exchange 
terms. Nevertheless, for Nc = 3 and Nf = 2 the contribution from the Fock terms can be regarded 
as a small correction to the result obtained only with the Hartree approximation. Consequently, 
one typically neglects the Fock terms arguing that this is allowed in the leading order of the 1/Nc 
expansion. A more rigorous treatment of the model in the 1/Nc expansion, which generalizes the 
results of this Section, will be presented in Chapter 9. 
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Neglecting the Fock terms, we find the gap equation in the form 


M = m- 4G{qq) = m + 2MGIi^yac{M^)- 


(5.28) 


In this case, we can again use our definitions ( |5.4| ) and (^.8|). We note that there is a factor AG in 
the gap equation (5.28), which is in contrast to our two previous cases where a factor 2G appeared. 



Figure 5.8: Diagrams used for the calculation of the quark self-energy in the case of the standard two-flavour 
model. 

Discussing the standard version of the model, we shall introduce yet another regularization 
procedure. This is so-called Schwinger proper-time method based on the following replacement 


/i,vac(M2) = 8W J 


d'^PE 1 

(27r)4 -h M2 


/ ds ew {-s\pI + M‘^]Y (5.29) 


(2vr) 


Introducing the exponential integral function defined as 


we can write 


En{x) = J dtt ‘^exp{—xt), 
1 




A2 


(5.30) 

(5.31) 
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Using the Schwinger proper-time regularization at finite temperature, it is not convenient to 
separate the vacuum and temperature parts. We only replace the integral over energy by the sum 
over the Matsubara frequencies. According to formula ( 4.18| ), we can generalize Eq. (5.31) to the 
finite temperature case writing 



Solving the gap equation ( 5.28 ) with h,vac{M^) replaced by T) we find again the 

characteristic T-dependence of M, similar to that shown in Fig. [5.3] or [5.6]. 


At this place we close our considerations of the Schwinger-Dyson equation in the different 
formulations of the model. The results derived by us in Sections 5.1 — 5.3 will be frequently used 
in the next Chapters. In particular, we shall return to the discussion of the two-flavour standard 
model and the Schwinger proper-time method in Chapter 9, where the calculations extending the 
HE approach will be presented. 
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Chapter 6 


Mesonic Excitations 


Mesons are quark-antiquark excitations of the true ground state of the theory. Their quantitative 
description requires introduction of the meson correlation functions, as we did it in the end of 
Chapter 4. In the NJL model, the overall scheme used to describe mesons is similar to the standard 
HF (Hartree-Fock) -|- RPA (random phase approximation) approach known from nuclear physics. 
The Hartree-Fock ground state (vacuum) contains the quark condensate. It is determined through 
a nontrivial solution of the gap equation. Mesons, being quark-antiquark excitations of the vacuum, 
behave very much like particle-hole states in a strongly interacting many-body system. 

In the first Section of this Chapter we present details of the calculation of the zeroth-order 
correlation functions at finite T and /x. The second Section contains further discussion concerning 
the regularization of the model. In the third Section we define the full correlation function in the 
RPA approximation and discuss the difference between the dynamic and screening masses. Finally, 
we show that at sufficiently large T or /r, the dynamic masses of chiral partners become equal. Such 
degeneracy is a characteristic feature of the chiral symmetry restoration phase transition. 


6.1 Zeroth-Order Correlation Functions 


Our investigation of the meson properties in the NJL model starts with the calculation of the 
meson correlation functions. At first, using Eqs. ( |4.25 ) and ( 4.261) , we construct the zeroth-order 
correlation functions. They correspond to the following quark loop diagram Q 


=*tr 


d'^p 


[Tat^G{p + Q)Tbt^G{p)]. 


( 6 . 1 ) 


Here = (a;,q) is the external (meson) momentum. The indices A and B take on the values P 
or S and Fp, F^ are the Dirac matrices in the spinor space: Fp = iys and Fg = 1. We note that 
the indices a and b in Eq. (4.26) correspond to two pairs of indices in Eq. (6.1), namely a = {A,i) 
and b = see also Fig. |5]. 

^Conventions introduced in this Chapter are suitable for the two-flavour symmetric model defined in Sections 3.2 
and 5.2. Nevertheless, as we have seen in Chapter 5, changing Nf and/or Gs one can easily obtain expressions whose 
form is valid for other versions of the model. 
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Figure 6.1; Quark loop diagram corresponding to the zeroth-order meson correlation function. 

fO) ij 

In the isospin symmetric case the matrix Xab diagonal in flavour indices and its all diagonal 
elements are the same. Consequently, we can write x^ab^ ~ ^^^Xab^ where 


X%{Q) = 2iN,Sp 


d^p 


[TAGip + Q)TBGip)] . 


( 6 . 2 ) 


Calculation of the trace over the Dirac indices, denoted here by Sp, shows that the nondiagonal 
terms {A ^ B) vanish. Therefore, we are left with two non-vanishing functions 


x'pUQ) = SiN, J ^ 


p2 _ ]\^2 _j_ p . Q 


(27r)^ [{p + QY — + ie ] [p2 — -|- ie ] 


and 


(°Vni = «-/V f _ p'^ + M‘^+p-Q _ 

XssW) j + Q^2 _ + ie]\p^ - + ie]' 


(6.3) 

(6.4) 


6.1.1 Vacuum Parts 


Using the imaginary-time formalism, discussed in detail in Chapter 4, we generalize expressions 
( |6.3D and (6.4) to the case of finite temperature (density). To do it, we replace the integration 
over energy in Eqs. (|6.3|) and (^) by the Matsubara sums (|4.18|), and convert these sums to 
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the contour integrals using expressions ( 4.23 ) and ( 4.24 ). In this case, the vacuum parts of the 
correlation functions can be represented in a compact form 


,(0) 


and 


.(0) 


Q^l2,vac(-^^) Q‘^) 

(6.5) 

-4M2)/2, vac (M2,Q2)_ 

(6.6) 


The function appearing in Eqs. (^) and (^), has been already defined by formula 

(5.7). The new function l 2 ,vac{M‘^, Q'^) is defined by the integral 


/2,vac(M2,-4) = -4iV, J 


d^PE 


{2eY [{pE + qE/2Y + M^][{pE - qE/2Y + M^] ’ 


(6.7) 


where p 4 = —ip^, q^ = —iuj and = dp^d^p. We note that frequency is purely imaginary 
now [Q^ = to ^ iqi), thus the “Euclidean” form of expression (§3) is natural here. 


6.1.2 Temperature Parts 

The form of the temperature part of x^ppiQ) follows from Eqs. (|6.3|) and (4.23) 


zoo+e 


/ ir 


dp^ f d^p 


—2cx)+e 


+ 


271 J {27rY [{p + QY — {p — QY — \ + 1 


zco+e 


—8iNr 


dp^ f d^p 


—zoo+e 


p-Q 


P-Q 


2nJ {27rY [[{p+QY — M"^] [{p—QY—M'^][p‘^ — M‘^]\ + 1 


( 6 . 8 ) 


Using again the fact that is the bosonic Matsubara frequency (Q^ = to ^ iujm = 2mEiT), we 
make the shifts of the integration variable p^, and simplify expression (|6.^) to the form|^ 


zoo+e 

X?i,tem(^) = J 

—zoo+e 


dp^ 


d?p 1 1 

(27r)^ p^ — eP °+ 1 


zoo+e 

+8iiVeQ2 J 

—zoo+e 


dp^ 


d^p 1 1 

(2+3 [(p + Q/2)2 - M2 ][(p - Q/2Y - M2] eiP°+i^rnl2)IT + i' 

(6.9) 


^At this place we write uj = iuim rather than uj = iq^ in order to point out the discrete character of the complex 
variable ui. This fact allows us to make the shifts of the argument of the exponent. 
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The last result indicates that the analogous decomposition to that find in vacuum, see Eq. (|6.5D 
holds also at finite temperature, namely 


Xpptem(^) ~ -^i,tem(-^ ) Q -^ 2 ,tern(-^ 5 Q)i 


( 6 . 10 ) 


where .fi^tem(-^"^) defined by Eq. ( 5.13 ) and the function ^Q) form 


ioo+e p 

4,tem(^^Q) = -8^iVc / 

— 200+6 


d^p 1 1 

(27r)3 [(p+Q/2)2-M2][(p-Q/ 2)2-M2] e(p°+i‘^-/2)/T + i' 

( 6 . 11 ) 


Since -^2,tem(-^^) Q) is a function of o;^ and we shall also use the notation -^ 2 ,tern 
In the special cases (i.e., either for g = 0 or a; = 0) one finds 


oe 

N f 

-^ 2 ,tern (-^ 5^ ) 0) = J 


and 


+tem(^^0,(?2) = _^ j 


Ne 

f dpp^ 1 

1 

(6.12) 

^2 

J Ep e^vIT + 1 El 

0 P 

1 , ,2 

- 4W 


00 



g7r2 

f dpp 1 

J Ep e^vIT + 1 

2p-q 

2p + q ' 

(6.13) 


We note that Eq. (3.12) has been obtained for a purely imaginary frequency u> = iuJm- Thus, the 
expression on the RHS of (|6.12| ) is a well defined quantity. Eurthermore, according to the rules 
discussed by us in Chapter 4, the result ( 6.12| ) can be analytically continued to real frequencies. 
This can be achieved by making the substitution icom —> w ± ie, where cu is real. 

In the similar way one does the calculations in the scalar channel finding the following decom¬ 
position 


^5S,tem(^) ~'^i.tem(-^^) iQ‘^ 4-^^)-^2,tem(-^^+)) 
which is the analog of Eq. (|6.6[). 


(6.14) 


6.1.3 Matter Parts 

In order to calculate the matter part of the correlation function Xpp(Q)) we use Eqs. (O) and 


(4.24), which yield 

^pp,mat(^) ^ J 


f dp^ 

r d^p 

J 27r J 

(2vr)3 


+ 


p-Q 


{p + QY-M"^ [(p+Q)2-M2][p2_M2] 


(6.15) 
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Making the shifts of the integration variable one can check that the decomposition of the form 
(|6.5D or (6.1C|) is valid at finite density. Indeed, after a few manipulations we obtain 


4p,mat(^) = A.mat(^') " g'/2,mat(^^ Q), (6-16) 

where Ii niat(-^^) defined by Eq. (5.14) and the function I 2 mat(-^^’ Q) given by expression 


^{M\Q)=AiN, I ^ 


' 2 ,ma 


d?p 1 

(27r)3 [{p + Q/2)2 - M2] [{p - Q/2)2 - M ^]' 


(6.17) 


The calculation of ^2 o;^, q^) proceeds in the same way as the calculation of the function 

/imat(-^^)- Now, as remarked above, the integral over p^ is evaluated for imaginary io and 
subsequently analytically continued to real frequencies. The final results (for a; > 0) are 




N /■ 

-^2,mat(-^^) 0) = J 


dpp^ 




27r2 


log ie) 


(6.18) 


and 


PF 

-^2,mat(^^^)0>9^) = “ 2 ^/’ 


Nc r 

2g7r2 J Ep 

0 

N, 


In 


2p — q 


2p + q 


2qTT- 


:G{M\q^), 


(6.19) 


where Ep = \Jpp + 


±ie) = - 


V4M2 — a;2 arctan ( - , 

KEfV^M^ 

Epy/uj"^ — 4M2 + ppoj 
Epy/uj"^ — 4M2 — ppijj 


- \lup' — 4M2 In 


±i—Q{uj- 2M) e {2Ep - Lo) \JuP - 4M2 
2a; 


(a; < 2M) 


(a; > 2M) 


( 6 . 20 ) 


and 
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G{M\q^) = 


' + 4M2 


— Ep In 


2pF + q 


yjq^ + 4M2 
H-In 


2 pi;’ - q 

2M^ + gpi? + Ep^Jq'^ + AM'^ 

2M2 - qpp + + 4^ 


— g'ln 


Pi? + Ep 


M 


( 6 . 21 ) 


It becomes clear now, that for the scalar channel we hnd the analog of Eq. (6.14), namely 


xSmat(^) = A,mat(^') " iQ" " 4M^)4.mat(^^ Q)- 


( 6 . 22 ) 


6.2 Regularization of the Zeroth-Order Correlation Functions 


In the last Section, we have reduced the problem of the evaluation of zeroth order correlation 
functions to the calculation of the functions Ii = /i,vac + med ^*4 h = l 2 ,Ya.c + I 2 med' 
theless, one can easily notice that the function l 2 ,Yac{M‘^, —qp) is defined by the divergent integral 
Moreover, although the functions /2 j^g(^(M^, g^) are finite, they should be regularized 




in the same way as the vacuum functions. In this way, see discussion in Section 5.1 following 
Eq. ( 5.15|) , we ensure the correct high temperature (density) behaviour of the complete expression 
h = l 2 ,Yac +12 med' in the next three Subsections, we shall present three different ways of 

the regularization of the function l 2 {M‘^, 00 “^,q^). We shall use the methods introduced in Chapter 
5. 


6.2.1 3-Dimensional Regularization 

At first, let us discuss the simplest method, which is based on the introduction of a 3-dimensional 
cutoff. Similarly as it was done in Section 5.1, we perform integration over energy in Eq. O' 
Setting q=0 and collecting the contributions from the poles one finds 


h,Yac{M‘^, —ql) 


Nc f dpp^ 1 

Ep Ej + lqX 


(6.23) 


Doing analytic continuation of expression ( 6.23 ), q^ —> —i{uj ± ie) 
p-integration we find {oj > 0) 


and introducing the cutoff for 


jR,3D 

-' 2 ,vac 




Nc f dpp^ 1 

2^1 E, E}-Hu,±kr 


(6.24) 


This integral can be calculated in the same way as expression ( 6.1 j ) — there is only an overall 
difference in sign between these two integrals and pp in (3.18) corresponds to A in (6.24). 
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The final formula for the regularized function l 2 {M‘^,u;‘^ ± 0) 


IS 


Nc f dpp 


h [M ,u; ±^e,0) = J 

for finite temperature (/r = 0 ), and 


.2 r 


En 


1 - 


gBp/T + iJ E^-l{uj±ie)^ 


(6.25) 


2 , • n^ f dpp‘^ 1 

I, (M.u, ±«,0) = -^y— 


(6.26) 


PF 


for finite density (T = 0). We note that integral (6.2(:) also can be done analytically, similarly as 
Eqs. (| 6 .lj ) or (H). 


Doing calculations for = 0 and q 7 ^ 0 one finds the analogs of (6.25) and (3.26), namely 


and 


jR,3D 

-'2 


(M2,0,(?2) 


A 

Nc f dpp ' 2 

27r2g J Ep + 1 
0 


- g 

2p + g 


jR,3D 

-'2 




Nc r d^ 2p - g 
27 r 2 g J Ep 2p + q 

Pf 


(6.27) 


(6.28) 


Although the 3-dimensional cutoff regularization is frequently used, for many purposes it is not 
satisfactory since it explicitly breaks the Lorentz invariance. We have collected the results based 
on this regularization in order to facilitate the comparison of different formulations of the model. 
Nevertheless, the major part of our further calculations will employ the covariant Pauli-Villars 
method. 


6.2.2 Pauli-Villars Subtraction Scheme 


In this Subsection we shall calculate the function I 2 {M‘^,Q) = / 2 ,vac(Af^, Q^) + I 2 j^g(^(M^, Q) 
using the Pauli-Villars regularization scheme. This calculation proceeds in the similar way as the 
evaluation of the function [M'^) in Section 5.2. The regularized form of (3.7) is obtained by 

making the substitution 


1\ 

h,vaciM‘^: —Qe) ~Qe) = ^ih,yacjEj, —q^)^ 


(6.29) 


i =0 


where the coefficients A* are defined by formula ( 5.25 ). To calculate 7^(M^, — g|.) one has to 
use first the identity 


1 

AB 




(6.30) 
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which allows us to simplify the form of the denominator in the integral on the RHS of Eq. B- 
Afterwards exchanging the order of integration over pE and u and introducing the new variable 
p'e ~ Pe + '^Qe/‘^ obtains the expressions which can be easily integrated. The final result is 


~^e) = 

i=0 




'1 + 


lE 

2Ai 




(6.31) 


The function , uJ^) can be obtained by performing the analytic continuation of the 


ex¬ 


pression given on the RHS of Eq. (3.31). The substitution q —> —i(a; ± ie) leads to the following 
formula (a; > 0) 


Nr 


N 


4 Vac ± ^ E 


i=0 


2 Aj I f io 


1- 


u 


\2Ai 


arcsm 


iO 


-|- InAj 


+ 0(w —2Aj) 




(6.32) 


The functions /2 med(-^^’should be calculated in the same way as the function nied^-^^)’ 
namely we make the replacement 


4,mediate) ^ = E V4,med(A?,l3), 

’ i=0 

where the function med(-^i ’ defined by Eqs. (5.11) or (|6.17|). 


(6.33) 


6.2.3 Schwinger Proper-Time Regularization 

Using the Schwinger proper-time method of regularization, we rewrite Eq. 
form 


(6.7) in the following 


tR,SPT 

-' 2 ,vac 


{M\ 


- 4 ) 


-4A^c 

_El 


d^PE 

(2vr) 


CXJ i 

J ds s J dttexp s[p|;-|--|-u(l — u)g|;]| 


du El 


A-2 0 

-j2+u{l-u)j^ 


(6.34) 


Replacing the integral over energy, p^, by the sum over Matsubara frequencies we can generalize 
Eq. (|6.34|) to the case of finite temperatures 
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d^p 


„ (27r)3 7 a-2 

0 


ds s du X 


exp s 
NrT 


M" + u{l - u){{E>^^Y + q") + [p - q(l - u)Y + - S^(l - u) 


27r2 A 


E duEi 


J 0 


p'f _ T?i> I 
?f Tpb I 


l2 


I n-K{i-u)Y 


where = 27rnT and E^- = (2j + l)7rT. 


(6.35) 


6.3 Meson Masses 


6.3.1 Definitions 


In the random phase approximation (RPA) the full correlation function has the form 


Xaa{Q) 


xaVQ) 

1 - Gsx^AiQ) 


(6.36) 


This is a nonperturbative expression, which follows from the infinite summation of a certain class 
of diagrams, see Fig. [6.2]. 



Figure 6.2: The class of diagrams used to evaluate the full correlation function in the RPA approximation. 
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The correlation function XAAiQ) depends on Q through the variables and q^. Therefore, 
in the following we shall use the notation xaa{^‘^,q‘^) rather than xaa{Q)- The dynamic mass is 
defined by the position of the nearest to zero pole of the function xaa((^^, 0)- In the NJL model it 
can be easily found by solving the equation 


1 - GsXAAi^Ln’O) = 0 


(6.37) 


This comes directly from Eq. ( 3.36 ). On the other hand, the screening mass is defined by the 
asymptotic behaviour of the static correlation function in space, namely 


m.o 


= — lim 


din XAAjr) 
dr 


(6.38) 


where 


XAA{r) = 


d^q 

(^ 


H-cxd 


XAA{0,q^) 6 *^ 1 - = 


dvr^ir 


dqq XAA{0,q‘^) e 


iqr 


(6.39) 


To get the last term in Eq. ( |6.39 ) we first performed the integration over the angles, and then 
extended the region of integration over q = | q | to the interval from minus to plus infinity. 

Let us now consider an arbitrary correlation function xaa{'^‘^i not necessarily obtained from 
the NJL model. One can ask the question when the dynamic mass is equal to the screening mass. 
At T = 0 the system is explicitly Lorentz invariant and XAAi<^‘^,0) = Xaa(0, —w^). Therefore, if 
the function XAAi^"^-,^) has a pole for u) = mdym the function xaa(0, q^) has a pole for q = im^yn- 
The latter gives the contribution to the integral (6.3£), which has the form ~ exp^—rndyni")- For 
very large r this contribution is the dominant one, since im^yn is the nearest pole. Consequently, 
we find that rndyn = mscr- Of course one has to be careful because the analytic structure of the 
correlation function can be complicated and there might be some other contributions to the integral 
( 5.3S ). Nevertheless at T = 0 we expect that mdyn = fnscr- A different situation takes place when 
T > 0. In this case it is difficult to find a simple relation between the functions xaa{^‘^,0) and 
Xaa{0^ The Lorentz invariance is implicitly broken by the existence of the preferable reference 
frame connected with the heat bath. The fact that XAAi<^'^,0) has a pole for ta = mdyn does not 
imply that X/ia( 0, g^) has a pole for q = imdyn- Eurthermore, at T > 0 the contribution to Eq. 
( |6.39| ) from the cuts of the correlation function may be important and it can substantially change 
the space asymptotics of the function XAAir). In consequence, at T > 0, because of these at least 
two reasons, the dynamic mass and the screening one are different. Of course, the same arguments 
hold for the case of finite density. 
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6.3.2 In-Medium Meson Dynamic Masses 


The dynamic masses of the pion and the sigma are obtained from Eq. (6.37), which has the form 


I-G5 




= ^ + Gs 0 ) = 0 


for the pseudoscalar (pion) channel, and 


(6.40) 


1 - Gs Ii(m2) - 


m 


= + {^dyn,a “ 4M ) Gs h{M , 0) = 0 


(6.41) 


for the scalar (sigma) channel. In Eqs. ( 6.4C| ) and ( 6.41j ), we have assumed that the parameters 
are such that the gap equation has a non-trivial solution, i.e., M = m + GsM/i(M^). 

Formula ( \ 6 . 4 Oi) is the cornerstone of the NJL model. It explicitly shows, that the spontaneous 
symmetry breaking leads to the appearance of the Goldstone bosons. If the current quark mass is 
zero, m = 0, the pion dynamic mass turns out to be zero. Moreover, at the Hartree level the mass 
of the sigma is simply 2M. We also note that the structure of Eqs. ( |6.40 ) and (6.41) is independent 
of the regularization scheme (as long as our regularization procedure modifies the functions Ii and 
h only.) 

We note that Eqs. ( 6.4C| ) and ( |6.41 ) are correct not only in vacuum but also at finite temperature 
or density. As follows from Eqs. ( |6.32 ), ( |6.12 ) and (|6.1^ ), the correlation functions have cuts for 
arguments larger than 2M. At high temperatures (densities) the meson poles merge with the 
quark-antiquark cuts, which means that in this model the mesons can decay into qq pairs. Thus, 
at these temperatures (densities) there are no isolated poles which can be identified with the mass 
of the pion or the sigma. We circumvent this difficulty by defining the mass as the zero of the real 


part of Eq. (|6.37|) , i.e.. 


l-GsRexS("^iyn,0) =0. (6.42) 

We note that by using Eq. ( |6.42|) we implicitly neglect the qq widths of the mesons |^. This seems 
reasonable, since these widths are non-zero only because the NJL model lacks confinement. 

Our result concerning the temperature dependence of the dynamic masses of mesons is shown 
in Eig. [6.3]. The calculations are based on the two-flavour symmetric model. We use now the 
same computational scheme as that used to calculate the temperature dependence of M in Section 
5.2, i.e., we make three subtractions with the regulating masses: Ai = 680 MeV, A 2 = 2.IAi and 

®In a general situation, the mass and the width of a meson are determined by the shape of its spectral function. 
In our case the latter has the form (2 Imx*'°^)/(1 — GsX^°'’) = (2 Gs Imx^°^)/( [1 — GsRex^°^]^ + [Gs Irax*-*^^]^). If 
the imaginary part of x^°^ (describing t he wi dth of a meson) is small, the spectral function is strongly peaked for the 
four-momenta satisfying the condition (|6.42|). 
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Figure 6.3: Temperature dependence of the dynamic masses of pion (solid line) and sigma (dashed line). 

A 3 = 2 .IA 2 . The coupling constant Gs = 0.75fm^ and the current quark mass m = 8.56 MeV. 
Using these values of the parameters we find that at zero temperature the constituent quark mass 
M = Mq = 376 MeV, the pion mass mdyn,-K = 138 MeV and the sigma mass mdyn,o- = 760 MeV. 
The pion and sigma masses are connected by relation 

mlyn,a ~ (6.43) 

In Fig. [6.3] one can see that the pion mass remains constant as long as the temperature is smaller 
than 250 MeV. Afterwards it increases suddenly. The temperature dependence of the sigma mass 
is rather complicated. In the interval 0 < T < 250 MeV it behaves in the similar way as the 
constituent quark mass, see Fig. [5.6]. Moreover, in this region relation ( |6.43| ) is still fulfilled. 
When the temperature reaches 300 MeV, the sigma mass stops decreasing, remains for a while 
constant and later increases. At very high temperature the sigma mass and the pion mass are with 
a good approximation the same. This fact signifies the chiral restoration. 

In Fig. [6.4] we show the density dependence of the dynamic masses of mesons. We use here 
the same values of the parameters as in the finite temperature calculations discussed above. Thus, 
our results for g. = 0 and T = 0 coincide with each other. One can see that with the increasing 
value of the chemical potential the dynamic masses of pion and sigma become degenerate. This 
fact indicates restoration of the chiral symmetry at high density. The overall //-dependence shown 
in Fig. [6.4] is similar to the T-dependence shown in Fig. [6.3]. 

The first calculation of the in-medium dynamic masses of the pion and the sigma in the NJL 
model was done by Hatsuda and Kunihiro [^ . Their calculation employed the 3-dimensional regu- 
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larization scheme. The results presented in Figs. [6.3] and [6.4] follow from the calculation based 
on the Panli-Villars method [^]. They qualitatively agree with those obtained in Q. The ad¬ 
vantage of using the Panli-Villars method in our case is that it will allow us to compare directly 
the temperature dependence of the dynamic and screening masses. As follows from the discussion 
in the next Chapters, the 3-dimensional regnlarization scheme introduces unphysical contributions 
to the Fourier transform defining the static correlation function. Thus, this way of regularization 
spoils the expected relations between the screening and dynamic masses. 



Fignre 6.4: Pion (solid line) and sigma (dashed line) dynamic masses plotted as functions of the quark 
chemical potential /i. 
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Chapter 7 


Static Meson Fields at Finite 
Temperature 


In this Chapter we study the damping of static meson correlation functions at finite temperature. 
Our interest in these investigations has been triggered by the results of the lattice simulations 
of QCD. They indicate that high-temperature hadronic correlation functions describe a system of 
(almost) non-interacting quarks. In the first Section of this Chapter we do the analytic calculations 
based on the perturbative QCD. This allows us to understand the high temperature behaviour of the 
correlation functions. In the second Section we come back to the study of the NJL model. We find 
that the screening masses obtained within the model have the desired temperature dependence. In 
particular, for stable mesons the screening masses at T = 0 are equal to the dynamic ones, whereas 
for extremely high temperatures they approach the values expected from the perturbative QCD 
calculation. 


7.1 Screening of Meson Fields in Hot QCD 


Lattice “measurements” of hadronic correlation functions in QCD show exponential behaviour at 


At high temperatures, T > Tc, the hadronic screening masses 


large distances 11, 1, 

of chiral partners are equal, but non-vanishing. This is consistent with the restoration of chiral 
symmetry. Eletskii and Ioffe [^] pointed out that the meson screening masses approached the 
value 27rT, whereas those of baryons were close to 3ttT and argued that such a behaviour is typical 
for a gas of non-interacting massless quarks. Their reasoning was the following: In the considered 
case the space-time correlation function is the product of two (mesons) or three (baryons) quark 
propagators. The space asymptotics of each propagator is determined by the lowest lying pole 
in the complex q = |q| plane. In the imaginary-time formalism this pole corresponds to the 
lowest fermionic Matsubara frequency vrT and, therefore, the screening masses are multiples of vrT. 
Shuryak has presented alternative arguments, based on the large-distance behaviour of the free 
quark propagator at finite temperatures, which lead to the same result |^]. 

We should emphasize that the screening masses obtained from the lattice calculations are not 
exactly the multiples of ttT. For baryons and (axial)vector mesons the relation of Eletskii and Ioffe 
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is indeed very well fulfilled. However, for (pseudo)scalar mesons one can see a deviation from 2'kT. 
This discrepancy indicates that there is a non-negligible interaction in this channel, present even 
at very high temperature. Nevertheless, in the limit T —> oo we expect that in all meson channels 
the screening masses would go to 2t:T. 

The arguments of Eletskii and Ioffe allow us to understand intuitively how the screening masses 
turn out to be 27rT and SvrT. We are of the opinion, however, that it is instructive to obtain 
this result from a more rigorous calculation. Actually, the aim of this Section is to show that 
starting from the regularized momentum representation of the static meson correlation function we 


can analytically derive its full spatial dependence [52|. Knowing the asymptotic behaviour of the 
correlation function, we can find the values of the screening masses which are in agreement with 
the result of Eletskii and Ioffe. 

The spatial dependence of the correlation function will be calculated by us in two ways, namely 




2 ^^ ^iqzZ 


or 






An'^ir 




(7.1) 

(7.2) 


The quantity x(0, g^) denotes the static meson correlation function obtained in the framework of 
the perturbative QCD. We shall restrict our considerations to the special case of the pseudoscalar 
channel. Eq. (^) describes the situation in which the correlation function is used to find the 


response of the system to the static plane-like perturbation. Here z is the distance from the plane 
in the transverse direction and Qz is identified with the momentum in this direction. Such a case 
is considered in the lattice simulations. Eq. Q, equivalent to the definition (|6.39|) , describes 
the situation when the perturbation is the static point. We note that the functions x{z) and x(r) 
depend differently on their arguments, however, for simplicity of notation we use for them the same 
symbol. 

In the perturbative QCD, the leading contribution to the meson correlation function is given by 
the quark loop diagram, see Eig. [6.1] (the diagrams with gluon lines contribute only in higher order 
terms). Thus, we can use our results from the previous Chapters. In particular, if we employ the 
Pauli-Villars regularization procedure, we can use expressions (|6.5| ) and (|6.10| ) with the functions 
A,vac(Ar2), tem(^^)> -q^) and 72,tem(^^> given by expressions (|5.24D , (^.13|) , 

( |6.31| ) and ( [6.13 ), respectively. 

In practice, we do not use Eqs. ( [7.1| ) and {7_^) directly in the form as they have been written 
down. To evaluate the Fourier transform, in both cases, we analytically continue 1^6 

whole complex q plane. Later we deform the contour of the integration, and represent the result as 


the integral around the singularities appearing in the upper half-plane. The Fourier integral (7.1) 
or 0 becomes in this case the sum of two terms (the first one represents the integral over the 
vacuum part and the second one is the integral over the temperature part). 
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7.1.1 Vacuum Part 


The vacuum part of the correlation functions Xvac(0, q^) is given by expression 


N ^ 

Xvac (0)9^) = X/ ^ 

2 = 0 



+ln Aj 


(7.3) 


which follows from Eqs. (6.5), (5.24) and (6.31). The analytic continuation of (7.3) has cuts for 
purely imaginary q = ik, starting at k = ±2M, see Fig. [7.1]. Consequently, the Fourier integral 
can be replaced by the integral around the cut stretching along the imaginary axis. In this way we 
obtain 


Xvac(2:) 


OO 



2M 


Xvac(O) 


—k'^ + ie) — Xvac(0, —k^ 



(7.4) 



Figure 7.1: The cuts of the vacuum part of the correlation function. They start at Imq = ±2M and go 
along the imaginary axis to ±oo. The integrations along the curves Ci and C 2 are equivalent. 


The real parts of Xvac(O) i ie) are the same on both sides of the cut, and the imaginary 
parts have opposite signs. Thus, the expression in the square brackets on the RHS of (^) equals 
2ilmxvac(0, —k‘^ + ie). As the final result one gets 
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(7.5) 


N ^ r /- 

Xvaciz) = J dkk^k^ - 4A| 

*=° 2Ai 


-kz 


N 


e = 


TT^Z 


J2AiA^K2{2Ai 

i=0 


where K 2 is the modified Bessel function [53|. The vacuum part of the correlation function 
Xvac(0,(?^) is divergent and must be regularized. Nevertheless, the Fourier transform remains 
finite for z > 0 also when one sends the cntoff masses to infinity. For Aj —> 00 (i > 0), only the 
first term in the above series is non-negligible and we find 


N, 


Xvac(^) = ^M^K 2 { 2 Mz). 

■K^Z 


Similar calculations starting from Eq. Q lead to the result 


XvacC?’) = 


Nr 


2TT^r^ 


[ 37 X 2 (2Mr) + 2MrKi{2Mr)] 


(7.6) 


(7.7) 


7.1.2 Temperature Part 

The temperature part of the correlation function is given by expression 

00 „ 

Nc f dpp"^ 1 


Xtem(6)? ) — 7]-^ y 


R eEp/T^l 


4 + - In 

2p — q ' 

L P 

2p + q \ 


(7.8) 


which follows from Eqs. dH), ( ^.131) and (|6.13D . In the following we shall make the snbstitution 

2p-q 


In 


2p + q 


_ 1 {2p - qf + 

2 ^ {2p + qY + ’ 


(7.9) 


In this way Xtem(9)9^) becomes an analytic function of q. Of course, at the final stage of the 
calcnlations we have to consider the limiting case e —> 0. 

The spatial dependence of the temperature part of the correlation function is obtained from 
Eqs. ( [7T|) and (^). Exchanging the order of the integration over q and p we get 


00 r, 

_ f dpp^ F{p,z) 

Xteml^J - 


27r^ J Ep eEvN -|- \ ’ 


(7.10) 


where the function F(p, z) is defined through the integral 


F{p,z) = 


dq 


pz 


— 00 
TT 


4++ 


^iqz 


2p {2p + q)'^ + 

(i - 2pz) - {i + 2pz) 


(7.11) 
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The calculation of F{p,z) proceeds as follows. First we analytically continue the integrand on 
the upper half-plane of complex q. Such a function has two cuts starting at g = ±2p -|- ie and 
going upward, parallel the imaginary axis (see Fig. [7.2]). Consequently, the initial integral can 
be replaced by the two integrals around the cuts. The values of the integrated function on both 
sides of the cuts are different only because the logarithm has a different phase there. However the 
difference of the phase, the quantity which enters our integrals, is constant. This simplifies our 
expressions and leads, in the limit e —> 0, to the result shown in the second line of Eq. (7.11). 



Figure 7.2: The analytic structure of the function being the integrand in Eq. (7.11). The integration along 
the real axis (contour C) can be replaced by the two integrals around the cuts (contours Ci and C 2 ). 


Substitution of Eq. ([7.11 ) into Eq. ( 7.10 ) leads to the following result 


N, 


Xtem(^) — o o ( „2 ^ 


1 d \ d f dp 


27r2 zdzj dzJ EpC^p/^ + l 


cos{2pz) 


(7.12) 


Now, using the series representation of the Fermi-Dirac distribution 


+ 00 


1 1 > O/ 

(21 + l)27r2 

/=—00 ^ ' 


e* -b 


(7.13) 


we can write the temperature part in the form 
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Nr f I 19 




=(tem/ 


(7.14) 


where 


=(tem/ \ _ ^ 




dz Jo 2E, 


dz) 1 3 

cos{2pz) = -—Ko{2Mz) = —MKi{2Mz), 


2dz 


(7.15) 


and 


= -TiZ Y: 


d 


H-C50 


dp cos{2pz) 


+ 00 


dz 


l = — OQ 
+ O0 


Jo {21 + + El 




dppe 


2ipz 


l=—oo 


/-oo {21 + I)27r2r2 + El 


= ttT Y. exp(-2zY/(2Z + l)27r2T2 + M2). 

/ = —OO 

In the 3-dimensional case one can derive an analogous expression 


(7.16) 


Xtem(^) ~ 


Nc I 2 


2d 1 92 


Q O r\ “I” 


I /^3 r^2 p 3z‘^ 


(j,tem(^) ^ (^iem(^) 


tem,- 


(7.17) 


7.1.3 Asymptotics 


Using Eqs. (7.6) and (7.14), and also employing the recursion relations for the Bessel functions 
[ p4[ we can check that in the full correlation function the vacuum part and the first piece of the 
temperature part (coming from the differentiation of the function G^®™) cancel each other. 0 We 
are left only with the term coming from the differentiation of the function Moreover, for 

large distances, it is enough to consider only the leading terms in the series which defines G^^™ 
{1 = 0 and I = —1). Consequently, we find for large z and finite M that 


X{z) 


2NrT 


'KZ 


\/7r2r2 M2 exp(-2zv'7r2r2 -hM2). 


(7.18) 


From Eq. ( [7.18 ) we read off the screening mass 
is obtained if we use Eqs. (TJT) and (7.17). 


= 2Vvr2r2 -kM 2. The same screening mass 


Let us also discuss in more detail the case of massless quarks. The expression for Xtem(^) 
Xtem(^) then given in a compact form because the series appearing in the definition of 

(^tem t)ecomes a geometric one and can be easily summed up. Using Eq. (7.14) we obtain 


_+ ^"r.coth(2^r.)]. (7.19) 

^Such cancellation takes place only because the vacuum part was calculated in the limit Ai ^ 00 (i > 0). If we 
kept Ai’s finite, we would have to regularize the temperature part as well. Otherwise, we could not get the correct 
result for the screening masses. 
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and similarly from Eq. ( 7.17 ) we get 


Xtem(^) 


M^o 


3iVc N.T 

27r^r^sinh(27rTr) 


1 + 2'nTr coth(27rTr) 


+27r2rV2 


sinm(27rrr) 



(7.20) 


Eqs. (7.6) and ([7.7|) in the limiting case M 


0 have the form 


and 


Xvac(2) 


M^o 


Nc 

2'ii'^z^ 


Xvcic{r) 


M^O 


3N, 

dvr^r® 


(7.21) 

(7.22) 


In this case we can easily see that the temperature parts of the correlation functions vanish if T = 0. 
On the other hand, in the expression for the complete correlation function, for each finite value of 
T, the temperature independent piece of the temperature part cancels exactly the vacuum part. 


7.1.4 Remarks on the Temporal Function 

Einally, let us briefly discuss the screening in the temporal direction. In this case we need to 
evaluate the correlation function in the imaginary time direction 


where 


x(")(T)=r^e—'*"x^(za;„,0). 


^R(a;,0) = ^^ dpp^ tanh + (reg.) 


is the regularized correlation function (“reg.” denotes the regularization terms). 
Using standard techniques we find 


(7.23) 

(7.2d) 






TT^ 


dpp' 


2 cosh [2£' p(t — (3/2)] 
cosh‘s {Ep(3/2) 


+ (reg.). 


(7.25) 


where (3 = 1/T. Since the integral converges for 0 < r < /3 we can again drop the regulating terms 
by letting Aj —> oo. For massless quarks, the integral can be done analytically 




NcT^ 52 A 2 t\ tt 

27r2 ci(r//3)^ V (3 ) sin(27rr//3) 


(7.26) 


The real time correlation function is obtained by analytically continuing to real times; for large t 
we find x^^^(t) ~ exp(—27rTt). Thus, for massless non-interacting quarks, the mesonic screening 
mass in the time direction is the same as in the spatial direction, ruscr = 2ttT. Consequently, 
we expect the temporal correlation function in mesonic channels to behave qualitatively in the 
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following manner; For a uniform source at zero temperature, where mesonic modes exist, 
oscillates in time with frequency m, the meson mass. On the other hand, at high temperatures, 
where the system consists of weakly interacting almost massless quarks, the correlation function is 
screened, with the screening mass 27rT. 


7.2 Screening in the NJL Model 

In this Section we analyze the spatial dependence of static meson correlation functions calculated 
within the NJL model. We use again the two-flavour symmetric model defined in Sections 3.2 
and 5.2. Thus, our considerations are restricted to the scalar and pseudoscalar channels. We find 
that the screening masses (obtained from the asymptotic behaviour of the correlation functions) 
differ from the dynamic masses (defined by a pole of the meson propagator in the energy plane, 
see discussions in Section 6.3). In the high temperature limit, the screening masses approach 27rT, 
which corresponds to a gas of non-interacting quarks. Nevertheless, the interaction effects remain 
well beyond the chiral transition temperature. The overall temperature dependence of the screening 
masses is in the agreement with lattice results. 


7.2.1 Analytic Structure in Complex Momentum Space 


Similarly to the QCD case considered in the previous Section, the Fourier transform ( 6.39| ) can be 
rewritten as a sum of a few contributions, which appear due to the singularities of the correlation 
function in the complex g-plane. Thus, at first we are going to study the analytic structure of the 
correlation functions. 


The results of Section 6.1 indicate that Xaa vac(0’ 9^) following we refer to these 

as the vacuum cuts) along the imaginary axis starting at q = ±2iM and going to ±ioo. In order to 
find the analytic structure of the function tejji(O) 9^) we study the function l 2 ,tem(-^^) Oj 9 ^) 
in the whole complex q = |q| plane. To this end we represent /2 tem (-^^5 *3) 9^) ^ ^wo 

functions, namely 


where 

OO 

(±) Nc r dpp 1 2p- z±ie 

2 ,tern 2zTr‘^ J Ep + 1 —2p — z Eie 

0 


(7.27) 

(7.28) 


At the end of the calculations we let the infinitesimal e go to zero. The functions have 

logarithmic cuts parallel to the real axis and stretching from minus to plus infinity. At the respective 
cut the imaginary part of the function 72^^gj^(z) is discontinuous. On the physical Riemann sheet, 
this amounts to a change in sign: 




NcT 

2qRTT 


In 


1 J- exp - 


\/m + gl/4' 


(7.29) 
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Here qR is real, 5+ = 2e and <5_ = 0. Thus, for the upper sign one is above the cut of and 

for the lower sign below. Obviously the imaginary part of ^^ove and below its cut is equal 

to that of .^ 2 ^^em below its cut. This can be summarized by the following equations 


Im/, 


(-) 

2 ,tern 


{z = qR-i5^) = Im (^ = lR + ^^±) 




(7.30) 


The cuts are arranged in such a way that for z on the real axis, the imaginary part of -^ 2 ^^em(^) 
cancels that of -^ 2 ^^em^^^ -^ 2 ,tem (-^^5 9 should be. Concluding, we 

state that the function x'^^\ tem^^’ following we refer to these as the temperature 

cuts) which are parallel to the real axis ^ Thus, the full correlation function XAyl(0,g^) defined 
by Eq. (6 .St ) has at least these two cuts. Moreover, it can have poles for imaginary arguments in 
between the cuts. The general structure of the singularities of the correlation function is shown in 
Fig. [7.3]. 



Figure 7.3: A general analytic structure of the static correlation function in the complex three-momentum 
space 9 = |q| (finite temperature). 

The position of the pole is given by the equation whose form is analogous to that of Eq. (|6.37 ), 
namely. 


_^ G5 xJJi(0 , -ml) = 0. (7.31) 

^We note that the Pauli-Villars regularization of the function ^2 tem(Tf^, 0, does not change its analytic 
properties. On the other hand, the 3-dimensional cutoff procedure makes the temperature cuts finite, which leads 
to unphysical oscillations rather than to screening. This fact is the main reason why we have decided to use the 
Pauli-Villars method in the present investigations. 
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The pole can appear only the interval 0 < m* < 2M because for m* = 2M the cut starts. At 
zero temperature, due to the Lorentz invariance, we have m* = mdyn- On the other hand, for 
T > 0 one finds that m* 7 ^ 'nidyn- Formally, the last property follows from the difference between 
the functions I 2 tem(-^^) 0 ) -^2 tem(-^^) 0 ; The latter is the analytic continuation of 

I 2 0 , q^) and has the form 


I. 


2 , mat 


2^ 


(M^0 ,-r) = - 


dpp 


1 


/cvr^ J Ep + 1 


vr — 2 arctg 


2p;j 


(7.32) 


The contribution to the Fourier transform from the pole has, of course, the form ~ exp(—m*r). 
The integral around the vacuum cut can be written as follows 


Xaa(^) = - 


1 




dk k 


'2M 


Xaa(0, -k'^ + ie) - XAAiO, -k^ - ie) 


^—kr 


(7.33) 


For large values of r, the exponential factor in ( 7.331) cuts down the integrand very quickly as a 
function of k. Thus, the remaining factors, which are slowly varying, can be approximated by their 
value at A: ~ 2M. Consequently, in the limit r ^ 00 the contribution of the vacuum cut is also of 
the exponential form, ~ exp(—2Mr). The contribution of the temperature cut is 


+ 00 


Xaa(^) = 


47r^z 


ir 


dq q XAa(0, [q + i5-f) - XAa{ 0, [q + i5+f) 


^iqr 


(7.34) 


where is defined below Eq. (|7.2S ). The imaginary shift of the argument, id±, denotes that 
the function should be calculated just above/below the temperature cut. We shall now discuss 
separately three different cases. The first two concern our analytic results for T = 0 and T —> 00 . 
The third one concerns intermediate temperatures for which we do numerical calculations. 


7.2.2 Results for T = 0 

At zero temperature the temperature cut vanishes so we are left only with one or two contributions. 
In the case of the pseudoscalar channel we have an isolated pole and a cut. This pole, because 
of the Lorentz invariance, coincides with that determining the dynamic mass. For large values 
of r the contribution from the pole, ~ exp(—md^n.-n-?’)) is larger than that from the vacuum cut, 
~ exp(—2Mr). Thus, the pion screening mass at T = 0 is equal to the dynamic one: mscr,n = 
^dyn,TT = 138 MeV. In the case of the scalar channel we do not have an isolated pole and the 
only contribution to the Fourier transform comes from the vacuum cut. In consequence the sigma 
screening mass equals 2M and it is smaller than the dynamic one: mscr,a = 752 MeV and mdyn,a = 
760 MeV. 

We note that the result mscr,<j = 2M is of course an artifact of the NJL model, which lacks 
confinement and of the approximation, which lacks the coupling of the sigma to the two-pion 
continuum. With this effect included, the screening mass in the sigma channel would equal 2mdyn,n- 
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7.2.3 The Case of Extremely High Temperature 


In the limit T —> oo the zeroth-order correlation functions vanish. It means that at sufficiently 
high temperature they are very small and, at the same time, they are a good approximation for 
the full correlation functions. In this case the expression Gs 9^)’ appearing in the definition 

(S.St), can be neglected when compared to unity. Using this fact we dehne the high temperature 
correlation function in space by equation PI 


= j dqqxfAi^,q")e 


2\ iqr 


(7.35) 


Here the integration of the vacuum part XAAvac(^’9^) temperature part 

defines the functions Xaa vac(^) ^AAmat^^^’ I'^spectively. Of course we have 


Xaa(^) = xZ,vac(^) + 


.(0) 


( 0 ) 


(7.36) 


The Fourier transform ( 7.35| ) can be calculated analytically. Using the same method as that 
developed by us in Section 7.1, we find that 


X?l.,vac(0 = ^ E [3i^2(2A,r) + 2A,riXi(2A,r)] 


N, 


N 


i=0 


On the other hand the temperature piece has the form 


(0) / ^ 
App,mat^’’) “ 


2d 1 (92 


dvr^r 


dr 


+ 


d^2 




tern / 


where 


and 


N 


Gf^{r) = -Y,AikiKi{2ki7 

i=0 


N +CO 


=7rTy^^Ai ^ exp{-2r^{21 + l)'^7r‘^T‘^ + kf). 

2=0 / = —OO 


(7.37) 


(7.38) 


(7.39) 


(7.40) 


Using the properties of the modified Bessel functions |54] we can check that in expression ( [7.36 ) 
the vacuum part ( 7.37 ) is exactly canceled by the first term of the temperature part (|7.3^ ), i.e., 
by the term coming from differentiation of the function G]j®™(r). We want to emphasize that this 
cancellation is independent of the number of subtractions N and of the values of the regulating 
masses Ap 

The function G^'^{r) has a simple asymptotic behaviour. For r ^ oo it is enough to consider 
only the leading terms in the series, i.e., these corresponding to I = 0,-1 and i = 0. The terms 

^The situation described in this Subsection resembles the case of the perturbative QCD calculation discussed 
thoroughly in Section 7.1. The basic difference is that now the cutoffs are finite and fixed by fitting physical 
observables at T = 0. 
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for i > 0 are suppressed since At > Aq = M (additionally we can neglect M with respect to T). 
Finally, we obtain the following asymptotic expressions 


and 


Gt®“(r) ~ 27rre-2-^’' 

(7.41) 

(0) ( 2T^Nc _27rTr 

(7.42) 


From Eq. ( 7.42|) we read off the screening mass in the pseudoscalar channel to be 2tiT. The 
calculations in the scalar channel look similarly. We can even immediately guess that the result 
will be the same: At very high temperature the constituent quark mass drops down and is equal 
to the current one. The latter is very small and the difference between the pseudoscalar and scalar 
channel practically disappears. 


Our present considerations showed again how important is the regularization of the temperature 
part. It leads to cancellations between the vacuum and temperature parts and, finally, to the 
asymptotic behaviour of the form ~ exp(—27rTr). If the temperature part were not regularized, 
the asymptotic behaviour of the correlation function would be governed by the lowest regulating 
mass, i.e., by Ai. Such a case could not be accepted from the physical point of view. 



Figure 7.4: The spatial dependence of the static correlation function in the pion channel. The numbers at 
the lines are the values of the screening masses at the corresponding temperatures (given in MeV). 
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7.2.4 Intermediate Temperatures 


For the temperature range 0 < T < 400 MeV we did numerical calculations and read off the 
values of the screening masses from the logarithmic plots representing the correlation functions in 
space. We did the calculations in two ways checking whether the results were the same. The hrst 
method was to calculate the Fourier transform directly from Eq. (3.39). The second method was to 
calculate the contributions from the singularities separately and later to sum them up. In the latter 
case we used Eqs. ( 7.33 ) and ( 7.34 ). Moreover, we numerically calculated the residue corresponding 
to the pole ( [f.31 ). Both methods encounter numerical difficulties. The direct calculation of the 
Eourier transform requires an integration of a slowly converging and oscillating function. On the 
other hand, the contributions from the cuts show large cancellations. Therefore, each contribution 
should be evaluated with a very high accuracy. Let us also note here that the parameters used to 
calculate the screening masses were the same as those used in the case of the dynamic masses, see 
Subsection 6.3.2. 


We had to restrict ourselves to rather small values of r, since the correlation functions de¬ 
crease very rapidly and for large r their numerical evaluation is difficult. In fact, they decrease 
exponentially which suggests that the screening masses can be read off already for small r’s. This 
exponential behaviour is shown in Figs. [7.4] and [7.5], where we plotted the functions rxpp(r) 
and rxss{r). 



r [fm] 


Figure 7.5: Same as Fig. [7.4] for the sigma channel. 

At zero temperature our numerical procedure gives the results which are in good agreement with 
our estimates based on the analytic considerations, see Subsection 7.2.2. The numerically calculated 
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pion screening mass agrees exactly with the analytic result, whereas in the sigma channel we find 
a small discrepancy. This difference is caused by the fact that we do not consider asymptotic 
distances in the numerical calculations. 

The temperature dependence of the screening and dynamic masses is shown in Fig. [7.6]. For 
0 < T < 200 MeV we observe that the screening masses are close to the dynamic ones. Such 
situation is expected, since for small temperatures the Lorentz invariance is only slightly broken. 
Interesting things can be observed for larger temperatures. Both screening and dynamic masses 
exhibit the restoration of chiral symmetry although they substantially differ from eaeh other. In 
the full interval, 0 < T < 400 MeV, the qualitative behaviour of the screening masses and of the 
dynamic ones is similar and resembles the results of the lattice simulations of QCD. 



Figure 7.6: The temperature dependence of the screening (dashed lines) and the dynamic (solid lines) 
masses. 


Our observation that the dynamic masses are different from the screening masses contradicts 
the results of Ref. |13| where it was argued that they are the same even at finite T. In Ref. |j^ the 
chiral limit, m = 0, is considered. Afterwards, the gap equation is used to find that the two masses 
are equal. The point is, however, that in Ref. |]^] the non-trivial solution to the gap equation is 
used. It exists only if T < Tc (in the chiral limit Tc is well defined) and for T > Tc we have to take 
into account just the trivial solution M = m = 0. The trivial solution does not allow us to simplify 
the denominator of the correlation function, what is essential for the proof of the equality of two 
masses in |]J|. Consequently, the arguments of Ref. [|^ do not apply at T > Tc and in this case we 
do not expect that rriscr = mdyn- The problem remains what happens at T < Tc- The structure of 
the poles suggests that mscr = Tndyn, but there exists a temperature cut whose contribution might 
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be not negligible. This situation deserves a separate study. 

Another interesting point is to check whether the screening masses approach the Eletskii-Ioffe 
limit, i. e., if they are equal to 2ttT for large T. Our numerical results for the pseudoscalar channel 
are shown in Fig. [7.7]. One can see that the ratio mscr,TT/T grows with the temperature for 
T > Tc but even at T = 450 MeV it is still smaller than 27r. (In the scalar channel the high 
temperature behaviour is the same as in the pseudoscalar one.) As we have already mentioned this 
type of behaviour is observed in the lattice simulations. In the (pseudo)scalar channel the screening 
masses are smaller than 2ttT, which is in contrast to the (axial)vector channel where they are equal 
to 27rT already for the temperatures slightly exceeding Tc. Our result suggests that the NJL model 
can explain the existence of the non-negligible (residual) interaction in the (pseudo) scalar channel. 
Nevertheless, for drawing more definite conclusions the calculations in the (axial)vector channel 
are required. If the difference between the behaviour of the screening masses in the (pseudo)scalar 
channel and the (axial)vector channel were found, this would act in favor of the model. 


6 


4 

Uh 

O 

1/3 

a ^ 

2 


Figure 7.7: The ratio niscr/T plotted as a function of the temperature. The dashed line corresponds to the 
Eletskii-Ioffe limit: macr = 2'kT. 
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Chapter 8 


Static Meson Fields at Finite Baryon 
Density 


In this Chapter we investigate the in-medium static meson correlation functions at finite baryon 
density. In contrast to the behaviour known from the high temperature calculations, studied by 
us in the previous Chapter, we find that the correlation functions at finite baryon density are not 
screened but oscillate in space |^|. This behaviour is analogous to the Friedel oscillations in a 
degenerate electron gas [ II |; in both cases the oscillations are caused by the existence of the sharp 
Fermi surface. 


In spite of the fact that it lacks confinement, the NJL model successfully describes many aspects 
561. However, a realistic description of nuclear matter at low densities as a 


of hadron structure 

collection of composite nucleons would be very complicated and outside the scope of our approach. 
Hence, we work again in the HF -|- RPA scheme, where the correlations that bind quarks in nucleons 
are neglected and consequently nuclear matter is described as a Fermi gas of quarks. Although 
our model for nuclear matter at low densities is not realistic, we believe that our results on the 
correlation functions at finite baryon density are generic and should be qualitatively the same 
in more realistic models. This conjecture is supported by studies of the effective nucleon-nucleon 
interaction in nuclear matter, which also indicates the presence of the Friedel oscillations 


Already at moderate densities we find new singularities in the correlation function, which prob¬ 
ably are artifacts of the NJL model in the HF -|- RPA approximation. This indicates that this 
approach to the correlation functions breaks down already at moderate densities. Hence, to com¬ 
plete the picture we compute the meson correlation function at higher densities in perturbative 
QCD. Keeping only the leading term, we find oscillations with a period 5r = tt/pf, where pf is 
the Fermi momentum of the quark sea. 


8.1 Low Density Theorem 

Before we analyze the properties of the correlation functions, let us discuss shortly the behaviour 
of the quark condensate at finite baryon density. In this way we shall supplement our discussion 
from Section 5.2. The calculation of the decrease of the condensate with increasing baryon density 
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in the framework of the NJL model is instructive, since there are model independent estimates of 
this quantity, which impose constraints on the in-medium behaviour of the condensate. 

At T = ^ = 0 the NJL model yields the Gell-Mann-0akes-Renner (GOR) relation 

P^^dyn,w = + md){uu + dd) = -2m{qq), (8.1) 

where LV is the pion decay constant. This relation is independent of the regularization scheme. 
One can use it to find the value of TV provided the values of m, m-,^ and (qq) are known. For the 
two-flavour symmetric model, with the values of the parameters the same as in Section 5.2, one 
finds TV = 94 MeV. 



Figure 8.1: Ratio of the in-medium quark condensate to the vacuum value, plotted as a function of baryon 
density. 


Using the Hellmann-Feynman theorem and the GOR relation, one finds |61, 62] to leading order 
in the density 


(gg) ^ 

{qq)o mlF^ 


( 8 . 2 ) 


where STrAf is the pion nucleon sigma commutator, (qq)Q is the vacuum value of the quark conden¬ 
sate, and p is the baryon density. The baryon density is for two flavours given by p = 2p|./37r^, 
where pp, as usual, is the quark Fermi momentum. We compute the pion nucleon sigma term by 
using the relation |62| 



= m 


dMo 

dm ’ 


(8.3) 
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where we assume that the ttN sigma term is simply the sum of the vrg sigma terms, like in the 
naive quark model. For our set of the parameters one gets = 18 MeV, which, using Eq. (8.2), 
implies that the condensate is reduced by 15 % at the saturation density of nuclear matter, po = 
0.17 fm“^. In Fig. [8.1] we show the numerical results for the quark condensate (normalized to 
its vacuum value) as a function of the baryon density (normalized to the saturation density) [|. 
One can see that our calculation agrees with the low-density theorem ( ^.21) , when we use the sigma 
commutator obtained within the model. Since the sigma term is much smaller than the empirical 
value of 45 MeV, the density dependence of the quark condensate is too weak. Nevertheless, the 
fact that Eq. (8.2) is satisfied, shows that the calculation is consistent. 


8.2 Priedel Oscillations in the NJL Model 


Let us come back to the discussion of the correlation functions. Our analysis of the static me¬ 
son correlation functions at finite baryon density is based on the two-flavour symmetric model. 
Moreover, we use again the same values of the parameters as in Sections 5.2, 6.3 and 7.2. Thus 
our results for zero baryon density agree with the results of Section 7.2 for the case T = 0. Our 
starting point are Eqs. (|6.39 ) an d (|6.36|) , where the zeroth-order correlation functions are given by 
formulae (|^), (^), (|6.16f) and (|6.22|) . Thus, the full correlation function (|6.36| ) can be expressed 
in terms of the functions mat(-^^); -^ 2 ,vac(~<?^) and that 

these four functions have been already calculated. In the case of the Pauli-Villars regularization 
scheme, which is adopted here, we can use Eqs. ( 5.24 ), ( 5.25 ), ( 6.31 ) and ( 5.33|) , respectively. 


Upper vacuum cut ' 

Imq 

2M 

-2Pf / 

2pp 


Re q 

matter cuts 



i 2M 



lower vacuum cut / 



Figure 8.2; A general analytic structure of the static correlation function in the complex three-momentum 
space g = |q| (finite baryon density). 

^We note that Fig. [8.1] corresponds exactly to Fig. [5.7] showing dependence of M on p. 
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The Fourier transform (|6.39| ), defining the static correlation function in space, can be again 
written as a sum of a few contributions connected with the appearance of the singnlarities in the 
complex momentum plane. At finite baryon density, the strnctnre of the singnlarities is very similar 
to that at finite temperature. It is shown in Fig. [8.2]. 

The basic difference between the finite temperature and finite density cases is that the matter 
cut has a finite range in q (the temperatnre cnt, shown in Fig. [7.3], is infinite). To understand 
this behavionr, we can analyze the function I 2 ysxq^\^{M ‘^in the same way as the function 
-^2 tern (-^^5 0) 9^) analyzed in Subsection 7.2.1. First, we write 


where 

pf 

(±) Nc r dpp 2p- z±ie 

2 ,mat J k ^-2p-z±ie' 

0 


(8.4) 

(8.5) 


The form of expression on the RHS of Eq. ( |8.5|) shows that the functions (z) have logarithmic 

cuts parallel to the real axis and stretching from — 2pi7’ ± ze to +2pp ± ie. The finite range of the 
cuts is caused by the finite range of momenta available in the integral ( ^.5|) . 

Let us now discuss the results of the finite density calculations. In the case of the pseudoscalar 
channel, we did the numerical calculations in the energy range Mq = 376 MeV < p < 415 MeV. In 
Fig. [8.3] we show the corresponding correlation for p = 376,400 and 410 MeV. In contrast to the 
exponential decay found in vacuum (solid line), the correlation function at finite density (dashed 
and dotted lines) oscillates, with a power-law decay of the amplitnde. Its period decreases with the 
increasing value of the chemical potential. From the physical point of view, such oscillations are 
caused by the existence of a sharp Fermi surface. 


The oscillatory behaviour of the correlation functions at finite density is well known for the 
non-relativistic degenerate electron gas, where the phenomenon is called the Friedel oscillations 
[M|. These long range oscillations lead to many interesting phenomena like, e.g, broadening of the 
nuclear magnetic resonance lines. A characteristic feature of the Friedel oscillations (at very large 
distances) is their period 5r = tt/pf- In onr case the situation is similar; for very large values 
of r, the contribution from the matter cut dominates over the other contributions. Since the cnt 
extends over a finite range in q (jgj < 2pf) the correlation function oscillates with the period Tr/pp 
(at r oo). The numerical study of the correlation functions at very large distances is difficult 
because the amplitudes decrease with r. Consequently, we restrict the numerical calculations to 
the interval 0.5 fm < r < 6 fm. 


So far we have discussed the results for p < 415 MeV. For larger values of the chemical poten¬ 
tial (corresponding ronghly to p > ^po) the pseudoscalar correlation fnnction acquires additional 
singnlarity on the real axis for q Pi 2pF- In the scalar channel, a similar singnlarity appears at even 
smaller densities (p > |po)- These singularities indicate instabilities of the ground state, which 
lead the system to states of lower energy. However, most likely the singularities are artifacts due to 
regularization procedure and do not correspond to physical instabilities. In order to check this point 
we have redone the calcnlation nsing a different regnlarization scheme with a 3-dimensional cutoff. 
Using the parameters obtained in Ref. ]^], we find no singularities on the real axis. This shows 
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Figure 8.3; Correlation function in the pseudoscalar channel plotted for three different values of the chemical 
potential: ji = 376 MeV (solid line), /i = 400 MeV (dashed line), and ji = 410 MeV (dotted line). 


that these singularities are unphysical, since their presence depends on details in the formulation 
of the model. Thus, our regularization scheme can be used only at low densities. 

On the other hand, the 3-dimensional cutoff regularization is unsatisfactory because it explicitly 
breaks Lorentz invariance. This obscures the relation between the correlation functions in time-like 
and space-like regions. Moreover, in the 3-dimensional cutoff scheme, one does not recover the 
well known screening of the correlation function at finite temperatures. Instead, the correlation 
function oscillates, much like the Friedel oscillations at T = 0, with a period vr/A, where A is the 
momentum cutoff. We stress that this behaviour is an artifact, due to the finite range of momenta 
(|g| < A) available in the Fourier transform. 

Consequently, we stick to the covariant Pauli-Villars method, in spite of its shortcomings. 
However, since the results are inconclusive at densities p ~ po and higher, we restrict the calculations 
to small densities p < \po- The physics at very large densities will be addressed in the following 
Section. 
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8.3 Oscillations in Perturbative QCD 


In a straightforward way, the formalism developed so far can be used to study finite-density corre¬ 
lations functions in perturbative QCD. Similarly to the high temperature case discussed in Section 
7.1, the leading term in the meson correlation function is the lowest order quark loop diagram (as 
in Fig. [6.1]). Thus, for the pseudoscalar channel we can write 


X{r) = Xvac(r’) +Xmat(0, 


where Xvac(^) is given by Eq. (7.7), and the matter piece has the form 


( 8 . 6 ) 


Xmat(^) 


Nc 

dvr^r \ 


2 d 1 92 





(8.7) 


Eq. ( |8.7|) is the analog of Eq. ([7.17|) for the case of finite density. One can check that = 

G']'em( 7 -), Similar cancellations to those described in Section 7.1 take place. Thus, finally only 
the function gives a non-vanishing contribution to the correlation function. A simple 

calculation yields 


(5^^*(r) 


d r°° 

dr Jpp 


dp 


cos(2pr). 


( 8 . 8 ) 


In the limit M —> 0 we find (5™^*(r) = cos{2ppr)/2r, which implies that the correlation function 
x{r) is not screened but oscillates in space with a period 5r = njpF. For massless quarks the 
behaviour in the scalar channel is identical. 
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Chapter 9 


Beyond the Hartree-Fock 
Approximation 


Calculations based on the NJL model show that the value of the quark condensate is strongly 
modified by the in-medium effects. It eventually goes to zero for T, // —> oo. The interesting 
question is whether the changes of the condensate (as described by the model) exhibit some universal 
features. At low density, we have seen (Section 8.1) that the NJL calculations agree with the model 
independent estimates. On the other hand, it turns out that the change of the condensate at small 
temperature is not compatible with the general requirements of the chiral perturbation theory: 
In the standard approach to the model (HF -|- RPA), the condensate remains completely flat at 
small T, whereas the chiral perturbation theory indicates that the absolute value of the condensate 
should decrease linearly with T^. 

As we shall see in this Chapter, the discrepancy between the NJL model and the chiral perturba¬ 
tion theory is an effect caused by the HF approximation. The standard approach is not sufficient to 
describe properly the thermodynamic properties of hadronic matter at small temperatures, which 
should be dominated by the lightest hadrons, i.e., pions. In order to describe properly the low tem¬ 
perature change of the condensate in the NJL model, one has to go beyond the HF approximation, 
including the so-called meson loops. In this Chapter, we outline the characteristic features of this 
approach. 


9.1 Effective Action 


Our present considerations will be based on the two-flavour standard model with scalar-isoscalar 
and pseudoscalar-isovector interactions. The corresponding Lagrangian is given by Eq. (3.16). In 
this Section, we shall apply the formalism of the effective action to this Lagrangian. Details of 
this procedure have been given in |64]. In the following, we shall describe only the basic ingredients 
of this method. 

At hrst we introduce the meson mean helds = (<I> 0 ) ^)) where $0 and $ correspond to the 
pion and the sigma, respectively. At the quark-loop level the effective action is 
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/(<!>) = J d^x — ^Tr \n{D^D), 


(9.1) 


where l/o^ = 2G and D is the Dirac operator, i.e., D = 9,- — io: • V + /3<ho + iP^ 5 T • For the 
moment, we consider the case T = 0 and work in Euclidean space-time (r, x). The generalization 
of our considerations to the case T > 0 will be described in the next Section. In Eq. o we have 
replaced the usual Tr InD term with ^Tr In(DlD), which is exact for SU(2). 


Meson loops bring an additional term to the effective action 


r($) = /($) + -Tr ln(iF-^), 

where K is the inverse meson propagator matrix defined as 


(9.2) 


Kab y) = 


52/($) 


5^aix)d^b{y) 


(9.3) 


In Eqs. (|9.lD and (^) Tr denotes the full trace, including functional space, isospin, and in addition 
colour and spinor trace for quarks. In the A^c-counting scheme, the quark loop term /($) is the 
leading contribution of order 0{Nc), and the meson loop term ^Tr ln(A'“^) is of order 0(1). Thus 
the one-meson-loop contributions give the first correction to the leading-results. 

Using standard methods. Green’s functions can be obtained from Eq. (9.2) via differentiation 
with respect to mean meson fields. Of particular importance is the one-point function, which gives 
the expectation value of the sigma field. The condition 


hr(^) 

(5$o(a;) 


<I>oU)=S' 


= a^{S — m) — ^Tr 

+2^" ( 


'(0.0)-.^) 


h$o(a;) J 
= 0 


$oU)=S' 


(9.4) 


4-oU)=S 


yields the equation for the vacuum expectation value of <l>o, which we denote by S. As shown in 
Ref. [^], introducing 


K.{S,ql) = 
KAS,ql) = 


—h,vaciS , —qE) {qE + 4 : 5 " ) -|- 


a m 


-1 


—h,vac{S‘^, —qE) qlk + 


am 


-1 


(9.5) 


and retaining terms up to order 0{N^), Eq. (p.4D can be written in the form 
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(S - m) - S h,yaciS^) 

_^ 

16 tt^ 

+ 3 


J I 2/2,vac(*S'^) 0) + (^/2,vac(*S'^5 — 

2/2,vac(<S'^;0) + -^j^l2,vaciS‘^, —Qe) Qe ^7T{S,q%)'^ = 0. 


Ka{S, q^) 


(9.6) 


The functions h,vac{S‘^) and l 2 ,vac{S^,—Q^) have been defined by expressions (^.21 ) and (|6.7|) . 
One can easily check that the leading-iVc term in (|9.6D coincides with the gap equation (5.28). 


It is important to observe that the one-meson-loop gap equation ( |9.6| ) requires introduction of 
an extra regulator for meson momenta. The regularization of the quark loop only is not sufficient 
to make meson loops finite, since at large one has / 2 ,vac('S'^, —gl;) ~ 1/^1;• Consequently, 
the structure of divergences is the same as, e.g., in the linear a model, in particular Eq. (^) is 
quadratically divergent without a cutoff in the qE integration. The need for an additional cutoff is 
not surprising, since a non-renormalizable theory requires new cutoffs at next levels of a perturbative 
expansions. Therefore we have to regularize the divergent integral over d^qE- In Ref. this was 
achieved by the substitution 


A2 



0 


where was the four-dimensional Euclidean meson momentum cutoff. At the same time, the 
quark loops were regularized with the Schwinger proper-time method. In the present study we 
employ the 3-dimensional cutoff procedure for the meson loops, i.e., we make the replacement 


J d'^qE —> dvr J dqi J dqq^, 


(9.8) 


where qE = (^ 4 , q) and q = |q|. The form (9.8) is convenient for the implementation of the boundary 
conditions satisfied by temperature Green’s functions. Similarly to Ref. |^, we regularize the 
quark loops (i.e., the functions /i,vac('S'^) and / 2 ,vac(>S'^, —<?|;)) using the Schwinger proper-time 
method. 


If chiral symmetry is broken, then Eq. 
sate and S are related by the formula 


has a nontrivial solution for S. The quark conden- 


{qq) = -a^{S - m) , 


which follows immediately from the fact that {qq) = 5T{^)/5m and Eq. (p.2p. 


(9.9) 
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9.2 Finite Temperature 


For calculations at T > 0 we shall again adopt the imaginary time formalism. In this situation, the 
finite-temperature analog of Eq. (^) can be written in the following form 


ST 


A(, 


a^{S-m)-S If’ T) - I diq 


21 . 


R, SPT, 


( 52 ,0,0, T) + ^ q, T) [{eI)^ + + 45^]) 


K^{S,n,q,T) 


+ 3 


2jR,SPT(^2^0,0,r) + 


d 

d^' 


tR, SPT, q2 


{S\n,q,T) (E^)2+q 


K^{S,n,q,T)\ =f) . 


(9.10) 


where Kfj{S,n, q,T) and KT^{S,n,q,T) are the generalized expression for the inverse meson prop¬ 
agators 


Ka{S,n,q,T) 

K^{S,n,q,T) 


f^^^{S\n,q,T) {Eif + c^^ + AS^ 

2 


+ a^m 


{ 

|-/^’^^^(S’^,re,g,r) +a^m/S^ 


/s}~\ 


(9.11) 


In Eqs. ( 9.10 ) and ( 9.11J ) E^ = 27rnT, and the functions ‘^^^(S'^, T) and ,n,q,T) are 

given by expressions ( 5.32| ) and (6.35) — since in this Chapter the Schwinger proper-time method is 
used to make the quark loops finite, we make the reference to the expressions giving the regularized 
form of/i(5^,r) and 12(5'^) Pj^T)- 


9.3 Low-Temperature Expausiou iu the Chiral Limit 


Before presenting our numerical results for {qq)T let us consider the low-temperature expansion. 


As shown by Gasser and Leutwyler |^], in the exact chiral limit the low-temperature expansion of 
the quark condensate has the form 


fe)T = fe)o (1 - ^ - 3^]^ + ■■■) ■ (912) 

Eirst, let us do the Nc counting in this formula. Since ~ subsequent terms in the 

expansion are suppressed by l/Nc- Since our one-meson-loop calculation accounts for first sub¬ 
leading effects in the 1/A"c expansion, we can hope for reproducing only the term in Eq. ( |9.12 ). 
Eurther terms would require more loops. 


Using techniques described in Chapter 4, the sum over the bosonic Matsubara frequencies in 
Eq. (^) can be converted to a contour integral in the complex energy plane. By deforming 
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this contour we collect all contributions from the singularities of the integrand, weighted with the 
thermal Bose distribution. At low temperatures, the dominant contribution comes from the lowest 
lying pion pole, and other singularities are negligible. Thus, the third term in (9.10) becomes 


n g 


dq 






1 + 


,glT _ 1 


(9.13) 


In Eq. (9.1C) we have approximated the function l 2 {S‘^,n,q,T), appearing in the pion propagator, 
by its value at n = g = 0. For sufficiently large cutoff A^, the integral over the thermal distribution 
function i n (|9.13| ) can be expressed by the Riemann zeta function ("(2) = Thus, the final 

result for (|9.13|) is 3 A^/ 47 r^ + r^/ 2 . Inserting the above result into the gap equation (|9.10D we find, 
with m = 0 , the following equality: 


h{S,T) = a^ -h{S^,T) + 


3A2 1, 


+ -T^ = 0. 


’ 2 

Eq. (|9.14) defines implicitly the function S{T), which satisfies the equation 


(9.14) 


dS dh{S,T)/dT‘^ 
^ “ dh{S,T)/dS 


, dh{S^T) 

^ dS 


-1 


(9.15) 


Here we have neglected the term dIi{S‘^,T)/dT‘^, since it is exponentially suppressed by the 
factor exp(—5/T). Furthermore, the RHS of ( |9.15| ) can be rewritten using the relations |64] 
dIi{S‘^,0)/dS = 4 S'/ 2 (S'^, 0,0, 0) and — 12 ( 5 ^, 0,0, 0) =F^/S'^, where is the leading-A'^c piece 
of the pion decay constant. Collecting these equalities we arrive at dS/dT'^ = — 5/(8F^), which 
finally gives 

r2 “ 


S{T) = ,9(0) 


1 - 


8F, 


(9.16) 


Proportionality (|9.9| ) implies that the above expression coincides (in the large Nc limit, where 
Ett —> En-) with Eq. (9.12). Hence our method is consistent with a basic requirement of chiral 
symmetry at the one-meson-loop level. 


9.4 Results of Numerical Calculations 

In the exact chiral limit the model has 3 parameters: a, A, and A;,. Since the objective of this 
calculation is to explore qualitatively the role of meson loops rather than to come out with accurate 
fits to many observables, we present the calculations for the arbitrary choice A;,/A = We have 
also done calculations for other values of A;,, with similar conclusions as long as it is not too small. 
The remaining two parameters are fixed by reproducing the physical value of E^- = 93 MeV and a 
chosen value for {qq)o. For the case of m 7 ^ 0 we have an extra parameter, m, which is fitted by 
requiring that the pion has its physical mass. We compare results with meson loops to results with 
the quark loop only (A^ = 0). Parameters for the two calculations are always adjusted in such a 
way, that the values of E^-, {qq)o, and are the same. 
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The calculation of with meson loops, although straightforward, is rather tedious, so we do 
not present it here. The method has been presented in detail in Ref. [^, |^. The only difference 
in our calculation is that the three-momentum cutoff ( |9.8D rather than the four-momentum cutoff 
of Ref. iQ is used. 



cr 

cr 



[GeV^] 


Figure 9.1: Dependence of the quark condensate on in the chiral limit tUtt = 0. The curves correspond 
to the calculation with meson loops (solid line), with quark loops only (dashed line), and the lowest-order 
chiral expansion (dotted line). The parameters for the solid line and dashed line are adjusted in such a way 
that Ftt = 93 MeV and (^< 7)0 = —(184 MeV)^. For the solid line a = 175 MeV, A = 723 MeV, and Af, = lA, 
whereas for the dashed line a = 201 MeV, A = 682 MeV, and At = 0. 

Figure [9.1] shows the dependence of (qq) on T^. The solid line represents the case with meson 
loops. We note that at low temperatures the curve has a finite slope, as requested by Eq. (ED- 
The slope is close to the leading-order Gasser-Leutwyler result (dotted curve). As explained earlier, 
the slopes would overlap in the large-Ac limit. This behaviour is radically different from the case 
with quark loops only (dashed curve). In this case at low temperatures {qq)T — {qq)o ~ 
where M is the mass of the constituent quark. All derivatives of this function vanish at T = 0, and 
{qq) is flat at the origin. We can also see from the figure that the fall-off of the condensate is faster 
when the meson loops are included. In fact, for the parameters of Fig. [9.1] we have an interesting 
phenomenon. At T = 162 MeV the condensate abruptly jumps to 0. There is a first-order phase 
transition, with a latent heat necessary to melt the quark condensate. Such a behaviour is not 
present in the case of calculations without meson loops [^, ^]. We note that with meson loops 
present the chiral restoration temperature is 162 MeV, i.e., about 10% less than 176 MeV of the 
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quark-loop-only case. 

Figure [9.2] shows the same study, but for the physical value of m.^-. We note that now l^q) 
(solid line) is also flat at the origin, since the pion is no more massless, and at low T we have 
{qq)t — {qq)o ~ Nevertheless, the region of this flatness is small, and at intermediate 

temperatures the curve remains close to the Gasser-Leutwyler expansion. We note again that 
meson loops considerably speed up the melting of the condensate compared to the case of quark 
loops only. However, there is no first-order phase transition such as in Fig. [9.1]. Instead, we 
observe a smooth cross-over typical for the case of m 7 ^ 0 . 


The faster change of the quark condensate in our study is not surprising. It is caused by the 
presence of light pions which are known to play a dominant role at low-temperatures. See, for 
example. Ref. [16| where the bulk thermodynamic quantities studied within the NJL model are 
determined by the pion quantum numbers. In this case, however, the decrease of the condensate 
does not agree with the requirements of the chiral perturbation theory, since the scalar-density 
equation is considered at the quark-loop level only and the slope of {qq) vs. vanishes at = 0 . 



Figure 9.2: Same as Fig. [9.1] for tUtt = 139 MeV, = 93 MeV, and {qq)o = —(174 MeV)^. For 
the solid line a = 164 MeV, A = 678 MeV, Af, = |A, and m = 15 MeV, whereas for the dashed 
line a = 175 MeV, A = 645 MeV, At = 0, and m = 15 MeV. 
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Chapter 10 


Real-Time Formalism 


In this Chapter we give a short introduction to the real-time formalism, presenting the most useful 
concepts and definitions. For our further considerations we do not need to know this method in 
more detail, nevertheless for a better understanding it is useful to discuss in a separate place the 
way in which the non-equilibrium Green’s functions are introduced. For more information about 
the real-time formalism we refer the reader to ||67||. 


Since we are interested in the non-equilibrium evolution of a system, we have to calculate the 
averages of the operators in arbitrary states. Thus, in analogy to (|4.9|) or (4.10), we define 


iG(ti,xi;t2,X2) = (T ['0(ti,xi)V^(t2,X2)]), 


( 10 . 1 ) 


where the symbol (...) denotes an expectation value with respect to the arbitrary initial state. The 
latter is specified at t = to by the density operator p, i.e., we have (...) = Tr(p...)/Tr(/5). 


The field operators appearing in (|10.1 ) are in the Heisenberg picture. In the situation when the 
perturbation expansion is appropriate, the following representation is frequently used 


0„{t) = U{to,t)0,{t)U{t,to), (10.2) 

where 0^{t) and Oj{t) are the operators in the Heisenberg and interaction pictures, respectively. 
The quantity U(t,to) is the evolution operator. For t > to we have 


and 


t t 

U{t,to) = J dh ...J dtnH]{h)...H]{tn) 

. to to 


n=0 


= T 


exp f dt'H^{t') 

\ to / 


U{to,t)=f 


to 


exp I —i / dt H^{t) 


(10.3) 

(10.4) 
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where H^{t) is the interaction Hamiltonian in the interaction picture, and T (T) is the chronological- 
ordering (antichronological-ordering) operator. 


Let us now consider the expectation value of the operator 0^(t) in the ground state of an 
interacting many-body system, denoted later by |0). Using a group property of the evolution 
operators U as well as the Gell-Mann theorem [41| we find 


(0|O,(t)|0) 


($|r 

exp / dt'HAt')^ Oj{t) 

|c^) 

($1 

T 

exp / dt'Hjit')^ 

[ch) 


(10.5) 


where | ) is the ground state of the noninteracting system. Eq. ( 10.5 ) can be generalized to the 

case of the Green’s functions, namely, we can write 


iG(ti,xi;t2,X2) = (0| T ['!/’(ti,xi)V^(t2,X2)] |0) 


($|r 

exp z 

OO \ _ 

I dt'i7|(t') V’7(U,Xi)V’^(t2,X2) 

-OO / 

|<I>) 

($|r 

exp / dt'H}{t')^ 

1$) 


Application of the Wick theorem to (|10.6| ) leads to the usual Feynman rules. In particular, the 
denominator cancels the disconnected diagrams appearing in the numerator. 


The method described above cannot be applied to the average in an arbitrary state. The 
reason is that such a state, in general, cannot be simply expressed in terms of the ground state 
of the noninteracting particles. Nevertheless, substituting (|10.3|) and ( 10.4 ) into (|10.2|) we obtain 
expression 


{Onit)) = {U{to,t)Oj{t)U{t,to)) 

to 


= {T 


exp I —i / dt'H^{t') 


OAt)T 


exp I —i J 
*0 


(10.7) 


which can be put into a form analogous to Eq. (10.5). This can be achieved by introduction of an 
integral contour JC running along the time axis and the associated operator ordering along this 


contour, see Fig. [10.1]. Using this procedure Eq. (10.7) can be rewritten as follows 


{0,{t)) = {T^ 


exp [ —i / dt'H^{t') j OjA) 


( 10 . 8 ) 
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Figure 10.1: Contour K. used to calculate the operator expectation value 0„{t). 


Generalizing the last result, we introduce the Green’s function on the contour 


tG^(ti,xi;t2,X2) = {T^ [V'(ti,xi)V^(t2,X2)]), (10.9) 

which has the following representation 


iG,^{ti,xi-,t2,X2) = (T^ 


exp -i / V’/(G,xi)V'^(t 2 ,X 2 ) 


( 10 . 10 ) 


Here the contour runs from to to the largest argument of the Green’s function and then goes back 
to to- The formal resemblance to (10.6) is recovered since the disconnected diagrams vanish 


1 = (T^ 





( 10 . 11 ) 


Assuming that the initial state allows for the Wick decomposition, G^(ti,xi;t 2 ;X 2 ) can be calcu¬ 
lated using the standard technique for Feynman diagrams. 

Very often it is convenient to divide the contour into the two branches — the upper one (stretch¬ 
ing from smaller to larger values of t) and the lower one (stretching from larger to smaller values of 
t). According to the time ordering T^, the time arguments along the upper branch are regarded as 
earlier than those situated along the lower branch. Depending on the position of its time arguments 
the contour Green’s function fits to one of the following categories: 
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,xi;t2 ,X2) = iG(ti,xi;t2,X2), (10.12) 

iG^(tj“,xi;t^,X2) = zG<(ti,xi;t2,X2) = -{'4>{t2,X2)'ip{ti,xi)), (10.13) 

iG^(t)^,xi;t^,X2) = zG>(ti,xi;t2,X2) = (V'(ti, xi)Vi(t2, X 2 )), (10.14) 

iG^{tf,xi;t^,X2) = iG(ti,xi;t2,X2). (10.15) 


The Green’s function (tj",xi; )X 2 ) is the usual causal (Feynman) function. On the other 
hand, (t)'', xi; X 2 ) is anti-causal. In addition, there are two extra functions: Gi^{tf ,:x.i;t 2 ,^ 2 ) 
and G^(t)'', xi; t 2 )X 2 ). Nevertheless, one can check that only two out of these four functions are 
independent. 



Figure 10.2; The second time argument of the function G<(ti, xi; ^ 2 , X 2 ) is later than the first one. Thus, 
the time ordering along the contour always changes the initial ordering of the field operators (in addition 
there is a change in sign for fermionic operators). 

The fact that at least two independent functions are required for description of non-equilibrium 
situations has a simple physical interpretation: besides the Green’s function describing the spectrum 
of the excitations in a medium, we need another function describing the evolution of the medium 
itself. In the classical limit such evolution is described by the kinetic equation fulfilled by the phase- 
space distribution functions. In the next Chapter we shall study how the information contained in 
G'^(ti,xi; t 2 ,X 2 ) allows us to derive the transport equation for the quark distribution function. 
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Chapter 11 


Mean-Field Transport Theory 


In this Chapter, we shall study how the concept of chiral symmetry can be explicitly included into 
the quantum and classical transport theory based on the NJL Lagrangian. In our investigations, 
we shall restrict ourselves to the mean-field approximation, since the interesting phenomena (like, 
e.g., spontaneous symmetry breaking) already take place at this level. Our approach is based 
on the spinor decomposition of the Wigner function and follows the treatment of and |69]. 
However, in contrast to |^] we do not neglect the spin degrees of freedom and take into account 
all the components of the spinor decomposition. On the other hand, we differ from |^, where the 


spin dynamics is fully discussed, in that we study a system that is governed by a different type of 
interaction. 


A transport theory for the NJL model has been initially formulated by Zhang and Wilets |]L], 
where the closed-time-path formalism (introduced in Chapter 10) combined with the effective-action 
method has been used to derive the kinetic equations. However, Ref. |^] does not discuss how 
the symmetry of the underlying theory (i.e., the chiral symmetry of the NJL Lagrangian) should 
reflect itself in the transport equations. Trying to clarify this point, we derive the general form of 
the transport equations (in the mean-field approximation) which are chirally invariant. We discuss 
the properties of these equations and then study their classical limit by performing an expansion 
in h for all the functions appearing in the spinor decomposition. Throughout our approach, we 
investigate and check the consistency of our equations under the assumptions of chiral symmetry. 
In particular, this implies that the axial current should always be conserved. For the spin evolution 
equation, as will be seen, this turns out to be a non-trivial requirement. 


11.1 Wigner Function and its Spinor Decomposition 

In general, for a non-equilibrium transport theory, one requires two independent Green’s functions 
to describe the system. However, it turns out that in the mean-field approximation (which neglects 
the collisions between the particles) only one function is sufficient. We thus choose the Green’s 
function 
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Gapix^y) = {lpp{y)^a{x)) 


( 11 . 1 ) 


to be the fundamental quantity in our description (for sake of simplicity we disregard here the 
imaginary unit appearing in the “canonical” definition (|10.13|)). 


In the following we shall restrict our considerations to the one-flavour model defined in Section 
3.1. Moreover, we shall work in the strict chiral limit setting m = 0. Only in the end of this 
Chapter (Section 11.7) the case m ^ 0 will be taken into consideration. One can easily notice that 
the Green’s function (|11.1 ) satisfies the same equation as the field V’(a^) does (compare Eq. ( 3.21 )) 


[i^- a{x) - i 75 vr(x)] G<(x,y) = 0. (11.2) 

The mean fields (t{x) and 7r(x) are determined by the Green’s function G''{x,y) through the 
relations 


cr{x) = —2G Sp G'^{x,x), Tr(x) = — 2G Sp i'j^G^ {x, x), 


(11.3) 


where the trace Sp is taken over the spinor indices. 

For the description of non-uniform systems, it is convenient to introduce Wigner transforms. 
This allows one to make an easier physical interpretation of the components of the Green’s function 
and facilitates performing the classical limit. In the case of the Green’s function G^{x,y), its 
Wigner transform is obtained by introducing the center-of-mass coordinate X = [x + y)/2, the 
relative coordinate u = x — y, and by taking the Fourier transform with respect to u. Following 
the notation of |7C], we define 


W^piX,p) = ^l d^ue^nP-G<f, 


X-^ 

2 


(11.4) 


In Eq. ( 11.4 ) we have explicitly incorporated the Planck constant h, since later on we wish to 
investigate the classical approximation. 

The Wigner transforms of the derivative of a two-point function, df{x,y)/dx^, and of the 
product of a one-point function with the two-point one, f{x)g{x,y), are given by the expressions 


df{x,y) 

dxP 

f{x)g{x,y) 


? 7 ? 

f{X)g{X,p) - -d^f{X)d^g{X,p), 


(11.5) 

( 11 . 6 ) 


with the notation = d/dXP,dP = d/dp^, and where in the second relation ( 11.6 ), only the 
derivatives of first order in dp and dx are kept. This approximation is appropriate only for the 
systems which are weakly inhomogeneous (see the discussion of this point in 


Using Eq. ( 11.2 ) and the properties of the Wigner transform ( 11.5 ) and ( 11. 6| ), we arrive at the 
following equation for the Wigner function 
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ih, 


ih 


pM + Q^^ + —df,a{X)d^ - i-f5Tr{X) - -j5d^7r{X)d^ 


W{X,p) = 0. (11.7) 


Eq. ( 11.7 ) together with Eqs. ( 11.3 ) and (|11.4 ) form the system of the coupled equations that 
we are going to study in detail in this Chapter. We note that there is no collision term in Eq. 
(|11.7|). This is a consequence of our mean-field approximation. We observe that under the chiral 


transformation (3.24), the Wigner function changes according to the prescription 


W = ex.p{-ij5^) W exp(-i75^). 


( 11 . 8 ) 


Hence, one can check that Eq. ( 11.7 ) is invariant under transformations ( 3.25 ), ( 3.26 ) and ( 11.8 ) 
and thus that the first order in the derivative expansion of the Green’s function is symmetry 
preserving. 


The Wigner function Wai 3 {X,p) is a 4 by 4 matrix in spinor indices and can be expanded in 
terms of 16 independent generators of the Clifford algebra. The conventional basis consists of 1, 
*75)7^)7^75) s-iid which we denote in what follows by Tj. In this basis, the Wigner function 
has the form 


W = x + 


(11.9) 


The properties = W and = Tj (the first one results from ( 11.4 ), whereas the 

second one is a simple consequence of the definition of Tj) indicate that the coefficients in the 
spinor decomposition X{X,p),V{X,p),Vf^{X,p),A^{X,p) and Sf^i^{X,p) are real functions. We 
note that the spinor decomposition technique has been already used in the formulation of the 
transport theory for QED ^ and for QHD (so-called quantum hadrodynamics) []^ . 

The vector and axial currents, see Eqs. (|3.60 and (|3.3|) , are related to the functions V^{X,p) 
and A^{X,p) by the simple relations 


and 


Vf^{X) = aJ 

At^{X) = -^j 


d'^p 

d^p 


vnx,p) 

(11.10) 

■A^{X,p). 

(11.11) 


Therefore V^{X,p) and A^{X,p) can be interpreted as the momentum densities of these two cur¬ 
rents. For the physical interpretation of the other components we refer the reader to Refs, f^ , 71]. 
Under the chiral transformation (|11. j ) the coefficients of the spinor decomposition (|11.9| ) change 
as follows 


JT ^ 7^' = .T cos X + T’ sin X) 

V R'= — JFsinx + T* cosX) 

V;. ^ v; = v^, (11.12) 
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A^ 

A! 

5. 

II 


S^u - 

- 5 ;, 

= Spy cos x +Spy siny. 


1 

- 5 ;, 

= -Spy sin X + Spy cos x 


where is the dual tensor to , namely 


11.2 Kinetic Equations 


Substituting expression ( ll.gj ) into Eq. ( 11.7 ) gives the following system of the coupled equations 
for the coefficients of the decomposition (11.9) 


- a:F + ttV = - 
-iK^Af,- cjV = - 

+ iK’'Syfj, - aV^ + iirA^ = - 

iK^^V - - aAf^ + iTrV^^ = - 

K^Afi +aS>^'' = 


(d^ad^T - d^7Td;v) , (11.14) 

{du(jd;v + d^,TTd;j ^), (11.15) 

(duCrdpV^ - idu-irdpAfj^ , (11.16) 

y [d.ad^Af^ - id, 7 Td!;v ^), (11.17) 

? ?? ~ 
^id^ad'^Sf^‘'-d^Trd^Sf^‘'). (11.18) 


In abbreviation we have introduced here the notation . Eqs. ( 11.14 ) - ( 11.18 ) 

should be supplemented by the formulae that determine the mean fields. They follow from (11.3) 
and have the form 


<r(.Y) = -2 g/^ Sp W(X.p) = -8G I -§^,HX.P) 




and 


7r(X) = -2gI^ Sp i^5W{X,p) = 8g| -^V{X,p). 


(27 

(Ap 


(11.19) 

( 11 . 20 ) 


Using the transformation properties ( p.25[ ), ( |3.26| ) and ( |11.12| ) we can check that this system 
of equations is chirally invariant. Strictly speaking, each of Eqs. ( |11.14| ), ( |11.15| ) and ( 11.18 ) is 
separately a chirally invariant equation. On the other hand, Eqs. ( 11.16| ) and ( |11.17] ), after a chiral 
transformation, form linear combinations of each other. These combinations imply, however, that 
each of their components must vanish separately. In consequence, we can treat the system of the 
two equations ( 11.16| ) and ( 11. IT)) as a chirally invariant expression. In the analogous way, one can 
check that the system of Eqs. (|11.19|) and (11.20|) is chirally invariant. 
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For further analysis, it is convenient to discuss separately the real and imaginary parts of Eqs. 


( 11.14 ) - ( 11.18 ). The real parts give 


p>^Vf,-aT + ttV = 0, (11.21) 

^d^A^-aV-irT = 0, (11.22) 

Pf.T aVf, = (11.23) 

_ (11.24) 

_ ^ (9MV^ _ a^V^) _ = 0. (11.25) 


One can notice that integrating of Eq. (11.22|) over the momentum gives 


^ / ■4:^9f,A>^{X,p) = 0, 


(27r 


(11.26) 


where the definitions (11.19|) and (11. 20]) have been used. This equation is nothing other than 


the statement that the axial current should be conserved (^). We thus see that after making the 
gradient expansion, this conservation law is still included in the transport equations. The imaginary 
parts of Eqs. ( 11.14 ) - ( 11. ij ) yield 



(d^adpX - duTTdpV'j , 

(11.27) 

II 

^ [duCrdpV + dy-ndpT^ , 

(11.28) 

+ p'^Si/fi + TrAfj, — 

-^dyad'^Vf,, 

(11.29) 

p^^V + ^TV^^ = 

-hyad;A>^, 

(11.30) 



(11.31) 


Although we have neglected the higher order gradients, Eqs. ( 11.21 ) - ( 11.31 ) are still quantum 
kinetic equations. In order to obtain the classical equations one makes an expansion of T in powers 
of h 


— -^( 0 ) + ^-^( 1 ) + ^^-^( 2 ) + 


(11.32) 


and similarly of V,V^,A^,S^’^,tt and a. Expansions of the form ( 11.32 ) are inserted into Eqs. 
( 11.21 ) - ( 11.311) and the expressions at the appropriate powers of h are compared. The detailed 
description of this procedure will be the subject of the next Sections. 
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11.3 Constraint Equations in the Leading Order of ?i 


Substituting expansions of the type ( [L1.32| ) into Eqs. (11.21) - (11.25), we find to leading (zeroth) 
order of h 


_|_ 7r(o)P(o) 

= 0, 

(11.33) 

O'(0)^(0) + ^(O)-^(O) 

= 0, 

(11.34) 

P^^io) - ^(o)Vfo) 

= 0, 

(11.35) 


= 0, 

(11.36) 

- - 7r(o)5(;'; 

Correspondingly, the zeroth order of Eqs. (|11.27) - (|11.31|) has 

= 0. 

the form 

(11.37) 


= 0 , 

p^4°)+vr(°)4°) = 0, 

pi^pio) + = 0 , 


(0) 


= 0 . 


(11.38) 

(11.39) 

(11.40) 

(11.41) 


These are only four equations, since Eq. (11.27) is already of first order in h. From Eqs. (11.34) 
and (11.351) one finds 


^(0) = -T 




(0) 


(0) 


^(0) 


and 




^0) - P 


^(0) 


(11.42) 

(11.43) 


Then Eqs. ( 11.40|) and (11.41) are automatically fulfilled. Moreover, Eqs. (11.33) - (11.35) lead to 
the mass-shell constraint for the function .7^(0) (-T, p), namely 


[p2 - m2(X)].T-(0)(X,p) = 0, M\X) = CT^0)(X) + 4)(X). 


From Eq. (11.37), we find the following expression for the spin tensor 


‘^(0) - 


pMro, 


^(0) 

F “^(0) 


M2 


and the dual spin tensor 


qtMU _ _ 
‘^(0) - 




c^(o) 

F -^(0) 


(11.44) 


(11.45) 


(11.46) 


90 
































Substituting now Eq. { llAt ) into Eq. ( ll.St ), and using ( 11.3S ) one obtains a mass-shell constraint 
for A^q^{X,p) also, i.e., 


\p^-M\X)]A'l^^{X,p)=Q. 


(11.47) 


Einally, using (|11.45D , (|11.381) and (|11.47|) we find that (|1 1.351) is satisfied. In summary, Eqs. (|11.330 
- (11.41) lead to the expressions for P(o) in terms of and for and in terms 

of . Eurthermore, they lead to the two mass-shell constraints. 

The mean fields appearing in Eqs. (|11.33 ) - (|11.37| ), (|11.39| ) and ( |1 1.401) are required to be 
calculated in a self-consistent way from Eqs. (11.19|) and (|11.2CI|). Doing so, we find 


= -8g/ 


(11.48) 

(11.49) 


Due to the relation (11.34), one finds that these two equations are not independent and can be 
reduced to a single equation which determines the invariant mass M{X). To see this relation more 
clearly, let us define the distribution function F{X,p) via the expression 


F{X,p) = 


^{p){X,p) _ V(Q){X,p) 


^( 0 )(^) 


vr(o)W 


(11.50) 


and the angle 4>(X) via the relations 


7r(o) = M{X) sin $(X), cj(o) = M{X) cos 4>(X). 


(11.51) 


F{X,p) is a chirally invariant quantity, and knowledge of it determines the value of M{X). In 
addition, it allows us to calculate the vector current density However, Eqs. ( [11.48| ) and 

do not determine separately the fields 7r(o) and cr(o), i.e., the angle <h(f^). This fact is in 
agreement with the requirements of the chiral symmetry of the problem: under the chiral transfor¬ 
mations 4>(X) —> 4>'(X) = <I>(X) + x and, consequently, the absolute value of 4>(X) has no physical 
significance. 


(11.49) 


With this result, we conclude the discussion of the expressions obtained to zeroth order in h. 
The derived equations all display the chirally invariant form. Let us turn now to the discussion of 
the equations which follow from Eqs. (11.27) - (|11.31j), considered to first order in h. 


11.4 Kinetic Equation for the Quark Distribution Functions 


Employing Eqs. ( 11.42 ) and ( 11.43 ) in Eq. ( 11.27 ), one immediately finds 


p^d^F{X,p) + M{X)d^M{X)d^F{X,p) = 0, (11.52) 
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where we have used the definition ( |ll.50| ). This is again a chirally invariant equation. Due to 
the mass-shell condition ( 11.44 ), we can express F{X,p) as the sum of the quark and antiquark 
distribution functions /"'■(X, p) and f~{X,p) defined by 


p) + V- -p)-1])■ (11-53) 

Here Ep{X) = + M‘^{X), and p is a three-momentum. Substituting expression (|ll.53|) into 

Eq. (11.52D, and integrating over gives 


p>^dpf^iX,p) + MiX)dpMiX)diif^{X,p) = 0 . 


(11.54) 


Since the last procedure involves non-trivial cancellations let us discuss it in more detail. Eq. 
(|11.52 ) can be shortly written as LF{X,p) = 0, where L is a linear differential operator, i.e., 
L = p^dp + MdpMd^. After some algebra, we hnd that LEp{X) = {p^Md^M)/Ep{X) and 
L5{p^^Ep{X)) = MdQM5'{p^^Ep{X))[l^p'^/Ep{X)]^ where the prime denotes the derivative with 
respect to p^. Since L is a linear operator, one also gets L{[5{p^ =F Ep{X))/Ep[X)]f^[X,±.p)} = 
L [5(p0 ^ Ep{X))/Ep{X)\ /±(A, ±p) + [5(p0 ^ Ep{X))/Ep{X)]Lf^{X, ±p). Of course, if L acts 
on f^{X,±p) its energy component dp produces zero. Now using our previous results we can 
check that dp^L =F Ep{X))/Ep{X)'\ = 0, where A± is an interval containing either Ep{X) 
or —Ep{X). Consequently, integrating the whole formula (11.52) over A±, one gets the desired 
equation ( |11.54 ). We note that in the case of antiquarks still the replacement p —> — p should be 
made. 


Eq. ( |11.54 ) displays a typical form for fermion mean-field theories. In our case, however, a 
characteristic feature of this equation is the appearance of the chirally invariant mass M{X). We 
emphasize that the four-momentum p^ occurring in ( 11.54 ) is to be taken on the mass shell, i.e., in 
this case p^ = Ep{X). At the same time the distribution functions are functions of the space- 
time coordinate X and the three-momentum p, with the derivative with respect to the energy 
defined to be zero. 


Using Eqs. ( 11.4£ ), ( 11.50|) and ( 11.53 ) we hnd that the NJL gap equation takes the form 


a(o)(A) = a(o)(A) iG J 


d^p 


(27r)3 Ep{X) 


[l-/+(A,p)-/-(A,p)], 


(11.55) 


which determines the mass M{X) appearing in Ep{X) in terms of the distribution functions 
f^{X,p). One can easily notice that using Eq. ( 11.49|) instead of ( 11.48 ) we would obtain the 


same equation as ( 11.55| ), with the appearance of 7r(o) instead of (T(o). We note that Eq. ( 11.55 ) is 


the generalization of formulae ) and ( 5.17 ) to the case of non-equilibrium quark distribution 
functions. 
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Using now Eqs. (Ill.lOl) and (|ll.53|) , we find the following expression for the baryon current 


V^{X) = 2 


cfp 

(27r)3 Ep{X) 


[f+iX,p)-f-{X,p) + l], 


(11.56) 


where p^ is again on the mass shell, i.e., p^ = Ep{X). The constant appearing at the end of the 
square bracket can be neglected, since it does not contribute to V(X), and for V^{X) it introduces 
a constant (although infinite) charge. 


We note that Eqs. ( |11.54 ) and ( 11.55| ) form a closed system of equations: the distribution 
functions f^{X,p) determine the invariant mass M{X) through Eq. ( |11.55| ), whereas the space 
dependence of M{X) determines in ( 11.54 ) the time evolution of f^{X, p). In the case tt^q) = 0 this 
system of equations was first derived by Zhang and Wilets El , whose calculations were based on 
the closed-time-path formalism. Our approach represents an alternative derivation and generalizes 
their results to the case 7r(o) / 0. The numerical solutions of Eqs. ( 11.54 ) and ( 11.55 ), describing 
the expansion of quark matter possibly created in the ultra-relativistic heavy-ion collisions, were 
found in pl[| . Another class of solutions of this system of equations will be discussed by us in the 
next Chapter. 


11.5 Spin Evolution 


In the classical limit the spin dynamics is described by the behaviour of the function A^Q^{X,p) 
[^, ^]. In order to derive the kinetic equation for A|)j^(A,p), we use Eq. ( 11. 301) and substitute 


into it the expression for obtained from Eq. ( 11.23|) . Such a procedure gives, to the first order 
in h, the following equation 

= 0. (11.57) 

Substituting the expression for the spin tensor ( |ll.45| ) and the dual spin tensor ( 11.46| ) into Eq. 


(11.57) leads to the equation determining the space-time evolution of A|fg^(A,p), namely 


+ Md^Md^A’;^^ + 


d,M r 


( 0 ) 


M 




( 0 ) 


- = 0 . 


(11.58) 


We note that Eq. ( 11.58| ) is again a chirally invariant equation, since both the mass M{X) and 
the gradient of the angle <h(A) are chirally invariant quantities. Due to the condition ( 11.38| ), only 
three out of four equations in ( |11.58| ) are independent. In fact, one can check that multiplication 
of (11.58) by p^ gives zero. 


The mass M{X) appearing in ( 11.5S ) has to be calculated from the system of equations ( 11.54 ) 
and ( 11.55| ). Hence, it can be treated in (11.55) as an externally prescribed function. On the other 
hand, the gradient d^^{X) is not known, and until we specify how to calculate it, Eq. ( 11.58 ) 
cannot be used to determine the time evolution of A|)j^(A,p). Moreover, one has to check whether 
the solutions of (p.l.58|) satisfy the requirement of axial current conservation (3.3). The last two 
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points will be discussed in more detail in the next Section, where we examine the consistency of 
our all equations up to the first order in h. 

In analogy to the QED calculations @| we introduce the spin up and spin down phase-space 
densities 


FUX,p) = F{X,p)±S^[X,p)- 


M{X) ' 


(11.59) 


where S^{X,p) is defined by 


Sf.iX,p) = - 




{X,p) 


-A^^\x,p)Al^^{X,p) 


(11.60) 


Due to the condition (|ll.38|) , the four-vector A^Q-^{X^p) is space-like, A^Q-^{X^p)A\i\x,p) < 0, 
therefore S^{X,p)S^{X,p) = —1. The quantity S^{X,p) describes the mean spin orientation of 
the classical particles placed at the space-time point X and having the four-momentum p. On the 
other hand, the magnitude of A^qJX,p), through Eq. ( |11.59| ), determines how many spin up and 


spin down particles are present in the system. Using Eq. (11.581), we find 


d^M 


p''duS>^ + MduMd!^S>^ + = 0 , 


(11.61) 


and consequently 


p"5,F±,(X,p) + Md,Md;F^s{X,p) = 0. (11.62) 

Hence, similarly as in transport equations for QED [^], we find that the spin components decou¬ 
ple. Since both functions F{X,p) and A^Q^{X,p) are the mass-shell, we can write the analogous 
decomposition as in (|11.53|), namely 


F±siX,p) = 2 tt 


' 6{p<^ - EpjX)) 
2EpiX) 


f±siX,p) + 


6{p^ + Ep{X)) 

2 Ep{X) 


/-,(X,-p)-l 


Substituting Eq. ( 11.63|) into Eq. ( 11.62|) , and integrating over p^ gives 


(11.63) 


p^dpfi,iX, p) + MdpMd^fi^iX, p) = 0. (11.64) 


By a direct comparison 
Therefore the functions 
given in agreement with 


of (111.531) and (|11.63|) , we find f^{X,p) = i ^f^^{X,p) + f^^{X,p) 
f^{X,p) simply represent the spin-averaged densities. This has been 
the normalization of ( 11.55 ) and ( |ll.56| ). 
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11.6 Consistency of the Classical Transport Equations 


The kinetic equations for the quark densities (11.52) and for the spin ( 11.581) follow from Eqs. 
( 11.27 ) and ( |11.30| ) considered to the first order in h. We have yet to examine (i) whether equations 
( |11.28 ), ( |11.29| ) and (11.31) are consistent with these to the same order and (ii) whether the 
additional information for the angle <h(-T) appearing in ( 11.5£ ) can be obtained. We note that this 
kind of checking has been done in the case of QED in 


First, in order to illustrate that there is no contradiction between ( 11.27 ) and (|11.29 ) we sub¬ 


stitute the quantities .4^ and obtained from Eqs. (|11.24|) and (|11.25D , respectively into Eq. 


(11.29). In this way, we find the following equation 


^ - 


p 

^(0) 

1 


du Pu- 


Ti 


( 0 ) 


O-(O) 


-duipf. 


Ti 


( 0 ) 


^(0), 


TT 


(0) 
^(0) 


+ - 


^(0) 




dtjPiQ) 


= 0 . 


(11.65) 


Since we are interested in retaining terms to the first order in h we can use the relations ( 11.34 ), 
(|11.35D , ( |11.44| ), and the kinetic equation ( 11.52|) to find that (|11.65|) , and consequently ( 11.29 ), is 
fulfilled. 


In an analogous fashion, using the kinetic equation for spin ( |11.58 ), the constraints (11.38) and 
( 111. 47 ), and the definition ( 11.51 ), we can show that Eq. ( |ll.3lj ) is fulfilled again to the first order 
in h. To do this, we multiply both sides of ( 11.31|) by Sapf^u- Using the formula for from ( 11.23 ), 
we find that the multiplication of the left-hand-side of ( 11.31 ) gives the expression 


^(0) 


- fT(o) (daAf' - } • 


On the other hand, the multiplication of the RHS of (11.31) yields 


( 11 . 66 ) 


+ ^7^(0) 


07 c(0) 


(11.67) 


Using now the expressions for S^J and sla ), after a rather lengthy calculation in which Eqs. ( 11.38 ), 
( 11.47 ), ( 11.511) and ( |11.58|) are used, we find that the two expressions above are equal to each other. 
Consequently, Eq. (11.31|) is fulfilled up to the first order in h as well. 


In this way, we have shown that Eqs. (|11.27| ), (|11.29|) - ( 11.31 ) are satisfied in the first order. 
To prove this fact, we have used Eqs. (|11.21j ) - ( 11.25|) , treating them as if they were satisfied up 
to the first order. Is this procedure admissible? The answer to this question is yes. Only T and 
A^ are independent functions (with A^ restricted by the axial current conservation ( 11.26| )). Eqs. 
( 11.22 ), ( 11.23 ) and ( 11.25 ) are just the definitions of V, and 5^*^, and therefore are in our 
approach always satisfied. On the other hand Eq. (|11.21J) gives in the first order 


(p^ — — 2 (7(o)iT(i)-|-7r(o)vr(i) T{q-^A 


ho) 


( 11 . 68 ) 
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and 


{p^ - M'^)V(i) — 2 (T(o)cr(i) +7r(o)7r(i) V(p) + 5iyCr(o)^^o)> 


(11.69) 


whereas using Eq. (11.24) one finds 


(p2 - = 2 


f^(0)O'(i) + ^(0)7r(i) 




(11.70) 


One can notice that Eq. ( 11.70| ) as it stands is a chirally invariant expression. On the other hand 
Eqs. ( 11.681) and (|11.69 ) form a system of chirally invariant equations in the same way as equations 
(IlLlel) and (|TlT^) , or equations (|11.19|) an d (|11.20D . One can also check, that using formula 
( 11.22 ) as the definition of V and Eq. ( 11.68|) , we can derive Eq. (|ll.69[) . Therefore, only one of 
equations ( 11.68|) and ( 11.69 ) is really independent. We observe that all the expressions ( 11.68 ) 
- (|ll.7C| ) are just the generalized mass-shell constraints for .7^(1), 'P(i) and and they do not 
influence the relations in the zeroth order. 


The last equation that we are still required to check is Eq. 
V, calculated from (|11.22|) into (|11.28|) , gives 


(|11.28). Substituting the quantity 


M‘^( X) 

-^-Ld^^X)d^F{X,p) =pMW(X,p). 


(11.71) 


We see that Eq. (|11.71|) defines the parallel part of in the first order. Multiplying ( 11.701) by 
and using Eqs. (|11.71j) and (11.44) we obtain zero. Hence Eqs. (|11.71) and (11.701) are consistent. 


After deriving all possible equations up to the first order in h and checking their consistency, 
we are still faced with two problems. Eirstly, the system of our equations in the leading order is 
not closed, since we do not have an equation for <h(A). Secondly, it is still not clear whether our 
kinetic equation for spin satisfies the requirement of the axial current conservation. 


One a priori possible situation is that the condition of the axial current conservation elim¬ 
inates the freedom connected with the choice of 5^<1>(A). However, by studying some simple 
situations, one can convince oneself that this is not the case. For example, assuming that M = 
const and that initially the spin distribution is homogeneous = 0,x,p) = B^{p) we find 

~ di^pj Bk{p). In this situation, the axial current conservation at t = 0 requires 
that di^ f (Pj/p^) Bk{p) = 0. We can see that this equation is not sufficient to determine 

the gradient of and consequently the derivative remains not well defined. Consequently, 

in order to obtain the closed system of equations in the classical limit we have to make some 
assumption on the form of d^^{X). The simplest chirally invariant choice is 


d^^{X) = 0. (11.72) 

In the case (|11.72D , our kinetic equation ( |11.58| ) allows us to determine the time evolution of 
the function A^Q-^{X,p) from its knowledge at some initial time. One can notice, however, that this 
equation allows for solutions which do not satisfy the requirement of the axial current conservation. 
One of these solutions, for the case M = const, has the form A^g^{X,p) = s^{p)sp{p)X''5^^\p — p), 
where s^{p) is a vector satisfying the conditions s^(p)p^ = 0 and s^(p)s^(p) = — 1, and p is some 
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given value of the four-momentum, with = M^. On the other hand, there also exist solutions 
of Eq. ( 11.581) which satisfy the requirement of the axial current conservation. These are, e.g., the 
static homogeneous solutions (for the case M = const) or the solutions which are odd functions of 
momentum. 


In consequence, we observe that Eq. (11.58) neither guarantees nor contradicts Eq. (3.3). 
Therefore the axial current conservation should be used as an external condition which selects the 
physical solutions of (|11.5^). 


11.7 Explicit Breaking of Chiral Symmetry 


In the case m 7 ^ 0, the chiral invariance of the Lagrangian (^]^) is explicitly broken by the expression 
—In this situation, the additional term in the Dirac equation leads to a simple modification 
of our kinetic equations ( 11.141) - ( 11.18| ); we have to replace the mean field a by the sum a + m. 
Of course, in practice, this change affects only the real parts of our equations, i.e., the formulae 
( |11.21 ) - (11.25), since the imaginary parts (11.27) - ( 11.311 ) depend only on the gradients of a. On 
the other hand, we observe that Eqs. ( 11.48 ) and ( 11.49 ) remain unchanged because they are just 
definitions of the mean fields. 

The important fact in the non-symmetric case is that the relations between the functions 
.F(o),'P(o) and look different. In particular, instead of Eqs. ( 11.42 ) and (p.l.43|) , we now 
find 


^( 0 ) = -TT, 




( 0 ); 


( 0 ) 


DO) - P 


(T(0) + m ' ‘ cr(0) + m 

Using the first of the relations ( 11.73|) in Eqs. ( 11.48|) and ( 11.49| ), we obtain 

d^p 


(11.73) 


^7(o)(X) + 8gJ ^y(o)(X,p) = 0, 


(11.74) 


and 


vr 


( 0 )(^) 


<^(o)(-^) + + 8G f -^^JF(o)(X,p) 


= 0 . 


A straightforward consequence of these two equations is the condition 


(11.75) 


7^(0) (^) = 0 , 


(11.76) 


which implies also that V(^Q'j{X,p) = 0 and 4>(A) = 0. Therefore, we observe that if the system is 
not chirally invariant, equations (|11.48) and (|11.49f) allow for the determination of the angle <h(X). 


For m 7 ^ 0, the mass-shell constraints (11.44) and (11.47), as well as the kinetic equations 
(11.52) and (|11.58) preserve their form, only in this case we should use the definitions 
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F{X,p) = 


Ti 


( 0 ) 


(T(o) + m 


, M(X) = cr(o)+m 


(11.77) 


The gap equation ( 11.55f) has now the form 

cfip 


M(X) =m + 4GM / [1 - P(X, p) - /-(X,p)] , (11.78) 

Let us now discuss the behaviour of the axial current in the case when m ^ 0. In this situation, 
we no longer have an axial current conservation law, and using Eq. (11.22|), we can write 


^ “ (^(0) + - ’^(i)-^(0)| = 0. 


(11.79) 


The integration of (11.79) over momentum gives 


d,AUX) = -^n^,)iX). 




G 


(11.80) 


In contrast to the m = 0, case we observe that such an integration does not lead to any constraint 
for itself. Therefore, without any restrictions on dfj,A^Q^ we can treat Eq. ( 11.79|) as a defining 


equation for 'P{i). One can still check that Eq. (11.79) leads to the mass-shell constraint for 'P(i) 
which agrees with ( 11.69 ). The latter formula, for the case m 7 ^ 0, should be rewritten in the form 
[p^ — = di^MA'^Qy We note that Eq. ( 11.80|) describes the partial conservation of the axial 

current (PCAC) in our case. 
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Chapter 12 


Large Time-Scale Fluctuations of the 
Quark Condensate 


In this Chapter we deal with the transport equations, again reducing our interest only to quark 
dynamics in the mean (Vlasov) field. We shall consider the systems which are not uniform but 
close to equilibrium. This allows us to linearize the kinetic equations and to use a combination of 
analytical and numerical methods to look for solutions. 

Our main aim is to study the properties of a hadronic system in which the energy density is 
very large but which still exhibits a chirally broken phase. One could think of a fireball formed 
in high-energy heavy-ion collisions but at energies below the threshold for production of a chirally 
restored phase. The temperature range of interest for us lies between 100 and 150 MeV, and we 
assume that the critical temperature is 190 MeV. The baryon density is zero in our case, therefore, 
our considerations are appropriate for the description of the central rapidity region. We cannot 
study the systems at temperatures close to Tc, because one of the ingredients of our approach is 
that the mean field contains a large constant component in addition to the fluctuating part. The 
former vanishes (becomes negligible) at Tc and this fact makes our approach inappropriate for very 
high temperatures. 

Because the energy density of such a fireball is large, we hope that the description of its space- 
time evolution in terms of quarks is meaningful. As we shall see later, our equations will not lead to 
the separation of colour, although the NJL model does not include confinement. In principle, one 
has to take into account the dynamics of mesons as well. However, in practice it is not completely 
clear how this can be done in the situation out of equilibrium and in the fully satisfactory way, i.e., 
protecting all the conservation laws and symmetries. 

It is important to notice that already a quantum fluctuation of the condensate, i.e., the sigma 
meson, exists. However, due to its large mass the typical time scale over which the sigma meson 
wave function changes is a fraction of fermi: At < l/m-o- ~ 0.2 fm. In the approach presented here 
we study the large time-scale (classical) fluctuations of the condensate and neglect such oscillating 
short time-scale changes, although if the latter are taken into account the interesting phenomena 
can occur from couplings of the fluctuations on both scales. Nevertheless, the inclusion of these 
short time-scale fluctuations is equivalent to taking into account the mesonic degrees of freedom 
and, as it was discussed above, represents a separate, as yet unsolved problem. 
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The non-uniformity of our system will lead to fluctuations of the mean field around its average 
value that is determined by the thermal background. Since in the NJL model a change in the mean 
field is proportional to a change in the condensate, the latter will also fluctuate. As a consequence, 
through the use of our transport equations for the quark-antiquark plasma we are able to predict 
the space-time development of such fluctuations. 


12.1 Linearized Kinetic Equations 


The starting point for our investigations are Eqs. ( |11.54| ) and ( 11.78 ). Their straightforward 
generalization to the case of arbitrary number of flavours (Nf) and colours (iVc) has the form 


p^df^fiX, p) + M{X)d^M{X)d^f{X, p) = 0 


and 


1 = 


” +4iV=iV/G/'^i[l-/(X,p)] 


M{X) 


( 12 . 1 ) 

( 12 . 2 ) 


Here = (p^, p) is the four-momentum of a quark, = (t, r) is its space-time position, m is the 
current quark mass, and G is the coupling constant. The function f{X,p) is the sum of the quark 
and antiquark distribution functions, i.e., f{X,p) = /+(A, p) -|- f~{X,p). Each of the latter two 
functions satisfies equation ( |11.54|) so we can add and subtract such two equations. Assuming that 
/■’“(A,p) = /“(A,p) the equation for the difference of two distributions is trivially fulfilled and 


we are left only with Eq. (12.1). In the case Nf = 2 the constituent quark mass M(A) is related 
to the quark condensate by the relation 


MiX)=m-4G{q{X)qiX)), 


(12.3) 


which is a generalization of formula (5.28) for non-equilibrium situations. 

In our approach quarks are always on the mass shell, i.e., p^ = Ep{X) = \/pi^ + M^(A), but 
their mass is a function of space-time position. Therefore, the distribution function depends only on 
space-time variables and three-momentum p. Consequently, the derivative with respect to energy, 
p^, in (12.1) vanishes. However, we have included it in our expression (12.1) to show the relativistic 
invariance of the kinetic equation explicitly. The gap equation ( [12.2 ) is also relativistically invariant, 
but the integral appearing there is divergent |^. Eor the sake of simplicity and also in order to make 
connection to other calculations done in the framework of the NJL model, we shall regularize ( 12.2 ) 
using a 3-dimensional cutoff A. This procedure, of course, introduces an explicit breaking of the 
Lorentz invariance. In the following, whenever the integration is restricted to the region p < A we 
shall indicate this explicitly in the expressions. 

It is interesting to observe that the equality of the quark and the antiquark distribution functions 
(for each colour separately) has some nice consequences for us; we can think that on the average 
colour is compensated by anticolour at each point and, therefore, the discussed states of the plasma 
are colourless. In other words, although we describe the evolution of our system in terms of 

^The quark distribution function vanishes fast for p 
difference 1 — /(A, p), does not. 


oo, however the integrand in Eq. (12.2), containing the 
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the motion of quarks our equations will not lead to the separation of colour. (We note that 
our assumption concerning the equality of the quark and the antiquark distribution functions is 
reasonable only for the systems with net baryon number density equal to zero.) 

In the following we are going to consider the situation in which the distribution function can 
be written as a sum of the isotropic background distribution and a small perturbation, namely, 


f{t, r, p) = /o(p^) + Sf{t, r, p). (12.4) 

Correspondingly, the mean field can be written as 


M{t,r) = Mc + dM{t,r), (12.5) 

where Me is the constant part and 5M{t, r) is the temporally and spatially varying fluctuation. We 
shall assume that the background distribution is a thermal one 


/o(/) = 


,y/pMMl/T ^ ■ 


( 12 . 6 ) 


Here T is the temperature and the factor 2 comes from summing up the quark and the antiquark 
distribution functions as was discussed above (the chemical potential /r = 0 in our case). According 
to Eq. ( 12. 2|) the fluctuating part of the condensate can be written as 6{q{X)q{X)) = —5M{X)/AG, 
therefore, from now on we shall discuss simply the disturbances of the mean field. 

Substituting expressions ( 12.4 ) and ( 12. 5|) into Eqs. ( 12.1 ) and (|12.2| ) and keeping only the 
zeroth and first order terms in the deviations 6f{t,r,p) and 5M{t,r) gives us the following three 
equations: 


m 


rA 


(Pp 1 


l = — + ANcNfG, , 

Me " W {271)3 Ep 


1 - /o(p^) 


k i/(t, r, p) + V. A 5/(t, r, p) - ^ r) ■ ^ /„(p^) 


= 0 , 


(12.7) 

( 12 . 8 ) 


5M{t, r) 


^ .AT ^ d^P Me . 2\\ 

— + meNfGj ^ (l-/o(p)) 


= -4NeNfG /^^^<5/(t,r,p). 


(27r)3 Ep 


(12.9) 


Here Ep and v are the leading order expressions for the energy of a quark, Ep = ^Jp^ + M|, and 
its velocity v = p/Ep. 


Equation ( [12.7 ) is the zeroth order gap equation connecting the constant part of the mean field 
Me with the background distribution fo{p^) and has a form known from the NJL calculations at 
finite temperature (see, e.g., Eq. (5.16)). Solving (12.7) we find Me as a function of the temperature 
T. Assuming that m = 0 and fitting other parameters as in Ref. |]^|, i.e., Ne = 3, Nf = 2,G = 5.01 
GeV and A = 650 MeV, we find that at zero temperature Me = 313 MeV and the corresponding 
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value of the quark condensate is (-250 MeV)^. We note that this set of parameters gives the same 
temperature dependence of the mean field as the parameters of the one-flavour model used in 
Section 5.1. 



Figure 12.1: Temperature dependence of the background mean field Me- The values of the parameters 
were taken in such a way as to get the agreement with the temperature dependence shown in Fig. [5.3]. The 
temperature range of interest (together with the corresponding range of Me) is framed by the box. 


12.2 Conservation Laws 

Because we neglect the processes of qq annihilation the total number of quarks and antiquarks 
should be conserved. We can write this conservation law formally as 

^^N^^{X) = 0 , ( 12 . 10 ) 


where the number current is defined as 


N^^{X) = 2NcNf 


cfp 

(27r)3 p^ 




( 12 . 11 ) 


In fact, starting from ( p.2.11| ) one can derive ( 12.10| ) using the kinetic equation (12.1) and integrating 
by parts. Similarly, we can show that the entropy of the system is conserved, namely. 


d^s^^ix) = 0 , 


( 12 . 12 ) 
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where 


5^(X) = -2NcNf J [f{X, p) In f{X, p) + [1 - f{X, p)] In [1 - f{X, p)]] 

is the entropy flux. The latter is conserved because we consider a collisionless plasma. 
Let us now define the energy-momentum tensor of quarks in the standard way by 


(12.13) 


Using again the kinetic equation ( 12.11) and integrating by parts one obtains 


d^p 1 


d,T^‘:^,,^,{X)=2N,NjMiX)d^MiX) J j^^f{X,p). 


(12.14) 


(12.15) 


The result ( 12.15 ) shows that the total energy and the total momentum of quarks alone are not 
conserved; this is due to the fact that ( 12.14D includes only the kinetic energy of particles. The 
question arises: What tensor, describing the potential energy, should be added to ( 12.141) in order 
to obtain a conserved quantity? In the thermodynamic approach of |15, 16| it was found that the 
expression 


^vac(^) — 


[M{X)-mf d^p 


4G 


- 2N,Nf 


(2vr) 


/ {p^ = Ep{X)) 


(12.16) 


represents the energy density of the vacuum (with the corresponding vacuum pressure given by 
—nvac(-^))- Because the vacuum is isotropic we write its energy-momentum tensor as 


T^cix) = g^^nyacix), 

where is the metric tensor. Using now the gap equation (|12.2) we can show that 


d^TZciX) = -2N,NfM{X)d^M{X) 


rA 


d^p 1 
(27r)^ p^ 


f{X,p). 


(12.17) 


(12.18) 


We see that for A —> oo the summed energy of the quarks and of the vacuum is conserved. This is 
an interesting result showing that one has to include the energy of the zero modes of quarks (the 
second term on the RHS of Eq. ( p.2.16| )) to have the conservation law. Since, A is finite in our 
approach, it means that the energy is not exactly conserved. However, as long as we look at the 
leading terms of the quantities in our expansion dehned by Eqs. (12.4) and ( |12.5 ) we hnd that this 
violation is very small and, therefore, can be disregarded. 

One can check that it is possible to obtain a strict conservation of energy if the gap equation 
( 12.2 ) is regularized in a slightly different way: The integral on the RHS of ( 12. 2| ) should be 
written as a sum of two terms. The first one, independent of f{X,p), has to be regularized with 
the cutoff, whereas the second term containing f{X, p) is converging and should not be additionally 
regularized. Such a scheme yields the exact conservation laws but at very high temperatures it can 
lead to negative and, therefore, non-acceptable solutions for the mean held — this is the situation 
discussed in Section 5.1 below Eq. (^.14 ). 
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12.3 Energy-Momentum Representations 


We now come back to the discussion of Eqs. ( 12.81) and (|12.g ). Since these equations are linear, it 
is useful to expand the functions 5M{t,r) and 5f{t,r,p) as Fourier integrals 


5fit,r,p) = 


p“s/(t.k,p)e-k.. 

(12.19) 

d^k 

(12.20) 


and substitute them into ( 12.8 ) and ( 12.9 ). In this way we find two coupled equations for each 
Fourier component separately 


d 


k, p) - ^ SM (t, k) k ■ ^ = 0 


Er, 


6M{t,k) 




5p 

^ d^p Me 


( 12 . 21 ) 


fl'^T) A// 

^^5/(t,k,p). (12.22) 


In the situation under consideration the system is perturbed at t = 0 and the perturbations 
propagate forward in time, i.e., (5M(t, k) and 5f{t,k,p) are only different from zero for t > 0. 
Consequently, to solve ( 12.21 ) and (|12.22D , we use a one-sided Fourier transform in time |^. For the 
fluctuating part of the mean field it is defined as 

roo 

5M{uj,'k) = / dt. (12.23) 

JO 

Assuming \6M{t,k)\ < M exp{at), where M and a are some positive numbers, we find that Fq. 


(12.23) is well defined for frequencies uj having an imaginary part ujj > a. Thus, the inverse 
transformation is given by 


/ OQ-\-ia , 

-oo+ia 27r 


(12.24) 


where the integration contour is taken to be a straight line parallel to the real axis. We note that 
the assumption discussed above means that all singularities of 6M{uj, k) lie below the contour. The 
one-sided Fourier transform of 5f{t,k,p) is defined in the analogous way. 


Multiplying (|12.21) by exp(itjt) and integrating over time from zero to infinity gives 


<^/(w,k,p) = T 


i(k ■ V — Lv) 


g{k, p) -|- 6M{io, k) k • 


Er, 


5p 


(12.25) 


^In what follows we adopt the Landau method ©El for studying the relaxation of the initial perturbations in 
electron-ion plasmas. 
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where ^(k, p) represents the initial condition for the distribution function 6f{t,r,p), namely, 


r 

= 0, r, p) = g{r, p) = J p)e*^ 


(12.26) 


Making the one-sided Fourier transformation of Eq. (12.22) and substituting into its RHS Eq, 


(12.25) we find the following expressions 


(5M(a;,k) 7 ( 0 ;, k) =G(a;,k), 


where 


TD 

7 (a;,k) = —+ 4iV,iV;G j 


(fp 


dfo 


1 


and 


G(u;, k) = AiNcNf G J 


^ cPp Me g(k,p) 


(12.27) 


(12.28) 


(12.29) 


{2ttY Ep k • V — cj 

Equation (|12.27|) is very much similar to that found in the case of investigations of the oscilla¬ 


tions of the electron-ion plasma |75]. Relying on this analogy, we shall call 7(0;, k) the permittivity 


of the quark-antiquark plasma. Substituting Eq. (12.27) into (12.24) we obtain 


J-oo+ia^TT 7(^^>k) 


(12.30) 


The analysis of this equation will be the subject of the next two Sections. 


12.4 Dispersion Relation and Analytic Properties of the Functions 
7 and G 


We already know that the integration in ( 12.301) must be performed above all singularities of the 
function G(a;, k)/ 7 (a;, k). Eq. ( 12.28 ) defines 7 ( 0 ;, k) as the function of the complex frequency, 
uj = LOR + iuJi, which is analytic everywhere except for the interval on the real axis where < 

< v^^k. Here = A/a/A^ -|- is the limiting maximal velocity of quarks, which because 
of the presence of the cutoff is slightly smaller than 1. The structure of the singularities of G(a;, k) 
is analogous to that of 7 ( 0 ;, k); they all he on the real axis in the interval defined above. The 
additional singularities which can appear in the integrand of ( |12.30| ) come from the zeros of the 
denominator, i.e., they are solutions to the equation 


7 ( 0 ;, k) = 0 . 


(12.31) 


Equation ( 12.31 ) is nothing else than the dispersion relation defining modes in the plasma, which 
have momentum k and energy uj{k) such that 'y[uj{k), k] = 0. 
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Let us now try to find a solution to Eq. ( 12.31 ). We first perform the angle integration in 


(12.28) and find the following representation of the function 7 ( 0 ;, k): 


, ,, m 2N^NfGM^ , p 

7(o;,fe) = —+ dp^ 


Mr 


TT^ 


/O 


En 




(12.32) 


where 


I / s 1 , 2 : + 1 

^(z) = -zhi - 1 . 

^ ^ 2 z-1 


(12.33) 


One can notice that $( 2 :) has a cut for real values of z contained in the interval — 1 < Zi? < 1. This 
is in agreement with an earlier remark that 7 ( 0 ;, h) has a cut for < (jOr < v'^^^k. 

For the arguments slightly shifted from the real axis we find 


$(z = Zij ± ie) = -ZR 


In |j =F *709(1 - |zr|) 

\zr ~ f | 


- 1 , 


(12.34) 


where 6{x) is the step function. One of the consequences of Eqs. ( 12.321) and ( 12.34 ) is that the 
dispersion relation (|12.31|) can have a solution for real frequencies {u! = ujr) only if |wij| > kv^ 


Otherwise, 7 ( 0 ;, k) has a non-vanishing imaginary part. A closer examination of the function ( 12.33 ) 
shows, however, that for |u;i?| > kv'^^^ the real part of ^{EpUjR/pk) is always positive and so is 
the real part of 'y{LV,k). To see the last point one can notice that the derivative of the thermal 
background distribution is always negative. Consequently, we conclude that the dispersion relation 


(12.31) has no solution for real frequencies. 

For purely imaginary arguments 4>(z) has the form 


4>(z = izi) = \zi\ — arctanlz/l^ — 1 , 


(12.35) 


hence it is a negative real function. In this case ( |12.32| ) is a difference of two real positive functions 
and one can numerically check if there exists a solution to Eq. ( 12.31| ). In the temperature range 
of interest, i.e., for 100 MeV <T< 150 MeV, we have found that 7 ( 0 ;,/c) is always larger than 
zero. This fact again means that our dispersion relation has no solution. 

To complete our investigation we have studied numerically the imaginary part of 7 ( 0 ;, k) for 
ujR,u}j > 0 again in the temperature range 100 MeV < T < 150 MeV. We have checked that Im 
7 ( 0 ;,/c) is always negative in the discussed region and vanishes at infinity, i.e., for |ti;| ^ 00 . 
Employing the symmetries Im 7 (a;ij — iuji,k) = —Im 7 (cui? -|- iuji,k) and lm.'^{—ujR + icoj^k) = 
—Im 7 (a;R -|- iu}j,k) we find that Im 7 (a;,A:) is different from zero also in other quadrants. This 
result, together with the two previous ones, indicates that there is no solution to Eq. ( 12.31[) for 
any finite uj. 


Because the only singularities of the integrand in ( 12.30|) are those coming from the cut discussed 
earlier (placed on the real axis) our parameter a can be taken to be an arbitrarily small positive 
number. For t < 0 we can close the contour of integration in the upper half-plane. Because the 
integrand in (|12.3C|) has no singularities in this region, we find in this case that 6M{t, k) = 0, 
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what is in agreement with the requirements of causality. For t > 0 we can close the contour of 
the integration in the lower half-plane and, later, we contract it to go just around the cut. As the 
result of this procedure we obtain the expression 


6M{t,k) 







'G(a; + ie, k) 
.7(0; -t- ie, k) 


G(a; — ie, k)' 
7(0; - ze,k)_ ■ 


(12.36) 


For further use, it will be convenient to represent 7 ( 0 ;/? ± ie, k) explicitly as the sum if its real 
and imaginary parts, i.e.. 


7 ( 0 ;^ ± ie, k) = a ± > 


(12.37) 


where 


m 2N^NfGMl 
«(^) = TT + ^ 


M, 


TT^ 


dpp 

En 


p L P 




^ f ^ dfo, \n+p/Ep 


[ dp^\n 

J op 


2 J ^ dp — p/Ep 


(12.38) 


and 


^ MiPMla 


TT 


/o(A) - /o 






(12.39) 


The reader may be surprised by our result, since we have a contribution from the cut but, on 
the other hand, we do not hnd a solution to the dispersion relation ( 12.31 ). In the standard case 
[73, 75] the situation is quite opposite: the permittivity is an entire function of the frequency and 
one does hnd solutions to (|12.31D . The point is that for a non-relativistic Maxwellian plasma one 
can analytically continue 7 (w, k) from the upper half-plane into the lower half-plane. In this way, 
one obtains a function 'y'{uj,k) such that (i) for w/ > 0 it coincides with 'y{uj,k), (ii) for w/ = 0 
it has no cut, (iii) for a;/ < 0 it is different from 'y{uj,k), (iv) it is an entire function of to, and 
hnally (v) there are solutions to the dispersion relation 7 ^( 0 ;, k) = 0. Because these solutions are for 
complex energies in the lower half-plane it means that the modes are damped — this is the famous 
Landau damping © caused by a transfer of the energy from the wave to the particles of 

the medium. 


In our case, although the mathematics is slightly different we can expect similar physics. Our 
contribution from the cut is a superposition of waves whieh ean exhibit destructive coherence phe¬ 
nomena, i.e., they ean he damped as in the standard ease. The main mathematical difference 
between our approach here and the typical non-relativistic one is that we do not perform the ana¬ 
lytic continuation of 7 (w, k) from the upper half-plane into the lower one because, in our case, this 
does not appear to be a straightforward procedure. 


12.5 Results 

As the initial conditions for our equations, we shall consider two extreme cases which lead to 
spectacularly quite different behaviour of the plasma. The first case corresponds to the perturbation 
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which is at t = 0 independent of the momentum p, whereas the second case describes the situation 
when the initial perturbation is strongly peaked at some value of p. The second case can be 
regarded as a more interesting one, because any initial perturbation (including the first one) can 
be treated as a superposition of such peaked elementary perturbations. Moreover, as we shall see 
later, the solutions are also superpositions of elementary solutions corresponding to such initial 
conditions. Let us, however, start our considerations with the discussion of the first case. 


12.5.1 Initial Perturbation Independent of Momentum 

Writing the initial perturbation of the quark distribution function in the form 


9 {Kp) = 51 (k), 

and using the above expression in (|12.29|) one can hnd that 


G(u;ij±z€,k) = 




where 


and 


A{n) = 


NcNfG 


27r 


, 2 

— - -A 

1 - 


A 


B(fi) = I (j^pp In 


vr^ 


_ |5^|) 

\n + p/Ep\ 

\i}-p/Ep\ ■ 


(12.40) 


(12.41) 


(12.42) 


(12.43) 


Equations ( 12.361 ), ( 12.37 ) - ( 12.3£ ), and ( 12.41 ) - ( 12.43 ) allow us to write the final expression 
for the fluctuating part of the mean field as 


5M{t,k)=giik)M,Ei{kt), 


where 


'U 

Fiikt) = J 


dn aA + f3B 

TT a2 + /32 


cos{Q,kt). 


(12.44) 

(12.45) 


One can obtain the space-time dependent function 5M{t, r) by calculation of the Fourier trans¬ 
form ( |12.20 ). The time evolution of the fluctuations is completely determined by the function 
Ei{kt). 

In Fig. [12.2] we show the function Fi{kt) calculated numerically for two different values of the 


temperature: T = 100 MeV and T = 150 MeV, respectively. The parameters were fitted as in j^j 


see also the discussion following Eqs. (12.7) - (12.9). For small time arguments Fi{kt) is negative 
as it should be — if we add particles to the system, gi (k) > 0, the mean field decreases. For t = 0 


the change of the mean held can be independently calculated either from Eq. ( 12.44 ) or directly 
from Eq. (12.9). The condition that two methods give the same result was used by us as a check 
of our numerical code. 
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Figure 12.2: Function Fi{kt) calculated for two values of the temperature: T = 100 MeV and T = 150 
MeV. The values of other parameters are: m = 0, Nc = 3, Nf = 2, G = 5.01 GeV“^, and A = 650 MeV. 


For both values of the temperature we observe that Fi{kt) oscillates in time but with the 
amplitude strongly decreasing. In our opinion, we observe here a phenomenon having features in 
common with the Landau damping in ordinary plasmas, i.e., the energy of the wave is transferred to 
the medium without entropy production, see Eq. ( 12.12| ). There exist, however, differences between 
the damping observed here and the standard Landau damping. First of all, one can check that the 
suppression of the oscillations is much weaker than exponential. Moreover, the damping seems to 
be independent of the temperature. In fact, our calculations indicate that it is a relativistic effect 
connected with the existence of the finite maximal velocity, i.e., the velocity of light. The latter 
determines the size of the integration interval in ( 12.45| ) — in practice we have always ~ 1. 
For the considered values of the temperatures one can check that the function {aA + PB)/{a‘^ + /3^) 
is very regular and, therefore, we expect that for kt = 27r the integral becomes negligible. This 
result is due to the fact that cos(27rn) takes on equally positive and negative values in the interval 
between -1 and 1. By inspection of Fig. [12.2], one can convince oneself that the initial perturbation 
of the mean field becomes strongly suppressed for kt > 2tt. 


It is perhaps also interesting to observe that our suppression of oscillations has some features in 
common with the Friedel oscillations analyzed in Chapter 8. The latter are caused by the existence 
of the sharp Fermi surface, in our case there exists a sharp maximal velocity Ri 1. 
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12.5.2 Initial Perturbation Strongly Peaked in Momentum 


Let us now turn to the discussion of the opposite case, where the initial perturbation is strongly 
peaked in momentum. Employing (12.29) and (12.36) and exchanging order of integration we 
obtain the formula 


5M{t,k)=4N,NfG / 

J (zvrj hjp 


where we used the definitions 


kv^ 


(5M(t,k;p)= J 



■G(cu + 7e,k; p) 

G(a; — ie, k; p)' 

J 27r 

. 7 ( 0 ; + ze,k) 

'y{uj-ie,k) . 


_^„max 


and 


G(a; + ie,k;p) = 5 r(k,p) 


iMr 


k • V — cu =F ie 


(12.46) 


(12.47) 


(12.48) 



Figure 12.3: The f unction F2c(0.5, kt) plotted for T = 100 MeV. The dashed line shows the complete result 
obtained from Eq. ( 12.50 ), whereas the solid line is the a pprox imation for F2c{0-5,kt) obtained by keeping 
only the pole contribution (the first term on the RHS of (|12.50)). Other parameters as in Fig. [12.2]. 


Equations (12.46) - (12.48) have a clear physical meaning, they say that the fluctuation of the 
mean field (12.46) can be treated as a superposition of the modes ( p.2.47| ) which correspond to 
elementary initial conditions ( 12.4^ ). Such elementary initial conditions are realized by adding to 
the system, at t = 0 , particles with well-defined momentum p. 

Using Eqs. ( 12.47| ), ( 12.37 ) - ( 12.39 ), and ( 12.48 ) we find that 


no 































SM{t,k;p) = g{k,p)Mc 


(12.49) 



vk 

~T~ 



iF2s 


vk 

~~k~ 



where 


F 2 cix,y) = -cos(xy) 


a(x) 


a^(x) + /3^(x) 


^ f dn , 
+ P / — cos{Qy) 

J TT 


Pin) 


a2(0) 


(12.50) 


and 


F 2 six,y) = -sin(xy) 


a{x) 


a^{x) + P‘^{x) 


+ 7 ? f sin(Qi/) 
J vr 


Pin) 


a'^in) + p'^in) X - n' 


(12.51) 



Figure 12.4: The same as Fig. [12.3] but for T = 150 MeV. 


In Figs. [12.3] and [12.4] we show the function F 2 c(v • k/k, kt) for two different values of the 
temperature, i.e., for T = 100 MeV and T = 150 MeV. All the parameters are fitted in the same 
way as in the previous cases. We fixed the projection of the velocity v on the wave vector k to 
be V • k/k = 0.5. The dashed curve represents in both cases the complete result obtained from 
( 12.50D, wh ereas the solid line shows only the pole contribution (only the first term on the RHS of 
Eg. (|12.^ ). Similarly, in Figs. [12.5] and [12.6] we show the function F 2 s(v • k/k,kt). 

One can see in Figs. [12.3] - [12.6] that the functions F 2 ci'v-k/k, kt) and F 2 si'v-k/k, kt) oscillate 
in time without any damping. It means that our elementary fluctuations (|12.49| ) will be undamped 
as well. This phenomenon is caused by the singularity of the initial quark distribution function. 


We want to point out that one finds similar waves in the theory of ordinary plasmas in the form 
of so-called Van Kampen modes [^|. Landau’s calculation emphasized the case where the 
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Figure 12.5: The function F2s(0.5, kt) plotted for T = 100 MeV. The dashed line shows the complete result 
obtained from Eq. (12.51), whereas the solid line is the approximation for F2s{0-5,kt) obtained by keeping 
only the pole contribution (the first term on the RHS of (12.51)). 


initial perturbations were entire functions of the momentum. The Van Kampen modes are excited 
by the perturbation being a 6 function in momentum. Of course, such a perturbation is not an 
entire function in the complex to space. It gives rise to singularities on the real axis and leads to 
undamped oscillations. 


In Figs. [12.3] - [12.6] one can also see that the functions F2c(v • k/A:, kt) and F2s(v • k//c, kt) 
can be quite well approximated just by the first terms on the RHS of Eqs. ( 12.501) and (|12.5lD . 
Using this fact we can find the following approximate form for the fluctuating part of the mean 
field: 


6M{t, r) Ri -4NcNfG J 

g{k-p) +g{-'k-p 


^ d^p Me f d?k a(v-k/A;) 

(27r)3 Ep J (27r)3 a'^(y ■ k/k) + /3^(v • k/k) 


^ w 1 ^ , ^(k • p) - ^(-k • p) . , 1 , 

cos(v • k t — k • r) H-^-sm(v • k t — k • r) 

2z 


(12.52) 


Equation ( 12.52| ) shows that an arbitrary fluctuation can be written as a sum of undamped waves. 
Nevertheless, such a superposition itself is very likely to be damped because of the destructive 
coherence phenomena among the elementary excitations; this is the reason why the fluctuations 
are damped in our first case. 
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Figure 12.6; The same as Fig. [12.5] but for T = 150 MeV. 
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Chapter 13 


Critical Scattering at the Chiral 
Phase Transition 


The calculations based on the NJL model indicate the presence of singularities in the elastic and 
hadronization cross sections This is a manifestation of criticality at a second order phase 

transition which favors the production of low-energy pions. In this Chapter we study the time 
evolution of the transition from partons to mesons and investigate to which degree the \ow-pT 
enhancement observed in ultra-relativistic heavy-ion collisions may reflect criticality. 


13.1 Critical Opalescence 

Critical scattering (called also a critical opalescence) is a rather general phenomenon observed in 
the vicinity of a second order phase transition |^: The cross section d?a/d(\ for the scattering of 
light, X-rays or neutrons on a medium is given by the expression 

d^r J ^(n(r )n(r')) 7 ’, (13. 

where q is the momentum transfer and (n(r )n(r' ))'r is the density-density correlation function for 
the particles in the medium at the temperature T. In the neighborhood of a phase transition a 
long-range ordering develops and the cross section for small values of q has the form 



where 


d^a _ 1 

(13.2) 

dq fe2(r)-|-q2’ 

e{T) (X11 - ^ r 

(13.3) 


is the inverse correlation length with 7 being the critical exponent. Criticality manifests itself by 
the fact that the system becomes opaque to an external observer. 
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An experiment, where an external test particle is scattered on the system in its critical state, 
cannot be performed for a quark-gluon plasma which only lives for 10“^^ s. Nevertheless, also 
particles which are created inside the plasma can experience critical scattering. One has to find, 
however, a suitable observable, which makes the internal critical scattering visible to an outside 
observer. In the following we shall argue that critical scattering leads to a faster production of low 
momentum pions and thus to an intermediate non-thermal distribution. The pions which escape 
to the detector, before the thermal equilibration of the hadron gas sets in, can be observed as an 
enhancement at low momenta of the pion distribution. 


A low-pT enhancement in pion and kaon spectra is indeed observed in ultra-relativistic heavy- 
ion collisions. At small values of the transverse mass, ttit = the observed distributions 

of pions and kaons, dN/duiT, considerably exceed the thermal distribution [^, ^]. Various 

effects have been proposed to explain this phenomenon; e.g., existence of the transverse flow |^, 8C] 
creation of small plasma droplets 
of chemical equilibrium K 


, decays of resonances |82], formation of the pion system out 
or the medium modification of the pion dispersion relation f^. 


To which degree the critical phenomena contribute to the observed phenomenon, is the subject of 
this Chapter. 


13.2 Kinetic Description of the Hadronization of the Quark Plasma 

The phase transition from a quark-gluon plasma to a hadron gas has two aspects: chiral restoration 
and deconfinement. The NJL model is a useful framework to study the chiral aspect. For a 
vanishing mass of the current quarks, a phase transition of the second order is predicted at a 
critical temperature Tc above which the constituent quarks have zero mass while below Tc the 
pions are the Goldstone particles with = 0. 

In explicit calculations |^, using the NJL model for a description of the quark-meson plasma 
(no gluons), it has been shown that singularities occur in the cross sections at small center-of-mass 
energies y/s of the colliding particles, when the phase transition is approached. In particular, the 
integrated elastic quark-antiquark cross-section aqg^qq{s, T) diverges like for T ^ Tc. A singu¬ 
larity is also found for the hadronization cross-section qq —> tttt. In both cases the singularity arises 
because the quark condensate {qq)T, which is the order parameter of the chiral phase transition, 
goes to zero at the critical temperature. 

Although these results have been obtained for a particular model they may be of more general 
relevance, provided the parton-hadron phase transition is second order. We use the results for the 
energy dependence of the cross sections in order to model the time evolution from quark to hadron 
matter: The cross sections are incorporated into the kinetic equations proposed for the description 
of the transition of a spatially homogeneous quark plasma (described by the single particle quark 
distribution function fq{p,t), p = |p |) into a pion gas (with a corresponding distribution function 

In order to keep the model as transparent as possible we work in the relaxation time approx¬ 
imation. However, in contrast to the usual approaches, the relaxation times r(p, t) depend on 
momenta via the energy dependence of the cross sections. The singularities in the cross sections 
due to criticality of the medium translate themselves into singularities of the inverse relaxation 
times l/r(p, t) in the momentum p. 
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The above ideas are translated into the following set of kinetic equations 


dfq{p,t) ^ fqip,t) - fq,th{p,t) fq{p,t) ^ U{p,t) 

dt nh{p,t) Thad{p,t) Tdec{p,t)' 


(13.4) 


df-K{p,t) 

dt 


fw{P,t) - f^,thiP,t) fq{p,t) _ fn{p,t) _ U{p,t) 

Tthip,t) Thad{P,t) Tdec{P,t) '^em (p) ’ 


(13.5) 


dfem{P,t) 

dt 


Tem{p) ' 


(13.6) 


Eq. ( 13.4 ) determines the time evolution of the quark distribution function fq{p, t) (it describes 
both qnarks and antiquarks, i.e., fq{p,t) = /quarks^?*’■^antiquarks^^*’term on 
the RHS of Eq. (13.4) is the collision term written in the relaxation time approximation; it 
is responsible for the thermalization of quarks since the distribution fnnction fq{p,t) is always 
attracted to the thermal one fq^th{p-,t)- The second (third) term describes the loss (gain) of quarks 
due to the hadronization (deconfinement) process where we limit ourselves to the reaction qq —> tttt 
{qq ^ tttt). 


In Eq. ( 13. 5| ), for the time evolntion of the distribution functions of pions, the first three terms 
on the RHS describe the thermalization of pions, appearance of pions due to hadronization of 
quarks, and the two-pion reaction into qnark-antiquark pairs, respectively. The last (fourth) term 
accounts for the emission of pions from the plasma into the detector. For t ^ oo the distribution 
function femip,t) describes the observed spectrum of pions. 

Eqs. ( 13. 4| ) - (|13.(: ) can be solved for any initial conditions (i.e., assuming some particular form 
of the distribution fnnctions at the initial time t = 0) provided the thermal distribution functions 
fq,th{p-,t) and f-K,th{p,'t) are known at all times. Since the first terms on the RHS of Eq. ( 13.4 ) 
and ( |13.5D represent the collision terms (in the relaxation time approximation) they must obey the 
symmetry leading to the energy conservation. This gives the following constraint for fg th{p^ i) and 
fn, th{p,t) at each time t 


d^p 

(27r)3 


7 "^ ^ fi{p, t) - fidhip, t) 

p^ + ^- = 0 


rth{p,t) 


(13.7) 


(here i = q or ir). Eq. (13.7) determines the temperatnre Ti(t) appearing in the thermal distribn- 
tions. For simplicity we assnme the low-density high-temperature form of these functions, namely 
a Boltzmann distribution 


fidh{p,t) = s-iexp 


_|_ ^2 
T^it) 


(13.8) 
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where gi are the degeneracy factors: gq = 24 (quarks and antiquarks having two different spin 
projections, 3 colours and 2 flavours) and g-^ = 3 (three different values of the isospin). One 
can notice that the distributions ( 13. 8| ) correspond to the case when the chemical potential fii is 
zero, consequently we allow in our approach for number changing processes like: qq 2(gg) or 
2(7r7r). 


TTTT 


13.3 Relaxation Times and Cross Sections 


The analysis of the exact collision term in the Boltzmann kinetic equation (analogous to that in 
[^) leads to an expression for the average time for hadronization of two quarks into two pions 


1 


1 


d^pi 


ThadiPj't) \Jpi + 


fqiPl-i t)F'qq{s) ^qq —^TTTT (^); 


(13.9) 


where cjqq^ 7 r 7 r('S) denotes the total cross section for this process. The relativistic flux factor of 
incoming quarks is Fqq{s) = ^ y^s(s — 4m^), with ^/s being the center-of-mass energy of the quarks 
with momenta p and pi, respectively. The expression giving the deconfinement relaxation time, 
TdeciPit), has the form similar to Eq. ( 13. 9| ). In the analogous way one also defines the relaxation 
time for thermalization Tth{p,t) which is mainly due to elastic processes {qq —> qq,qTT —> qir) with 
the respective cross sections. 

In the chiral limit and at the temperature of the phase transition, only for which the singularities 
in the cross sections are observed, the masses of quarks and pions have to be set to zero. We shall 
also assume that all cross sections: (Tqq^^ 7 r(s), cjgg^gq(s), cj 7 r 7 r^ 7 r 7 r(s) and aq-jr^gTris) are given by the 
generic expression 


a{s) = fJo 



(13.10) 


where (Tq represents a constant contribution to the cross sections, and the appearance of the singular 
term sq/s accounts for the critical phenomena. The hadronization and deconfinement cross sections, 
(Tgg^^ 7 r(s) and (T^ 7 r^qq(s), are related to each other by the principle of detailed balance 


(13.11) 


This relation guarantees that in the absence of emission (rgm oo), the kinetic equations lead to 
the chemical equilibrium n^/n^ = gqjg-,r- 


The possibility of the emission of preequilibrium pions, i.e., the emission of pions at each stage 
of the hadronization process is crucial for our description. In the model defined by the kinetic 
equations ( 13.4 ) - ( 13.6 ) this emission is only treated in a very global and rather crude way. The 
emission time Tem is a free parameter of the theory, which allows us to decouple a certain fraction 
of pions from the interacting, non-equilibrium system. We have no good model for the momentum 
dependence of Tem- Therefore we test two assumptions: 
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7 ^ 


(13.12) 


and 


1 

Tth{p = Tc,t = 0) 


Ip TZ 

Tc Tthip = Tc,t = 0)' 


(13.13) 


In the first case the emission time is independent of momentum, whereas in the second case the 
inverse emission time depends linearly on p (the emission of low energetic pions is suppressed). 
The magnitude of TZ indicates how much faster the emission is, in comparison to the rate of the 
thermalization processes inside the system. 


13.4 Results of Numerical Calculations 


In order to reduce the number of free parameters in our calculation we have set all the cross sections: 
o’qq^Tnris), CTqq^qq^s), and ('?) to be equal (within a factor of 3 such a result is 

supported by calculations within the NJL model) and use the form (13.10) which is appropriate 
for T = Tf.. Then, we are left with three parameters: do, sq and TZ. One can notice that (Tq sets the 
overall time-scale and is unimportant if we are interested in the final {t — > oo) pion distributions. 

We start solving our kinetic equations by assuming that the initial quark distribution is a 


thermal one and that there are no pions in the system. After integrating Eqs. ( 13.4 ) - ( 13.6 ) till 
the time when all quarks are hadronized and all pions are emitted, one obtains the distribution 
function fem{p,t oo) of the observed pions for a set of parameters sq and TZ. 


In Fig. [13.1] we show our results for the initial temperature of quarks Tq{t = 0) = Tc = 140 
MeV, and for sq = 0.5 GeV^ and TZ = 10. The two dashed lines represent the distribution functions 
of the emitted (observed) pions. The upper one (a) corresponds to the emission independent of 


momentum (according to Eq. ( 13.12 )), and the lower one (b) corresponds to p-dependent emission 
(according to Eq. ( |13.13| )). One can notice that the resulting pion distributions are non-equilibrium 
ones. We can also see that they coincide for large values of momenta. The solid line represents the 
thermal distribution fitted to these two curves in the high momentum region. One can see that 
for small momenta, 25 MeV < p < 200 MeV, both pion distributions exceed the thermal one. In 
the case (a) the pion distribution is very much peaked in the limit p —> 0. This behaviour reflects 
directly the criticality of the hadronization cross section. On the other hand, in the case (b) the 
pion distribution function drops to zero for p —> 0 since the effect of critical hadronization is, in 
this situation, partially destroyed by the slow emission of low energetic pions. 


In consequence, our model calculation indicates that the criticality in the hadronization cross 
sections leads to the excess in the production of low energetic pions. The effect becomes particularly 
strong if the produced pions, independently of their momenta, can manage to decouple fast from 
the interacting system. The question whether the discussed scenario is really the realistic one and 
responsible for the observed low-py enhancement can only be checked by more detailed calculations. 
They should include, in particular, the dynamics of the space-time evolution of the system and of 
the freeze-out process. As we showed in our simple approach it is, in general, very important to 
take into account the possibility that the cross sections exhibit non-trivial energy dependence. This 
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can lead to interesting phenomena connected with the existence of the phase transition and even 
serve as a diagnose of its occurrence. 



Figure 13.1: The distribution functions of the emitted pions (dashed lines) together with the thermal 
distribution (solid line) fitted to them in the hig h mom entu m regio n. The cases (a) and (b) correspond to 
the emission times calculated according to Eqs. ( 13.12|) and (13.13), respectively. 
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Part IV 

SUMMARY 
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1 In general the vacuum correlation functions in the meson channels are screened and, due 
to Lorentz invariance, the screening mass of a stable meson equals its dynamic mass. For 
unstable mesons the screening mass equals the mass of the branch point, since the asymptotic 
form of the correlation function picks out the lowest lying singularity. As demonstrated by 
our work, these results are reproduced by the NJL model: the pion screening mass is equal 
to its dynamic mass, whereas the sigma screening mass is 2Mo, where Mq is the constituent 
quark mass in vacuum. 

At finite temperature and zero density, the correlation functions are again exponentially 
damped. Moreover the screening and dynamic masses differ — this fact indicates that at 
finite temperature it is impossible to obtain the exact information about time-like excitations 
by studying space-like correlation functions. The overall temperature dependence of the 
screening masses in the NJL model agrees with the corresponding lattice results. 

We carefully explored the analytic structure of the correlation functions. In the complex 
momentum plane they have two cuts (scalar channel) or two cuts and a pole (pseudoscalar 
channel). At sufficiently high temperature the pion pole disappears and we have to deal, in 
both channels, with two cuts. The contributions from these cuts show large cancellations 
which lead to the exponential decay of the correlation functions. 

2 We have demonstrated that the correlation functions at T = 0 and finite density differ qual¬ 
itatively from those in vacuum and those at T > 0: they exhibit long ranged oscillations, 
of the Friedel type, rather than exponential damping. We found that the appearance of the 
oscillations is connected with the existence of a cut of finite range in the complex momen¬ 
tum plane. This cut is responsible for the leading contribution to the correlation function 
at large distances. The length of the cut, which is proportional to the Fermi momentum of 
the constituent quarks, is reflected in the oscillation period at large distances 5r = ttIpf- 
Consequently this form of the correlation function is quite general and is expected in all nor¬ 
mal Fermi liquids. In particular, the existence of the oscillations is independent of whether 
the basic fermionic degrees of freedom are quarks or nucleons. Therefore, we feel the oscil¬ 
latory behaviour of the correlation function at finite density will not change qualitatively by 
confinement. 

The fact that the correlation function exhibits Friedel-type oscillations indicates again that 
in medium it is impossible to obtain the information about time-like excitations by studying 
space-like correlations. Their long distance behaviour is dominated by low lying particle-hole 
excitations, and consequently not connected with the dynamic mass. 

3 Our NJL studies concerning static meson correlation functions were supplemented by the 
calculation within the perturbative QCD. Keeping only the lowest order term we obtained 
compact analytic results. This work extends some of the earlier results, which were concen¬ 
trated only on the asymptotics of the correlation functions. 

4 Temperature dependence of the quark condensate was studied in the NJL model with meson 
loops. Substantial differences were found compared to the standard results obtained in the 
Hartree-Fock approximation. In particular, we found finite slope of the condensate vs. 

for T ^ 0 in the chiral limit, faster melting of the condensate, and a lower chiral restoration 
temperature. The importance of pions is connected with the smallness of their mass, much 
smaller than the masses of other hadrons. In the particular case of the NJL model, the pion 
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mass is much smaller than the mass of the non-confined constituent quarks. The behaviour of 
l^q) is affected by the coupling of the quarks with pions. This happens in such a way that the 
presence of pions leads to the destruction of the condensate. Therefore, including meson loops 
we speed up the decrease of the condensate with increasing temperature. We demonstrated 
consistency of our results with the chiral perturbation theory. Our calculation represents the 
first self-consistent treatment of the NJL model with meson loops at finite temperature. 

5 The chiral symmetry concepts have been explicitly included in the construction of the trans¬ 
port equations for quark matter. Our starting point was the chirally invariant Lagrangian 
of the NJL model and we derived the transport equations via a spinor decomposition of the 
Wigner function and a gradient expansion. Our calculation was restricted to the mean-field 
approximation. We have taken into account the spin dynamics and discussed the possibility 
of having a non-zero pseudoscalar condensate. In this aspect, our results are a generalization 
of some earlier results. 

The classical quark distribution functions satisfy the kinetic equations of the standard form 
with the effective chirally invariant mass M‘^{X) = + (T^g^(X), where 7r(o)(J^) and 

(T(o)(^) leading terms in the classical expansion of the pseudoscalar and scalar con¬ 
densates. Furthermore, the angle ‘L(X), defined through the relation 7r(o)(X)/cj(o)(^) = 
tan<L(X), must be a constant. However, its value remains undefined, which reflects the chiral 
symmetry of the problem. 

The classical equation for the spin evolution, which has been derived for the first time by us, 
is also invariant under chiral transformations. However, its solutions are constrained by the 
additional condition of the axial current conservation, which is not simple to incorporate. 

The inclusion of the quark mass term into the Lagrangian explicitly breaks the chiral invari¬ 
ance. Consequently, we find that <h(X) as well as 7r(o)(J^) must be zero. Nevertheless, the 
general form of the kinetic equations does not change and we can still use them with the 
substitution M{X) = (T(o)(Jf) + w-, where m is the current quark mass. Moreover, in this 
case the requirement of the axial current conservation reduces to a form familiar from PC AC. 

6 The chirally invariant kinetic theory for quark matter was used to study large time-scale 
fluctuations of the quark condensate. Our results showed that the general features of the 
propagation of the excitations in the system were very similar to those characterizing the 
wave propagation in a non-relativistic Maxwellian plasma. The elementary fluctuations of the 
condensate (Van Kampen modes) are not damped; they correspond to initial perturbations of 
the quark distribution function, which are strongly peaked for some given value of momentum. 
On the other hand, the superpositions of such elementary fluctuations are in practice always 
damped due to a phenomenon of destructive coherence (Landau damping). The overall 
results of this study showed the stability of the in-medium gap equation with respect to small 
perturbations of the quark distribution functions. 

7 The system of kinetic equations describing hadronization of the quark-antiquark plasma was 
analyzed. It was found the the singularities in the elastic and hadronization cross sections, 
as described by the NJL model in connection with the critical behaviour at a second order 
phase transition, favor production of low-energetic mesons. We found that this effect can give 
a contribution to the observed low-p^ enhancement of pions observed in the ultra-relativistic 
heavy-ion collisions. 
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